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Preface 


As  a  flight  control  engineer  in  the  Flight  Dynamics 
Laboratory  of  the  Air  Force  Wright  Aeronautical  Labora¬ 
tories,  I  had  the  opportunity  to  work  with  advanced  air¬ 
craft  control  objectives  and  designs.  I  was  impressed 
with  the  difficulty  of  achieving  ultimate  aircraft  control 
behavior  when  employing  conventional  design  techniques 
derived  from  designs  for  single-input  single-output  (SISO) 
systems.  There  is  a  substantial  gulf  between  techniques 
and  design  methodologies  of  the  so-called  modem  control 
theory  and  the  designs  that  are  actually  pursued  by  prac¬ 
ticing  flight  control  engineers.  This  gulf  is  becoming 
still  broader  in  the  realm  of  controller  design  for 
sampled-data  implementation. 

I  have  held  an  interest  for  several  years  in 
optimal  model -following  design  as  a  method  for  achieving 
advanced  flight  control  systems  for  sampled-data  implementa¬ 
tions.  The  Command  Generator  Tracker  control  system  con¬ 
sidered  in  this  study  is  a  new  development  of  the  model¬ 
following  design  method,  and  offers  several  advantages  over 
earlier  such  designs. 

I  wish  to  thank  ray  thesis  advisor,  Dr.  Peter  S. 

Maybeck  of  the  Air  Force  Institute  of  Technology,  for  his 
interest  in  this  study.  His  consistent  concern  for 
thorough  and  accurate  research  and  reporting  is  the  model 
for  my  own  efforts  in  this  and  future  studies  I  may  undertake. 
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Abstract 


This  study  develops  a  computer  program  for  inter¬ 
active  execution  to  aid  in  the  design  of  Command  Generator 
Tracker  control  systems  employing  Proportional-plus-integral 
inner-loop  controllers  and  Kalman  Filters  for  state  estima¬ 
tion  (CGT/PI/KF  controllers) .  Design  parameters  are  speci¬ 
fied  in  the  continuous-time  domain  and  the  computer  program 
obtains  the  corresponding  discrete-time  parameters  and  deter¬ 
mines  a  direct  digital  design  for  sampled-data  implementa¬ 
tion.  Designs  are  based  upon  the  Linear  system  model, 
Quadratic  cost,  and  Gaussian  noise  process  (LQG)  assumptions 
of  optimal  control  theory. 

The  report  discusses  the  theoretical  background  and 
applications  of  optimal  model -following  designs  which  pre¬ 
ceded  the  CGT  theory.  A  development  of  the  CGT/PI/KF  con¬ 
troller  theory  is  presented,  and  performance  evaluation 
tools  for  the  controller  design  are  discussed.  Following  a 
brief  description  of  the  computer  program  developed,  results 
of  applying  it  to  example  aircraft-related  controller  design 
problems  are  presented  and  discussed.  ^ Among  the  designs 
presented  are  controllers  for  conventional  pitch  rate  and 
decoupled  pitch-pointing  control  for  an  Aircraft  system 
model  representative  of  modern  aircraft  longitudinal 
dynamics.  The  CGT/PI/KF  controller  is  found  to  be  a 


x 


technique  particularly  well  suited  to  the  typical  aircraft 
control  design  problem  wherein  a  multi-input  multi-output 
(MIMO)  system  is  to  have  specified  output  response  behavior 
to  commanded  inputs  ("handling  qualities")  while  simul¬ 
taneously  rejecting  disturbances  of  specifiable  characteris¬ 
tics. 

The  computer  program  is  fully  documented  in  the 
appendices  of  the  report.  Included  are  a  "Programmer's 
Manual,"  a  "User's  Manual,"  sample  program  input  and  output 
a  program  listing,  and  a  listing  of  job  control  language 
required  to  obtain  an  executable  object  file.  These  per¬ 
tain  to  the  computer  program  as  implemented  on  a  Control 
Data  Corporation  CYBER  computer  system  and  interactive 
execution  under  INTERCOM. 
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DESIGN  OF  ADVANCED  DIGITAL  FLIGHT  CONTROL  SYSTEMS 
VIA  COMMAND  GENERATOR  TRACKER  (CGT) 

SYNTHESIS  METHODS 

I.  Introduction 

1.1  Background 

Modern  aircraft  designs  entail  increasingly  strin¬ 
gent  and  complex  control  requirements.  Newer  aircraft 
employ  digital  flight  control  systems,  utilize  multiple 
control  surfaces  in  each  axis,  and  demand  highly  refined 
control  characteristics — both  in  coupled  and  decoupled  com¬ 
mand  modes. 

Early  flight  control  designs  initially  served 
pilot  relief  functions  exclusively.  Later  control  systems 
were  designed  to  improve  aircraft  stability  and  control 
characteristics,  but  had  limited  control  authority.  Such 
systems  were  designed  in  an  evolutionary  fashion,  with  con¬ 
siderable  trial  and  error,  and  generally  dealt  with  multiple 
control  surfaces  in  an  ad  hoc  fashion. 

While  modern  control  design  techniques,  as  exempli¬ 
fied  by  optimal  control  theory,  showed  considerable  promise 
in  application  to  flight  control  problems  during  the  early 
and  mid-1960s,  they  have  not  been  adopted  by  the  aircraft 
designers.  The  optimal  control  techniques  have  been  seen 
as  suffering  from  among  the  following  deficiencies: 
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traditional  design  criteria  are  not  readily  specified 
directly  in  the  performance  index;  it  is  difficult  to 
select  appropriate  weighting  matrices  for  the  cost  func¬ 
tionals  to  achieve  the  desired  response;  the  resulting  con¬ 
trollers  require  full-state  feedback  but  measurements  of 
all  states  are  generally  not  available,  so  the  designs 
must  be  reduced  or  approximated  for  implementation  (or 
filters  or  observers  must  be  added) ;  and  the  typical  formu¬ 
lation  of  the  optimal  controller  solves  only  the  regulator 
problem  and  not  the  required  tracking  response  to  a  command 
input. 

To  achieve  greater  utility  in  the  application  of 
optimal  control  techniques  to  flight  control  problems,  a 
control  synthesis  technique  known  as  "model -following"  has 
been  used.  Essentially,  the  goal  of  model -following  is  to 
control  a  given  system  so  that  its  outputs  "follow"  those 
generated  by  a  "model"  system  which  represents  the  desired 
dynamic  behavior.  By  formulating  the  optimization  to 
achieve  model -following ,  the  difficulties  associated  with 
defining  an  appropriate  performance  index  and  selecting 
quadratic  weighting  matrices  have  been  alleviated.  But 
other  problems  have  persisted:  the  resulting  designs  still 
require  full-state  feedback;  it  is  difficult  to  achieve 
desired  response  behavior  to  disturbances  affecting  the 
system;  and  the  model-following  achieved  is  only  for  system 
response  to  initial  conditions  with  no  inputs.  Extensions 
of  the  theory  to  allow  model -following  for  forced  response 
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to  command  inputs  have  been  developed  but  require  that  the 
input  itself  be  specified  by  a  dynamic  model. 

Recent  developments  in  modern  control  theory  have 
been  unified  into  a  new  synthesis  technique  for  model- 
following  referred  to  as  Command  Generator  Tracking.  The 
controller  is  designed  so  that  the  system  outputs  follow 
the  trajectories  prescribed  by  a  command  generator,  while 
simultaneously  rejecting  disturbances  with  specified  charac¬ 
teristics.  The  model  generating  commands  as  a  prefilter 
to  the  pilot  inputs  may  incorporate  the  desired  closed- 
loop  dynamics  and  tracking  characteristics,  and  the  con¬ 
troller  can  take  full  advantage  of  all  available  control 
surfaces  appropriate  to  the  control  task  with  systemati¬ 
cally  determined  crossfeed  gains  as  well  as  single-channel 
gains. 

As  discussed  in  this  thesis,  the  Command  Generator 
Tracking  control  system  is  designed  as  a  digital  controller 
composed  of  three  elements:  a  Command  Generator  Tracker 
(CGT)  processing  command  inputs  to  define  inputs  to  the 
system;  a  Proportional-plus-integral  (PI)  regulator  acting 
as  an  inner-loop  controller  to  drive  the  system  to  follow 
the  CGT  inputs;  and  a  Kalman  Filter  (KF)  providing  esti¬ 
mates  of  the  system  and  disturbance  states  needed  by  the 
controllers.  The  availability  of  digital  computers  of  a 
size,  speed,  capacity,  and  reliability  appropriate  to 
flight  control  tasks  makes  digital  controller  designs 
feasible.  The  CGT/PI/KF  controller  is  a  direct  digital 
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design  from  continuous- time  specifications,  especially 
suited  to  the  modem  digital  flight  control  capability. 

The  design  technique  for  the  CGT/PI/KF  controller 
retains  the  desirable  qualities  of  the  earlier  model¬ 
following  techniques.  Furthermore,  it  provides  directly 
for  prescribed  response  to  disturbances,  does  not  require 
that  the  command  input  have  specified  dynamics,  and 
inherently  incorporates  the  state  estimation  needed  to 
implement  the  control  in  the  face  of  only  incomplete  and 
noise-corrupted  measurements  being  available  instead  of 
all  states. 

1.2  Problem 

The  primary  objectives  of  this  thesis  are: 

1.  To  develop  an  interactive,  user-oriented  com¬ 
puter  program  to  aid  in  the  design  of  CGT/PI/KF  controllers. 
It  is  to  be  applicable  to  arbitrary  systems  of  varying 
dimensions.  Systems  and  design  parameters  are  specified 
for  a  continuous- time  problem  representation  but  the  con¬ 
troller  is  a  direct  digital  design. 

2.  To  apply  the  design  program  to  an  aircraft 
flight  control  design  problem  in  order  to  evaluate  charac¬ 
teristics  of  the  CGT/PI/KF  design  technique  and  qualities 
of  the  resulting  controller  designs. 

1.3  Sequence  of  Presentation 

The  results  of  this  effort  are  fully  documented  in 
the  body  of  the  thesis  and  in  the  appendices.  The  thesis 
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and  its  appendices ,  while  integral ,  are  intended  to  serve 
as  separable  entities.  The  main  chapters  of  the  thesis 
consider  the  theoretical  aspects  of  the  CGT/PI/KF  design 
technique  and  demonstrate  practical  application  of  the 
technique  using  the  computer  program  developed  in  this 
thesis  effort.  The  appendices  specifically  document  the 
design  computer  program  to  allow  understanding  of  its  code 
and  operation,  as  well  as  its  successful  application  to 
design  problems  different  from  those  considered  in  the 
thesis. 

The  body  of  the  thesis  report,  composed  of  Chapters 
II  to  IV,  discusses  the  background  of  model-following  con¬ 
trol  designs,  the  theoretical  development  of  the  CGT/PI/KF 
controller,  and  the  evaluation  tools  employed  for  the 
design.  These  chapters  are  followed  by  Chapter  V  which 
presents  a  general  description  of  the  computer  program 
which  was  developed,  and  Chapter  VI  discussing  the 
CGT/PI/KF  controller  designs  achieved  for  several  design 
problems  through  use  of  the  program.  A  final  chapter 
offers  conclusions  and  recommendations  for  further  research. 

The  first  two  appendices  are  guides  to  the  under¬ 
standing  and  use  of  the  program:  the  first  is  a  detailed 
description  of  the  computations  performed  by  the  program 
and  the  computer  source  code;  the  second  describes  the  opera¬ 
tion  of  the  program  from  a  user's  perspective,  including 
discussion  of  specific  items  of  input/output.  These  are 
followed  by  Appendix  C  which  illustrates  sample  input  and 
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output  from  an  execution  of  the  program,  Appendix  D  which 
contains  a  full  source  listing  of  the  computer  code,  and  a 
final  appendix  which  shows  the  job  control  language  needed 
to  obtain  an  executable  program  file. 


II.  Model-Following  Control 


2.1  Introduction 

The  control  philosophy  employed  in  the  CGT/PI/KF 
controller  is  not  new.  Since  the  early  1960s  there  has 
been  work  on  the  class  of  optimal  controllers  referred  to 
as  "model-following  systems."  Early  work  by  Kalman  was 
unified  by  Tyler  in  1964  (Ref  34)  and  presented  as  a 
design  method  appropriate  to  aircraft  control  problems. 
From  1964  through  about  1977,  various  articles  in  the  tech 
nical  journals  developed  extensions  to  the  methods  dis¬ 
cussed  by  Tyler  and  likewise  applied  them  to  aircraft  con¬ 
trol  design.  The  CGT/PI/KF  design  method  was  reported 
first  in  1978  (Ref  8) .  While  it  is  clearly  in  the  class 
of  model -following  controllers  and  has  characteristics 
closely  related  to  those  typical  of  the  earlier  model¬ 
following  designs,  the  theory  from  which  it  is  derived  is 
distinctly  not  in  a  lineal  path  with  the  work  of  the  early 
1960s.  In  one  consistent  development  it  incorporates  all 
the  capabilities  of  the  various  model-following  designs  of 
the  1964-1977  time  frame,  provides  new  capabilities,  and 
does  so  in  a  single  unified  controller/ filter  structure. 
Before  presenting  the  theoretical  development  of  the 
CGT/PI/KF  controller,  it  is  appropriate  to  outline  the 
theoretical  bases  of  the  earlier  model- following  designs 
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briefly  and  to  discuss  some  of  the  various  extensions  and 
applications  which  have  occurred. 

2.2  Model-Following  Control 

in  Aircraft  Design 

The  design  objective  in  employing  model-following 
control  is  to  achieve  a  control  design  which  forces  the 
output  behavior  of  a  given  system  to  be  like  that  of  a 
"model"  system.  This  approach  to  design  of  optimal  con¬ 
trollers  has  been  previously  pursued  in  aircraft  control 
applications  to  design  controllers  yielding  either 
"desirable"  characteristics  or  characteristics  like  those 
of  another  aircraft  (e.g.,  in-flight  simulators  such  as 
the  Calspan  TIFS  aircraft)  (Ref  1) . 

It  has  been  common  to  consider  aircraft  dynamic 
motion  as  decoupled  between  the  longitudinal  and  the 
lateral  axes.  In  each  axis  the  dominant  dynamic  modes  of 
motion  have  been  characterized  by  second-order  response 
models — the  "short-period"  and  "dutch-roll"  modes  of  the 
longitudinal  and  lateral  axes,  respectively.  Correspond¬ 
ing  to  these  characterizations,  it  has  been  common  to 
develop  standards  of  desirable  aircraft  control  behavior 
according  to  specifications  on  a  time-response  modeled  as 
second-order  (Ref  12) .  With  the  advent  of  so-called 
"decoupled"  modes  of  motion,  it  has  also  become  common  to 
define  the  desired  decoupled  responses  in  terms  of  first- 
order  characteristics.  Thus,  while  the  aircraft  itself  may 
entail  dynamics  adequately  depicted  by  a  model  of  relatively 
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high  order,  design  specifications  and  criteria  are  framed 
in  terms  of  low-order  models.  In  parallel  with  this  use  of 
low-order  models  for  control  synthesis,  it  has  been  common 
to  develop  the  specifications  and  "handling  criteria"  for 
aircraft  control  from  experiments  in  which  a  given  aircraft 
is  constrained  to  exhibit  dynamics  either  like  those  of 
another  aircraft  or  like  those  to  be  examined  for  goodness 
or  badness  of  control  quality  (Refs  34,  35). 

Thus  the  model -following  techniques  find  natural 
application  in  the  control  problems  typical  of  aircraft. 

For  control  synthesis,  the  designer  seeks  to  achieve  a  con¬ 
trolled  response  like  that  of  a  first-  or  second-order  sys¬ 
tem  with  specific  attributes.  For  evaluation  and  develop¬ 
ment  of  controllers  and  control  specifications,  the  experi¬ 
menter  seeks  behavior  of  a  type  postulated  as  different 
from  that  inherent  for  the  test  aircraft. 

2.3  Types  of  Model-Following 

Controllers  (Refs  26,  34) 

Two  different  techniques  for  achieving  and  imple¬ 
menting  model- following  controllers  were  developed  in  the 
early  1960s  and  were  presented  formally  by  Tyler  (Ref  34) 
in  1964.  One  technique  is  referred  to  as  a  "model- in- the- 
performance- index"  controller,  or  more  simply  as  an 
"implicit"  model- following  controller.  The  other  is 
referred  to  as  a  "model-in-the-system"  controller,  or  more 
simply  (and  in  contrast  to  the  first  type)  as  an  "explicit" 
model-following  controller. 
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For  the  implicit  scheme,  the  model  is  employed 
only  in  design  to  determine  feedback  gains  for  the  system 
states.  For  the  explicit  scheme,  the  model  is  used  not 
only  to  determine  necessary  feedback  gains  but  also  feed¬ 
forward  gains  on  the  model  states  themselves,  making  it 
necessary  to  include  a  simulation  of  the  model  dynamics 
in  the  controller.  Although  it  is  derived  differently,  the 
CGT/PI/KF  controller  uses  an  explicit  type  of  implementa¬ 
tion. 

As  originally  developed,  both  the  implicit  and 
explicit  model-following  controllers  used  models  with  state 
dynamics  but  no  inputs  (Refs  23,  34).  The  responses  which 
were  matched  for  system  and  model  thus  were  the  responses 
to  initial  conditions.  Later  work  sought  to  include 
response  to  inputs,  but  this  generally  then  further 
required  specification  of  the  input  for  which  matching  was 
desired  (Refs  1,  20,  26,  27). 

2.3.1  Implicit  Model-Following  Controller 
(Refs  23,  26,  34).  In  implicit  model-following,  the  model¬ 
following  objective  is  pursued  by  employing  the  feedback 
gains  of  the  optimal  controller  to  modify  the  coefficients 
of  the  open-loop  system  matrix  so  that  they  approach  those 
of  the  model.  This  is  achieved  by  including  the  model  in 
a  performance  index  of  a  form  suggested  by  Kalman  as  an 
alternative  to  the  usual  index  that  weights  system  state 
or  output  deviations  only  (Ref  22) . 
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For  the  system 


*  «=  Ax  +  Bu  (2-1) 

with  output 

y  =  Cx  (2-2) 

A  standard  performance  index  is  of  the  form 


J 


uTRu] dt 


(2-3) 


which  weights  state  and  input  deviations  and  for  which  the 
matrix  Q  may  be  due  to  weights  directly  on  the  states  or 
derived  from  weightings  on  the  outputs  and  given  by 

Q  =  CTQyC  (2-4) 


The  infinite  terminal  time  plus  time-invariance  of  the  sys¬ 
tem  and  cost  weighting  matrices  are  used  to  generate  a 
constant-gain  steady-state  controller. 

Instead  of  the  index  of  equation  (2-3),  define  a 
performance  index  which  weights  the  error  between  system 
output  derivatives  and  the  model  dynamics  where  the  model  is 


*  =  Ax 

— m  — m — m 

and  the  corresponding  performance  index  is 


(2-5) 


/OO 

I  (^-Amy)TQI(^-Amy)-mTRu]dt  (2-6) 
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and  the  dimension  of  the  output  vector  y  and  the  model 
state  vector  are  the  same  and  the  weighting  matrix 
weights  errors  between  the  output  and  model  dynamics. 
This  index  can  be  rewritten  using  equations  (2-1)  and 
(2-2)  as 

/oo 

[ (CAx+CBu-A  Cx) TQ  (CAx+CBu-A  Cx) 

-  -  — m —  —I  - - m — 

*0 

+uTRu]dt  (2-7) 

Carrying  through  the  various  matrix  multiplications 
and  collecting  terms  leads  to  a  performance  index  similar 
to  that  of  equation  (2-3)  but  with  a  cross-weighting  term 
relating  deviations  in  x  and  u  (Ref  26) : 


^  /  TV  T  ^ 

J!  •  J  [x  Q-J.X+2U  Sx+uARu]  dt 


(2-8) 


where 

=  (CA-A  C)TQ  (CA-A  C) 
I  —  —in—  —I  —  — m— 

S  =  BTCTQt(CA-A  C) 

— - *T - m- 

and 


( 2-9a) 
(2-9b) 


(2-9c) 


Such  an  index  is  appealing  since,  for  a  definition 
of  quadratic  weights  and  R  with  simple  and  direct  mean¬ 
ing  to  the  designer  in  terms  of  the  desired  aircraft 
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response,  potentially  complex  and  not  particularly  obvious 
cost  weightings  result  to  which  the  usual  optimization  tech¬ 
niques  may  be  applied. 


This  directness  of  specification  of  weighting 
matrices  appropriate  to  the  design  objectives,  along  with 
its  inherent  simplicity  in  implementation,  has  made  the 
implicit  model-following  technique  attractive  to  aircraft 
control  system  designers  (Refs  20,  26,  27,  34) .  By  model¬ 
ing  either  specific  inputs  or  classes  of  inputs,  it  has 
been  possible  to  extend  the  technique  to  provide  response 
matching  for  systems  driven  by  inputs  (Refs  20,  26,  27). 

2.3.2  Explicit  Model -Following  Controller  (Refs 
1,  20,  26,  27) .  In  explicit  model-following,  the  model¬ 
following  is  pursued  by  employing  feedback  gains  around  the 
system  to  make  it  behave  as  a  "tight"  tracker,  with  feed¬ 
forward  gains  on  the  model's  states  providing  the  reference 
input.  The  design  is  essentially  of  an  optimal  tracker 
with  the  reference  being  the  output  of  the  linear  system 
defined  by  the  model.  The  controller  gains  are  determined 
from  optimization  employing  a  performance  index  which 
weights  the  differences  between  the  system  and  model  out¬ 
puts. 

For  the  system  defined  by  equations  (2-1)  and  (2-2) 
and  the  model  of  equation  (2-5),  an  appropriate  performance 
index  is  (Ref  26) 
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(2-10) 


£ 


t  <2-^1 


'2Eli-!!ml+!l  sal  at 


with  the  system  output  vector  y  and  the  model  state  vector 
xm  having  the  same  dimension  and  the  weighting  matrix 
weighting  the  differences  between  the  system  outputs  and 
the  model  states.  This  is  a  useful  special  case  of  match¬ 
ing  model  outputs  (here  the  entire  model  state  is  con¬ 
sidered  the  output) . 

To  modify  the  performance  index  to  achieve  the  form 
of  equation  (2-3),  define  an  augmented  system 


with 


*  A  X 
—a— a 


+  B  u 
— a— 


and 


A  = 
—a 


A  [  0 


0  I  Am 

“  i  -1Bj 


(2-11) 


(2-12a) 


(2-12b) 


( 2-12c) 


For  this  augmented  system  the  corresponding  perform¬ 
ance  index  is 


00 

JE  =  /  a^E-a+— T— ^ dt 


(2-13) 
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with 


(2-14) 


The  optimal  control  input  for  the  augmented  system 

then  is 

u*  =  -G£xa  (2-15) 

where  G£  represents  the  optimal  feedback  gain  matrix.  This 
can  be  rewritten  in  terms  of  the  original  system  and  model 
states: 


u*  =  -[G*  x+G*  x  ] 


(2-16) 


where  G*  and  G*  are,  respectively,  the  feedback  gains  on 
"^1  2 

the  system  states  and  the  feedforward  gains  on  the  model 

states.  Tyler  (Ref  34)  demonstrates  that  the  gains  G*  are 

1 

independent  of  the  model  to  be  followed  while  the  G*  gains 

_c2 

depend  both  on  the  system  and  the  model . 

Because  the  quality  of  the  model-following  achieved 
depends  on  the  tightness  of  the  inner-loop  tracker,  the 
explicit  model- following  design  often  entails  rather 
high  feedback  gains.  These  high  feedback  gains  along  with 
the  greater  complexity  due  to  incorporation  of  the  model 
within  the  control  system  has  made  the  explicit  scheme 
less  popular  for  controller  synthesis  than  the  implicit 


scheme,  although  it  has  found  application  in  the  cases  in 
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which  an  experimental  aircraft  is  to  behave  with  specific 
dynamic  characteristics  (Refs  34,  35). 

2.4  Comparison  of  Implicit  and 

Explicit  Model -Following 

For  perfect  model  following,  i.e.  system  output 
and  model  output  being  equal  for  all  time,  Erzberger  (Ref  17) 
showed  in  1968  that  the  implicit  and  explicit  model - 
following  controllers  are  of  equivalent  capability  if  the 
system  is  perfectly  modeled  and  no  disturbances  impinge 
upon  it  in  actual  operation.  For  either  controller,  the 
following  must  be  true  (Ref  17) : 

[ (CB) (CB) +-i]  [A^C-CA]  =0  (2-17) 

where  the  superscript  t  denotes  the  natrix  pseudo- inverse , 
and  where  A,  B,  C,  and  Am  are  as  defined  in  equations 
(2-1),  (2-2),  and  (2-5).  This  relation  was  derived  from 
algebraic  conditions  ensuing  from  the  requirement  that 
system  and  model  outputs  be  exactly  equal.  For  implicit 
model-following  this  exact  equality  leads  to  a  condition 
on  the  range  spaces  of  (CB)  and  [A^C-CA] .  For  explicit 
model-following  the  requirement  of  equation  (2-17)  follows 
from  constraining  the  time  derivatives  of  the  system  and 
model  outputs  to  be  equal  at  the  initial  time,  and  the  con¬ 
dition  that  the  control  input  must  be  bounded.  If  equation 
(2-17)  is  not  satisfied,  then  neither  model- follower  can 
achieve  perfect  tracking  with  bounded  inputs  (Ref  17) ,  and 
increases  in  the  quadratic  weights  Q^.  or  Q£  of  equations 
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(2-6)  or  (2-10)  will  not  improve  the  model-following 
beyond  some  minimum  for  the  system-model  pair. 

But  in  the  general  case,  and  particularly  when  the 
system  is  not  perfectly  modeled  and/or  disturbances  act  on 
it,  the  implicit  and  explicit  model-followers  have  differ¬ 
ent  characteristics.  Following  is  a  summary  of  some  of 
those  characteristics  discussed  in  the  works  of  the  1964 
to  1977  era  and  referenced  previously  in  this  chapter. 

By  its  very  nature  of  weighting  deviations  on 
rates,  the  implicit  scheme  places  primary  emphasis  on  the 
controlled  system's  transient  behavior.  It  does  not 
guarantee  matching  of  steady-state  behavior  even  for  the 
nominal  parameter  values.  The  feedback  gains  on  the  system 
tend  to  be  lower  than  for  the  explicit  scheme,  but  their 
values  and  the  ultimate  success  achievable  for  the  design 
depend  on  the  initial  disparity  in  the  open-loop  system 
matrix  and  the  model  matrix  (A  and  Am,  respectively) . 

The  relatively  low  gains  and  the  inherent  sensitivity  of 
those  gains  to  the  specific  values  of  the  matrix  coeffi¬ 
cients  for  the  system  make  the  implicit  controller  more 
sensitive  to  model  inaccuracies  and  parameter  variation 
than  the  explicit  controller.  Also,  since  the  model  out¬ 
puts  and  sys  :em  outputs  are  not  actually  compared  in  the 
controller  implementation,  the  controller  is  blind  to  any 
errors  between  them  and  can  take  no  corrective  action  to 
align  them.  However,  in  contrast  to  the  explicit  scheme, 
since  the  system  dynamics  are  modified  to  match  those  of 
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the  model,  response  of  the  closed-loop  system  to  random 
zero-mean  disturbances  (such  as  clear-air  turbulence)  can 
be  made  to  have  desirable  characteristics  without  the  need 
for  disturbance  states  in  the  controller.  Finally,  the 
most  significant  advantage  of  the  implicit  model -follower 
is  its  simplicity  in  implementation,  and  it  is  primarily 
this  characteristic  which  motivated  earlier  strong  interest 
in  it  as  a  control  synthesis  technique. 

Due  to  the  weighting  of  the  difference  between 
system  outputs  and  model  outputs,  the  explicit  scheme 
places  relatively  greater  emphasis  on  the  steady-state 
behavior  than  on  the  transient  behavior  of  the  model¬ 
following.  Thus,  for  nominal  parameter  values,  a  perfect 
system  model,  and  no  disturbances  acting  on  the  actual 
system,  the  explicit  controller  will  achieve  model¬ 
following  in  the  steady-state.  But  for  gains  of  magnitude 
comparable  to  those  of  the  implicit  controller,  the  quality 
of  the  model-following  in  the  transient  phase  of  response 
will  often  be  worse  for  the  explicit  controller.  Since 
the  inner-loop  feedback  must  achieve  tight  tracking  perform¬ 
ance,  the  gains  required  generally  are  greater  than  those 
needed  for  the  implicit  scheme.  On  the  other  hand,  these 
higher  gains  along  with  the  fact  that  the  system  outputs 
and  model  states  are  actually  compared  by  the  controller 
make  the  explicit  controller  less  sensitive  to  parameter 
variation,  modeling  errors,  and  errors  in  initial  condi¬ 
tions.  Model -following  is  achieved  only  with  respect  to 
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states  of  the  model  which  are  driven  through  defined  com¬ 
mand  input  channels.  Thus,  unless  disturbances  are  spe¬ 
cifically  incorporated  into  the  model,  the  closed-loop 
system  response  to  disturbances  will  be  like  that  of  the 
inner-loop  system,  which  may  not  be  what  the  designer 
wishes.  Finally,  since  the  command  inputs  must  first  be 
processed  through  the  model  dynamics  before  becoming  the 
reference  input  to  be  tracked  there  may  be  some  time  delay 
in  the  system's  response  to  inputs. 

2.5  Introduction  to  the 

CGT/P I /KF~Cont roller 

In  the  previous  section  it  was  seen  that  the 
explicit  model-following  controller  structure  potentially 
offers  important  advantages  due  to  its  inherently  superior 
performance  in  real  implemen tat ions,  i.e.  controllers  for 
systems  not  perfectly  known  and  subject  to  parameter 
variation.  Its  primary  disadvantages  are  its  greater  com¬ 
plexity  and  its  potentially  undesirable  performance  when 
the  system  is  subject  to  random  disturbances. 

With  the  readily  available  digital  computers  for 
implementation  of  real-time  digital  control  laws,  the  com¬ 
plexity  of  the  controller  becomes  a  less  significant  issue. 
Accepting  that  a  greater  degree  of  controller  complexity 
can  be  accommodated  in  aircraft  designs  of  the  future,  the 
CGT/PI/KF  controller's  characteristics  make  it  attractive 
for  achieving  design  objectives. 
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Like  the  explicit  model- follower  discussed  in  this 
chapter,  the  CGT/PI/KF  controller  employs  feedforward  gains 
providing  a  reference  input,  but  to  a  proportional-plus- 
integral  inner  loop  regulator  which  seeks  to  maintain  the 
difference  between  the  reference  and  the  system  outputs  at 
zero.  Both  the  feedforward  and  the  feedback  gains  are  com¬ 
puted  independently  without  resort  to  an  augmented  system 
description  and  corresponding  augmented  Riccati  equation. 
The  final  gain  matrices  for  the  closed-loop  controller  are 
then  obtained  from  the  two  independent  solutions  and  char¬ 
acteristics  of  the  system,  using  simple  matrix  multiplica¬ 
tion. 

Unlike  the  explicit  model-follower,  the  CGT/PI/KF 
can  readily  include  models  of  random  disturbances  affecting 
the  system.  The  resulting  controller  then  can  be  made  to 
reject  the  corresponding  real  disturbances,  in  the  stochas¬ 
tic  zero-mean  sense.  Thus,  disturbance  rejection  can  be 
concentrated  in  those  frequency  bands  where  predominant 
disturbances  are  expected  in  actual  use. 

Since  the  inner-loop  of  the  CGT/PI/KF  controller 
consists  of  a  PI  regulator,  the  actual  system  controlled 
response  will  achieve  model-following  in  steady-state 
despite  errors  in  system  modeling,  parameter  variation, 
or  unmodeled  constant  disturbances.  Thus  the  CGT/PI/KF 
controller  is  less  sensitive  to  model  definition  errors 
than  the  earlier  explicit  scheme. 
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As  is  described  in  Chapter  III,  the  determination 
of  the  feedforward  gains  for  the  CGT  controller  makes  a 
single  assumption  for  the  command  input:  the  command  input 
is  assumed  to  vary  slowly  in  comparison  to  the  system  and 
model  dynamics  and  thus  is  approximately  constant  during 
the  controller  sampling  interval.  With  this  minor  restric¬ 
tion  on  the  command  input,  feedforward  gains  on  the  input, 
the  model  states,  and  the  disturbance  states  can  be  deter¬ 
mined.  Moreover,  the  controller  will  employ  same-cycle 
feedthrough  of  the  command  input  both  directly  to  the  sys¬ 
tem  input  and  to  the  model  dynamics  update.  In  this  way 
the  potential  for  delay  in  system  response  to  inputs  is 
greatly  alleviated  and  transient  response  is  enhanced  sig¬ 
nificantly. 

Finally,  since  the  controller  inherently  incorpo¬ 
rates  a  Kalman  filter  for  estimation  of  modeled  system  and 
disturbance  states,  all  variables  needed  for  controller 
implementation  are  available.  Moreover,  account  is  thus 
properly  taken  of  the  true  stochastic  nature  of  the  problem, 
wherein  the  system  and  disturbances  may  be  driven  by  noises 
modelable  as  zero-mean  Gaussian  random  processes,  and  the 
measurements  of  the  system  may  be  both  incomplete  and  cor¬ 
rupted  by  noises  modelable  as  zero-mean  Gaussian  random 
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III.  CGT/PI/KF  Theoretical  Development 

3.1  Overview  of  the  Theory 

The  design  objective  in  employing  the  Command 
Generator  Tracker  control  system  is  to  constrain  a  given 
system  so  that  its  output  response  to  commanded  inputs 
follows  a  model  trajectory  while  rejecting  modeled  dis¬ 
turbances.  Both  the  model  trajectory  and  the  disturbances 
are  derived  as  outputs  of  linear  system  models,  with  the 
command  model  driven  by  the  command  inputs. 

Although  the  CGT  design  does  not  ensure  perfect 
tracking  of  the  model  outputs  during  the  transient  phase 
of  response,  by  formulating  the  design  equations  based  upon 
such  an  idealization,  the  necessary  feedforward  gains  from 
the  command  model  states,  command  inputs,  and  disturbance 
states  to  the  system  inputs  can  be  readily  derived. 

The  CGT  design  solution  can  be  formulated  as  an 
open-loop  design,  depending  only  upon  the  system,  command, 
and  disturbance  models.  However,  if  the  system  is  mar¬ 
ginally  stable  or  unstable,  is  sluggish  or  otherwise  ill- 
behaved  in  response  to  inputs,  then  an  inner-loop  regulator 
may  be  employed.  In  addition,  the  system  model  is  often 
not  known  with  certainty  or  not  modeled  in  full  detail  in 
the  design,  and  unmodeled  disturbances  may  also  impinge 
upon  the  system.  Therefore,  the  preferred  implementation 
consists  of  an  inner-loop  controller  employing  state 
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feedback  to  act  as  a  regulator  with  Proportional-plus- 
integral  control  action  (to  follow  nonzero  commands  with 
zero  steady-state  error) ,  and  with  the  feedforward  gains 
of  the  CGT  providing  the  translation  from  command  inputs 
to  system  inputs. 

Moreover,  if  the  system  and  disturbance  states  are 
not  all  available  as  needed  by  the  controller,  or  if  noise 
corrupts  the  available  measurements,  then  a  Kalman  filter 
may  be  employed  for  system  and  disturbance  state  estimation 
in  the  overall  controller  implementation. 

Thus,  the  controller  design  to  be  developed  in  this 
thesis  consists  of  a  Command  Generator  Tracker  (CGT)  pro¬ 
viding  inputs  to  the  system,  a  regulator  with  Proportional- 
plus-integral  control  action  (PI)  operating  on  those  inputs 
and  the  system  states  so  as  to  drive  the  system  outputs 
along  the  model  output  trajectory,  and  a  Kalman  filter  (KF) 
providing  disturbance  state  estimates  to  the  PI  controller. 
A  general  block  diagram  showing  the  resulting  CGT/PI/KF 
controller  structure  is  shown  in  Figure  3-1. 

While  the  CGT/PI/KF  controller  can  be  developed 
as  a  continuous-time  controller,  it  is  developed  here  in 
its  discrete-time  form.  With  the  increasing  availability 
of  small,  fast,  and  rugged  computers,  many  controller 
designs  are  being  implemented  as  discrete-time  algorithms 
operating  on  sampled  system  measurements  and  providing 
discrete-valued  inputs  to  the  system  through  sample-and- 
hold  devices .  Advantages  of  such  an  implementation  include 
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Fig.  3-1.  CGT/PI/KF  General  Block  Diagram 


greater  operating  reliability,  easier  modification  of  the 
control  laws,  and  the  opportunity  to  select  from  among 
various  different  controllers  each  of  which  may  be  more 
sophisticated  than  feasible  for  a  continuous-time  con¬ 
troller. 

Although  the  controller  design  determines  a 
discrete-time  control  law,  the  design  itself  proceeds  from 
continuous-time  specifications.  The  system  to  be  con¬ 
trolled  is  generally  a  continuous-time  process  and  so  it 
is  appropriately  defined  by  a  continuous- time  model. 
Similarly,  the  desired  system  performance  is  best  repre¬ 
sented  by  a  continuous- time  model.  Also,  most  designers 
have  acquired  their  design  experience  from  problems  which 
were  posed  and  solved  entirely  in  the  continuous-time 
domain.  To  take  advantage  of  this  experience  and  in  recog¬ 
nition  of  the  fact  that  the  design  objectives  will  be 
formulated  as  requirements  on  the  behavior  of  the  system 
for  all  time,  and  not  just  the  sampling  instants,  the  param¬ 
eters  affecting  the  sampled-data  control  law  determination 
are  derived  from  the  corresponding  continuous-time  design 
parameters  provided  by  the  designer. 

The  development  here  proceeds  from  the  relevant 
problem  description  formulated  in  the  continuous- time 
domain  to  solution  in  the  discrete-time  domain,  i.e.,  a 
sampled-data  controller  is  synthesized.  The  elements  of 
this  development  are  derived  from  the  work  presented  in 
References  5,  6,  8,  9,  and  32. 
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3.2 


System  Models 


The  system  for  which  a  CGT/P1/KF  controller  is  to 
be  developed  is  assumed  to  be  well-represented  by  a  set  of 
linear,  time- invariant,  stochastic  state  differential  equa¬ 
tions  with  zero-mean  white  Gaussian  noise  driving  the  sys¬ 
tem  and/or  disturbance  states,  and  corrupting  the  system 
measurements.  Such  equations  generally  are  derived  as 
linearized  perturbation  equations  for  a  non-linear  system 
about  a  nominal  operating  point.  The  assumed  time- 
invariance  in  this  context  is  associated  with  the  derived 
perturbation  model  which  may  be  slowly  varying  but  is 
treated  quasi- statically  for  design.  Thus,  the  system  actu¬ 
ally  is  represented  by  sets  of  such  models,  each  approxi¬ 
mately  valid  near  the  operating  point  at  which  defined,  and 
control  laws  are  designed  for  each.  The  designer  may 
finally  formulate  the  controller  gains  as  functions  of  vari¬ 
ous  parameters  which  serve  to  define  the  various  design 
point  nominal  conditions  and  may  implement  the  control  design 
as  a  controller  of  fixed  structure  but  variable  gains. 

Two  system  models  are  employed:  a  "truth  model" 
and  a  "design  model."  The  truth  model  is  a  model  of  the 
system  which  is  as  complete  and  accurate  as  possible  for 
the  control  task  under  consideration.  Since  the  truth 
model  may  be  of  high  order,  may  include  states  the  designer 
would  prefer  not  to  employ  for  feedback  (such  as  control 
actuator  states) ,  or  may  include  effects  the  designer  con¬ 
siders  relatively  insignificant,  it  is  often  desirable  to 


26 


use  a  different  (generally  simpler)  model  for  design.  This 
design  model  is  the  basis  for  the  controller  and  filter 
gains  and  provides  the  set  of  states  to  which  these  gains 
are  applied. 

While  the  quality  of  the  resulting  designs  is  often 
developed  by  initially  evaluating  their  performance  with 
respect  to  the  design  model,  their  fine-tuning  to  final 
solutions  must  ultimately  be  with  respect  to  the  system 
truth  model.  In  addition,  even  if  the  designer  employs 
the  same  model  states  and  parameters  in  both  the  design 
and  truth  models,  the  truth  model  is  useful.  With  the 
truth  model  available,  the  effects  of  parameter  variation 
may  be  evaluated  since  it  can  be  modified  while  retaining 
the  controller  as  designed  for  the  nominal  parameter  values. 

3.2.1  Design  Model .  The  design  model  consists  of 
a  system  state  differential  equation,  a  disturbance  state 
differential  equation,  an  output  equation,  and  a  measurement 
equation . 

The  system  state  differential  equation  is  given  by 

x  ( t )  =  Ax  ( t)  +  Bj  (t)  +  E  n,  ( t)  +  Gw(t)  (3-1) 

where  the  under-tilde  denotes  that  the  variable  is  modeled 
as  a  random  process,  x  and  n^  are  the  Gaussian  system  and 
disturbance  state  vectors  respectively,  and  w  is  a  zero- 
mean  white  Gaussian  noise  with  covariance 
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(3-2) 


E{w(t)wT(t+T)}  =  Q6(t) 

where  6(0  is  the  delta  function  and  u  is  the  system 
input . 

The  disturbance  state  differential  equation  is 

A^t)  =  A^tt)  +  GnWj(t)  (3-3) 

where  w^  is  a  zero-mean  white  Gaussian  noise  independent 
of  w  and  with  covariance 

E{wd(t)  w^(t+T  )  }  =  Qn<$(T)  (3-4) 

The  output  equation  is 


£(t)  =  Cx  (t)  +  Dyu(t)  +  E^n^t)  (3-5) 

Finally,  the  measurement  equation  is  assumed  to 
provide  discrete-time  measurements  given  by 


z(t.)  -  Hx ( t . )  +  H  n,(t.)  +  v(t.)  (3-6) 

—  i  — —  l  — n-^d  i  —  l 

where  the  measurement  noise  v  is  zero-mean  white  Gaussian 
discrete-time  noise  independent  of  both  w  and  w^,  and  of 
covariance 


E{ v(ti) v^ ( t j ) }  =  R6ij 


with  6 . .  the  Kronecker  delta  function  defined  as 


SO  when  i  t  j 
1  when  i  =  j 


(3-7) 


(3-8) 
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The  dimensionalities  for  the  design  model  are, 
n  =  number  of  system  states 
r  =  number  of  system  inputs 


p  =  number  of  system  outputs 
m  =  number  of  system  measurements 
d  =  number  of  disturbance  states 
w  =  number  of  independent  system  noises 
wD  =  number  of  independent  disturbance  noises 
A  constraint  to  be  imposed  on  the  dimensionalities 
is  that  p=r.  While  a  design  can  be  achieved  in  the  more 
general  case  by  employing  the  matrix  pseudo- inverse,  in 
general  it  is  convenient  to  consider  equal  numbers  of  out¬ 
puts  and  inputs. 

The  dimensions  of  the  various  matrices  can  be 
inferred  from  the  equations  and  the  above  dimensionalities. 
However,  since  explicit  and  clear  knowledge  of  the  sizes 
of  all  matrices  is  needed  for  computer  coding  of  the  equa¬ 
tions  developed  in  this  chapter,  appendices  dealing  spe¬ 
cifically  with  the  design  computer  program  fully  delineate 
the  sizes  of  all  arrays. 

The  design  model  is  discretized  at  a  specific  fixed 
controller/ filter  sampling  period  T  as  follows: 

The  disturbance  state  description  is  augmented  to  the  sys¬ 
tem  state  model  to  form 

x  ( t )  =  A  x  ( t)  +  B  u (t)  +  G  w  ( t)  (3-9) 

—a  — a~a  —a—  —a— a 

where 
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X  (t) 

u 


x(t) 

,  W  (t)  * 

w(t) 

1 — 
V 

ft 

1  w 

~a 

*d(t> 

and,  from  equations  (3-1)  and  3-3) , 


—a  " 


r  i  1 
A  E 

--4--S- 


2  i  A 


■n 


B.  = 


B 

0 


G  = 
—a 


i 


G  j  0 


2  i  G 


L  i  _J 


and,  from  equations  (3-2)  and  (3-4) , 


E{w*(t>Wa(t+T)}  =  Qa6(x) 


where 


2a 


fQ  !  0 


0  1  Qn 

—  i  — n 


(3-9a) 
( 3— 9b) 


(3-10a) 


(3-10b) 


(3-10c) 


(3-10d) 


(3-10e) 


The  corresponding  discrete-time  state  transition  equation 
is 


VW  =  4.5. (ti>  *  ^  u(ti)  +  w,  (ti)  (3-11) 

*v  -  a  a 


where ,  assuming  that  u  is  constant  over  a  sample  period, 


j 
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Although  the  matrix  of  (3-13c)  is  not  generally 
block  diagonal,  this  does  not  impact  the  feedback  control 
that  assumes  full  state  knowledge.  Invoking  the  "Certainty 
Equivalence"  property  discussed  in  Section  3.6,  the  discrete 
time  models  to  be  used  for  the  deterministic  controller 
design  are  the  system  state  transition  equation  defined  by 


x<t.+1>  =  *x(t.)  +  BjUtt.)  +  E.  nj(t.)  (3-14) 

a 


and  the  disturbance  state  transition  equation 


^d(ti+l)  -  -n2d(ti) 


(3-15) 


The  corresponding  output  equation  is 


y(t±)  =  Cx(ti)  +  Dyu(ti)  +  Eynd(ti)  (3-16) 

with  matrices  C,  Dy,  and  Ey  as  in  equation  (3-5) . 

Henceforth,  all  equations  relating  to  the  con¬ 
troller  design  are  considered  deterministic.  As  discussed 
in  Section  3.6,  the  resulting  deterministic  optimal  con¬ 
troller  gains  are  identical  to  the  gains  in  an  LQG  stochas¬ 
tic  optimal  controller  in  which  a  Kalman  filter  provides 
state  estimates  from  incomplete  noise-corrupted  measure¬ 
ments  rather  than  assume  perfect  state  knowledge  (Refs 
2,  32). 


3.2.2  Truth  Model .  For  this  development,  the 
truth  model  consists  of  a  state  differential  equation,  a 
measurement  equation,  and  two  equations  relating  the 


32 


system  and  disturbance  states  of  the  design  model  to  the 
truth  model  states. 

Any  disturbances  which  may  be  considered  to  impinge 
upon  the  system  are  incorporated  into  system  states  for  the 
truth  model.  The  system  state  model  then  becomes 

xt  =  Afcxt(t)  +  Btut(t)  +  Gtwt(t)  (3-17) 

where  xt  is  the  system  state  and  modeled  as  a  Gaussian 
random  process,  wfc  is  a  zero-mean  white  Gaussian  noise 
with  covariance 

E{ wt ( t ) w^ ( t+T ) }  =Qt6(x)  (3-18) 

and  ufc  is  the  system  input. 

The  measurement  equation  is  assumed  to  provide 
discrete-time  measurements  given  by 

-t(ti}  =  -t-t(ti)  +  ^t(ti)  (3-19) 

where  the  measurement  noise  vfc  is  zero-mean  white  Gaussian 
discrete-time  noise  independent  of  with  covariance 

E(vt(ti)v^(tj)}  =  Rt6±j  (3-20) 

with  as  defined  by  equation  (3-8) . 

The  dimensionalities  for  the  truth  model  are, 
n,p  =  number  of  system  states 

rT  =  number  of  system  inputs 

=  number  of  system  measurements 

wT  =  number  of  independent  noises 
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Dimensional  compatibility  for  computations  necessi¬ 


tates  that  the  number  of  measurements  and  inputs  be  equal 
for  both  design  and  truth  models:  m  «  and  r  » 

The  additional  equations  relating  the  system  and 
disturbance  states  of  the  design  model  to  the  system  states 
of  the  truth  model  are,  using  a  prime  to  distinguish  these 
from  the  states  of  equations  (3-1)  and  (3-3) 


^  "  -DTrt 


(3-21a) 


and 


?a  -  ImSt 


(,3-21b) 


The  truth  model  is  discretized  for  a  specific 
controller/filter  sampling  period  T,  yielding 


St'W  +  +  2ta(ti'  (3-22) 

where,  for  u  constant  over  a  sample  period. 


*  -  ^ 
*t  “  a 


(2-23a) 


lt(x)BtdT 


(3-23b) 


and  w  is  zero-mean  white  Gaussian  discrete-time  noise 
-td 

with  covariance 


It(T)^tSt-t^t(T)dT 


(3-23c) 


[ 
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3.3  Open-Loop  Command  Generator 
Tracker  (CGT)  (Refs  5,  32) 

The  command  generator  to  be  used  here  is  assumed 
to  be  given  by  a  linear,  time-invariant  state  differential 
equation  and  an  associated  output  equation  which  is 
referred  to  as  the  command  model.  After  defining  the  com¬ 
mand  model  and  obtaining  its  discrete-time  equivalent,  the 
concept  of  the  "ideal  state  trajectory"  is  introduced  and 
employed  to  achieve  necessary  equations  for  the  CGT  feed¬ 
forward  gains.  These  gains  assume  open-loop  implementation. 
Although  such  is  generally  an  unsatisfactory  implementation, 
the  feedforward  gains  thus  computed  are  identical  to  those 
needed  for  a  system  using  inner-loop  feedback,  as  is 
shown  subsequently. 

3.3.1  Command  Model .  The  command  model  is  defined 
by  a  linear,  time -invariant  state  differential  equation  and 
an  output  equation  as, 

(t)  =  Ax  (t)  +  B  u  (t)  (3-24) 

— m  — m- m  — m— m 

and, 

=  ^(t)  =  ^(t)  (3-25) 

The  dimensionalities  of  the  command  model  are, 
n„  =  number  of  model  states 

n 

r„  =  number  of  model  inputs 
PM  =  number  of  model  outputs 
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Since  it  is  desired  to  cause  the  system  outputs  to 

follow  those  of  the  command  model,  it  is  necessary  that 

the  number  of  outputs  for  each  be  equal:  p  «  p. 

M 

The  discretized  command  model  becomes 


*m(*W 


m  m  1 


(V 


(3-26) 


Y_(t.)  “  Cx  (t.)  +  D  u  (t.) 
^-m  1  — m-m  l  — m— m  i 


(3-27) 


where,  for  um  constant  over  a  sample  period 

A  T 
*  -  “in 

— m  e  (3-28a) 

T 

B  —  f  t(i)Bd:  (3-28b) 

— md  *'0  — m 

and  C  and  D  are  as  before  in  equation  (3-25) . 

— m  — m 

3.3.2  Ideal  Trajectory.  The  design  objective  for 
the  CGT  controller  is  to  force  the  system  and  command  model 
outputs  to  be  equal, 


(3-29) 

The  error  in  so  doing  at  time  t^  is 

e(ti)  =  y(tA)  -  ^(t^  (3-30a) 

or,  from  equations  (3-16)  and  (3-27) , 
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To  aid  in  deriving  the  feedforward  gains  for  the 
CGT  controller,  it  is  useful  to  formulate  an  idealization 
of  the  solution  in  which  the  "ideal  state  trajectory"  and 
the  command  model  state  trajectory  are  identical  for  all 
time.  To  do  so,  the  deterministic  ideal  system  state  and 
output  must  be  defined  so  as  to  satisfy  the  original  system 
state  equation  given  by  equation  (3-14) 

+  SoSl'V  +  5xdSd(ti'  ,3-31> 

and  the  corresponding  output  equation,  equation  (3-16) : 

*I(ti)  =  CXjtt.)  +  DyUjlt.)  +  ^(t.)  (3-32) 

where  x ^  and  are  the  ideal  state  and  output  vectors, 
respectively.  By  definition,  the  ideal  state  trajectory 
must  also  be  such  as  to  maintain  zero  error  between  the 
system  and  command  model  outputs: 

e(t.)  =  Zl(t.)  -  ym(t.)  =  0  (3-33) 

so  that 

<3-34» 

or,  substituting  expressions  for  ^  and  Zm  from  equations 
(3-32)  and  (3-27) ,  respectively,  gives 


Since  a  feedforward  control  law  operating  on  the 


command  model  states  and  inputs,  and  the  modeled  disturbance 
states,  is  required  for  the  CGT  controller,  an  additional 
constraint  is  imposed  for  tractability :  that  the  ideal 
state  and  input  vectors  be  a  linear  function  of  those 
three  vectors.  Representing  this  linear  function  by  par¬ 
titioned  submatrices  gives 


-11  —12  -13 

^21  -22  — 23_ 

x  (t.) 
— m  x 

W 

w 


(3-36) 


3.3.3  Open-Loop  CGT.  The  appropriate  driving 
input  to  the  system  is  the  u^.  vector,  so  solution  for  the 
values  of  the  constant  matrices  through  A^  gives  the 

necessary  feedforward  gain  matrices  for  the  open-loop  CGT 
controller.  Equations  allowing  solution  for  these  matrices 
are  developed  as  follows:  using  equation  (3-31)  and  augment¬ 
ing  the  forward  difference  expression  for  the  x^Jt^)  with 
the  output  equation  for  ^(t^)  of  equation  (3-32)  results 


*I<V 


(£-1) 

2d 

+ 

Ex 

xd 

c 

IryJ 

n,(t.) 

— Q  1 


(3-37) 


Now  substituting  the  assumed  form  for  xI(ti)  and  Ujtt^ 
given  by  equation  (3-36 )  yields 


*X(ti> 


(*-I) 

-XI 


-11  -12  -13 
-21  -22  -23 


x(t.) 
— m  l 

w 

Sd(ti> 


E 

—x . 
c 

E 

-y 


w 

(3-38) 


The  forward  difference  of  xT(t.)  may  also  be  obtained 

A  A 

directly  from  equation  (3-36)  as 


[£l(ti+l)_*I(ti)] 


[-ll  -12  -131 


x  (t .  . ) —x  ( t • ) 
— m  i+l  — m  l 

2d(ti+l,~2d,til 


(3-39) 


Assuming  to  be  constant  (or  slowly  varying  in  comparison 
to  the  sample  period) ,  then 


u  ( t .  i )  -  u(t.)  1  0  (3-40) 

— m  i+l  — m  i  — 

Using  the  state  models  for  x^  and  n^  of  equations  (3-26)  and 
(3-15)  respectively  with  driving  noises  deleted  leads  to 
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=  [-ll  -12  -13 ] 

(4-1)  B  0 

-m  -  -md 

0  0  0 

x  (t.) 
— m  i 

W 

0  0  (4n-I) 

(3-41) 


An  additional  equation  for  ^(t^)  is  obtained  from  equa¬ 
tion  (3-27),  wherein  yT(t.)  =  v  (t.),  so 

■*-1  i  ‘‘-m  i 


yT(t.)  =  (C  D  ] 
*■1  i  — m  — m 


x  (t.) 
— m  i 

u  (t.) 
— m  i 


(3-42) 


Performing  the  interior  matrix  product  of  equation  (3-41) 
and  augmenting  that  result  with  equation  (3-42)  gives 


-I^i+l*  "-l(ti) 


*I(ti> 


A.  .  (4  -I)  A, ,  B  A.  (4  -I) 
—11  —m  —  —11  — md  —13  — n  — 


— m 


D 

— m 


x  (t  ) 
— m  i 

W 


(3-43) 


Equations  (3-38)  and  (3-43)  now  give  two  different  expres¬ 
sions  for  the  forward  difference  of  Xjft^)  and  the  output 
^(tf),  both  involving  the  desired  feedforward  matrices. 
Setting  these  two  expressions  equal  yields 
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Partitioning  the  indicated  matrix  inverse  into  arrays 
of  the  same  dimensions  as  the  component  arrays  gives 


n  = 

-nj 

-21  j 

^12" 

’^22 

= 

(£-1) 

c 

^d 

D 

_ 

-xJ 

(3-47) 


where  for  this  development  the  right-hand  side  matrix  is 
square  since  y  and  u  have  the  same  dimension,  and  the 
inverse  is  assumed  to  exist.  It  is  possible  to  generalize 
this  result  using  matrix  pseudo- inverses,  but  that  is  not 
pursued  herein  (Ref  32) .  Thus  equation  (3-46)  can  be 
rewritten  as 


-11 

-12 

-12 

-21 

—22 

—23 

-11  -12 

-21  -22 


A.  .  <*  -I)  A.  .  B  A  _  (i  -I )  -E 
11  ™id  11  mH  13  n  x^ 


C 

— m 


-E 


(3-48) 


and  the  explicit  partitioned  sets  of  equations  are 


*11 

+  7T..  -C 
— 1 2— m 

(3-49a) 

*12 

=  211*11%  + 

— 1 2— m 

(3-49b) 

*13 

“  -11-13 

•  *11%  -  *12^ 

(3-49c) 

*21 

“  -21-11 ’In'-1 

+  l22£m 

( 3-49d) 
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(3-49e) 


-22  -21-1  l-md  +  —22— in 


—23  =  *21*13 "  ^2l4d  -  — 22^y  (3"49f) 


Equations  (3-49a)  and  (3-49c)  define  solutions  for  and 
A13  respectively,  while  the  other  equations  give  solutions 
for  the  remaining  A^^  matrices  which  follow  from  computa¬ 
tions  involving  matrices  of  known  value.  These  two  equa¬ 
tions  are  of  the  form 


X  =  AXB  +  C  (3-50) 

for  which  an  algorithm  for  solution  has  been  reported 
(Ref  4)  and  implemented  previously  in  applications  of  the 
CGT  design  technique  (Ref  10) .  Several  conditions  must 
be  satisfied  in  order  for  a  solution  to  exist.  These  con¬ 
ditions  are  (Refs  4,  5) 

K  *  *  X  *  1.  (3-51a) 

TT,i  m,j 

and 

Vi  •  \,.j  * (3-51b' 

For  equations  (3-49a)  and  (3-49c)  respectively,  where  X  . 

are  the  eigenvalues  of  the  7T^  partition  and  X^  ^  and 

X  .  are  the  eigenvalues  of  the  matrices  (4>  -I)  and 
n, j  — m  — 

(£n-I)  respectively.  Since  the  eigenvalues  of  tt^  are 
related  to  the  inverses  of  the  system  transmission  zeros, 
this  constraint  can  be  formulated  as  between  the  eigen¬ 
values  of  the  command  and  disturbance  models  and  the  system 
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transmission  zeros  (Refs  5,  14) .  In  addition,  no  discrete¬ 
time  transmission  zero  of  the  system  may  equal  one. 

With  A^  and  A^  determined,  and  thence  all  remain¬ 
ing  A^j  partitions,  the  open-loop  Command  Generator  Tracker 
control  law  is  obtained  from  the  lower  partitioned  equation 
of  equation  (3-36) : 


-I^i*  "  -21-m(ti)  +  -22-m(ti)  +  -23£d(ti* 


(3-52) 


The  open-loop  CGT  is  implemented  as  in  Figure  3-2. 


3.4  Optimal  Regula tor/PI 

Controller  (Ref  32l 

The  design  goal  in  employing  a  Proportional-plus- 
integral  (PI)  controller  is  to  generate  a  feedback  con¬ 
troller  which  will  maintain  the  system  output  defined  in 
equation  (3-5)  at  a  nonzero  commanded  value  with  zero 
steady  state  error  despite  unmodeled  constant  disturbances 
which  may  also  drive  the  system.  The  idea  is  well-known 
in  conventional  control  theory  wherein,  for  a  unity  feed¬ 
back  configuration,  the  designer  seeks  to  achieve  a  forward 
path  transfer  function  which  includes  a  pole  at  the  origin. 
This  is  often  achieved  by  employing  a  controller  with 
integration  of  the  error  in  the  control  variable  as  a 
feedback.  The  resulting  design  is  referred  to  as  a 
Proportional-plus-integral  controller  and  the  feedback 
system  is  described  as  having  the  "Type-1  property"  (Ref 
19)  due  to  the  e1  factor  in  the  forward  path  characteristic 
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polynomial  (where  a  is  the  Laplace  operator) . 

In  discrete-time  state  feedback  control  systems 
integration  is  generally  approximated  based  upon  the 
simple  Euler  integration  technique.  The  PI  action  can  be 
achieved  by  performing  such  pseudo-integration  on  the 
regulation  error  (the  differencing  of  the  system  and  desired 
outputs)  or  on  the  control  input  rate.  The  PI  controller 
developed  here  employs  pseudo-integration  of  the  control 
input  rate. 

Also,  the  controller  may  be  formulated  for  imple¬ 
mentation  in  either  of  two  forms:  the  "position  form"  or 
the  "incremental  form."  The  position  form  represents  the 
current  input  in  its  entirety  and  does  so  in  terms  of  the 
total  values  of  the  feedback  variables.  In  the  incre¬ 
mental  form  only  the  change  in  control  input  is  computed, 
and  it  is  in  terms  of  changes  in  the  values  of  the  feedback 
variables  since  the  preceding  sample  period.  The  incre¬ 
mental  form  for  the  controller  has  certain  advantages  over 
the  position  form  (Ref  7)  and  is  the  method  for  implementa¬ 
tion  which  is  developed  here. 

The  optimal  PI  controller  is  first  developed  from  a 
discrete-time  problem  formulation.  Subsequently,  the  tech¬ 
nique  for  translating  a  continuous-time  quadratic  cost 
formulation  to  the  appropriate  discrete- time  cost  function 
is  demonstrated. 
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3.4.1  Control  Difference  PI  Controller  (Ref  32) . 

The  control  difference  PI  controller  achieves  its  integral 
quality  based  upon  a  control  rate  pseudo- integration. 

Define  the  perturbation  control  variable  as 

6u(ti>  =  u(t±)  -  uQ  (3-53) 

where  u^  is  the  nominal  control  input  to  maintain  the  sys¬ 
tem  at  its  equilibrium  operating  point.  Then  the  forward 
difference  for  this  perturbation  control  variable  is. 

“  ^(t^  =  (“(ti+1) -uQ)-(  u(ti)-uQ) 

(3-54) 

or, 

6u(ti+1)  »  6u(ti)  +  (u(ti+1)  -utt.^ )  (3-55) 

This  can  be  thought  of  as  an  update  relation  for  6u(ti+1) 
based  on  Euler  integration  of  control  rate,  since  for 
Euler  integration  the  change  in  u  is 

Au(ti)  =  u(ti+1)  -  u(ti)  ~  T  •  jl(t^)  (3-56) 

where  T  is  the  controller  sample  period  and  il(t^)  is  the 
time  rate-of -change  of  the  input  u  at  the  sample  time  t^. 

This  Au  then  is  defined  to  be  the  control  differ¬ 
ence  or  "pseudo-rate."  Thus  equation  (3-55)  becomes 

<5u(ti+1)  =  Suft^)  +  Auft^  (3-57) 
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Defining  the  perturbation  state  flx(t^)  as 

6x (ti)  *  x(ti)  -  xQ  (3-58) 

where  Xq  is  the  nominal  state  vector  at  the  desired  oper¬ 
ating  point,  and  noting  that  the  perturbation  state  satis¬ 
fies  the  same  transition  equation  as  in  equation  (3-14) 
with  the  disturbance  state  and  the  noise  vector  deleted 
and  the  perturbation  control  variable  as  the  input  vector, 
an  augmented  state  description  can  be  formed.  The  control 
pseudo-rate  is  considered  the  input  to  an  augmented  pertur¬ 
bation  regulator  control  problem  defined  for  the  augmented 
perturbation  state  equation 


i 

*d 

6x(ti) 

0 

1 

6u(ti) 

Au(t±) 


(3-59) 


For  the  optimal  regulator  the  cost  criterion  to  be 
minimized  is 


N 

6x(ti) 

T 

in 

—12 

*1 

6x(t±) 

j  = 

6u(ti) 

ill 

-22 

Icn 

K> 

6u(ti) 

i=-l 

Au(ti) 

*1 

—2 

u 

Au(t±) 

Si(tN+l> 

T 

xf 

*11 

0 

ia'Vi’. 

0 

a. 

.‘s'ViL 
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where  X^  weights  state  deviations  from  the  nominal  Xq, 

X22  weights  control  deviations  from  the  nominal  Uq,  and 

U  weights  control  pseudo-rates.  The  weight  X,  applies 

rll 

to  the  state  deviation  at  terminal  time,  and  will  not  be 
used  further  since  the  regulator  to  be  used  here  is  based 
on  the  infinite-time  steady-state  regulator  problem.  As  is 
discussed  in  the  next  section,  the  weights  and  S2 
will  arise  as  the  continuous-time  cost  formulation  is  con¬ 
verted  to  a  discrete-time  cost.  Finally,  the  cross  weights 
X12  between  the  state  and  input  vectors  will  be  non-zero 
if  the  system  output  equation  includes  a  non-zero  or 
may  arise  due  to  the  discretization  of  the  continuous-time 
cost.  Note  the  index  for  the  summation  begins  at  -Is  this 
serves  to  weight  the  potentially  large  control  difference 
which  may  occur  at  the  initial  time  due  to  a  change  in 
setpoint  (Ref  32)  (note  Au(t_1)  =  u(tQ)  -  u(t_1)  by  equa¬ 
tion  (3-56)  . 

3.4.2  Continuous-Time  Cost  (Refs  15,  32).  Since 
the  system  to  be  controlled  is  a  continuous-time  system, 
and  since  its  behavior  is  important  at  all  time  and  not 
merely  at  the  controller  sample  times,  the  cost  function 
appropriate  to  the  regulator  design  is 


J  *  h 


f 


'N+l 


X(t) 

T 

u(t) 

W 

—XX 

T 

*xu 


w 

—xu 

W 

-uu 


X  ( t ) 

u(t) 


dt 


(3-61) 
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where. 


x 


(3-62a) 


u 


Au 


(3-6 2b) 


(3-62c) 


(3-62d) 


and,  for  this  development  the  cross  weight  Wxu  is  assumed 
to  be  zero,  since  its  inclusion  is  rarely  necessary  to 
achieve  control  objectives.  Note  also  the  problem  is 
posed  as  deterministic  since  its  true  stochastic  nature 
does  not  impinge  on  the  optimal  regulator  design  due  to 
certainty  equivalence.  Terms  in  the  residual  cost  associ¬ 
ated  with  the  stochastic  driving  noise  in  the  stochastic 
optimal  controller  cost  formulation  are  independent  of  the 
choice  of  control  function  applied  and  thus  do  not  affect 
the  choice  of  optimal  control  (Refs  2,  32) . 

Furthermore,  since  the  CGT  design  objective  is  to 
drive  the  system  so  that  its  output  tracks  the  model  output, 
it  is  appropriate  that  the  quadratic  weighting  matrices 
specified  should  apply  to  the  system  outputs  and  inputs. 

Thus  defining  the  weight  on  output  deviations  as  Y  and  on 


50 


input  magnitudes  as  U^,  the  components  of  the 
are  obtained  as  (Ref  15) 

X_  *  CTYC 


^*11 

X  =  U  +  DTYD 

~°22  -y  -y—y 

x  =  ctyd 
~ci2  ~  —  y 


matrix 


(3-63a) 

(3-63b) 

(3-63c) 


where  C  and  Dy  are  as  defined  in  equation  (3-5)  and  Y  and 
Uy  are  positive  semidefinite  and  positive  definite,  respec¬ 
tively.  Therefore,  W  is  positive  semidefinite  while  the 

aA 

U  weighting  matrix  is  required  to  be  positive  definite. 

w 

After  forming  W^_  according  to  equation  (3-63)  the  designer 
may  then  modify  any  elements  to  achieve,  for  example,  weights 
on  some  states  not  included  directly  in  the  output  defini¬ 
tion  of  the  CGT  design  objective. 

In  order  to  employ  this  continuous-time  cost  func¬ 
tion  of  equation  (3-61)  for  solution  of  the  discrete-time 
optimal  regulator  it  is  necessary  to  obtain  the  correspond¬ 
ing  discrete-time  cost  function.  Begin  by  conceptually 
dividing  the  control  interval  t q  to  tfJ+^  into  (N+l)  control 
intervals  of  duration  equal  to  the  intended  controller 
sample  period  T .  The  cost  can  then  be  expressed  as 


s  /  -t.+T 

■e  i  r  ■ 

i=0  \i 


h [xT ( t ) Wxxx ( t ) +uT ( t ) Wuuu ( t ) ] dt 


) 

/(3-64) 
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where  u(t)  is  assumed  constant  over  a  sampling  period  and 
x(t)  satisfies 


x(t)  *  l^xlt..)  + 


(3-65) 


for  all  te[ti,t^+T)  and  where,  from  equation  (3-59), 


$  ^ 

$  = 

”  [£  !  I 


(3-66a) 


—6  =  ” 


(3-66b) 


The  discrete-time  cost  function  then  becomes 


N 

J  =  lxT(ti)X6x(ti)+uT(ti)U6u(ti) 


+  2xT(ti)S(5u(ti)] 


(3-67) 


where , 


f  ti+T 

s«  -  Jt  1  iI(T>Mxxi6(T)dT 


(3-68a) 


/ti+T 

[bTw  B.+W  ]dx 
— 6— xx— o  — uu 


/t^T 


(3-68b) 


(3-68c) 


with 


*(x)  1  B.(t) 

»  (,)  .  -- — — — 

2  !  i 


(3-68d) 


Bj(t) 


=  f  $  (o)Bc 

J  0 


(3-68e) 


Equation  (3-67)  may  be  expressed  more  compactly: 


!  i 


j  _  y>  S<ti>  T  x6  i5  xttA) 
J  "  u6  ult,) 


(3-69) 


Note  that  the  cross  weighting  matrix  S  has  been 
introduced  into  the  cost  function  by  the  discretization 
process.  In  order  to  obtain  an  equivalent  discrete-time 
cost  function  with  no  cross  weighting  (to  allow  use  of 
standard  Riccati  equation  solvers  that  assume  such  a  form) , 
define  a  new  system  (Ref  29) 


2£<ti+l)  =  *6*'  j.>  +  ®6^'(ti.) 


(3-70) 


for  which 


-1  T 

Is  - 


(3-71a) 


u'Ct^  =  ult^  +  U(j1S^x(ti) 


(3-71b) 


and  for  which  the  corresponding  cost  function  to  be  mini¬ 


mized  is 


.  -  ■*  »  r>- 


with 


N 

*5  +  u"T(ti)0'u"(ti)l  (3-72) 

i=0 


-1  T 

S«2«  5* 


(3-73) 


If  the  system  of  equation  (3-70)  is  either  con¬ 
trollable  or  stabilizable,  letting  N-*»  leads  to  a  steady 
state  solution  of  the  discrete  Riccati  equation  represented 
as  K^,  where  satisfies  (Ref  15) 


H 


*  SjSRtj 


(3-74) 


and  the  optimal  feedback  control  is 

u  *  "*  ( t . )  =  -G*'x(t.)  (3-75a) 

1  v  JL 

where 

-  [“6+^K'Bj]-1  B^K't'  (3-75b) 

The  corresponding  optimal  feedback  gain  matrix  for  the 
original  state  system  (x,  u)  is  (Ref  29) 

G*  =  G*'  +  U"1  sT  (3-76) 

— c  — c  —o  —o 

Remembering  the  definition  of  x  of  equation  (3-62a)  G£ 
can  be  written  in  partitioned  form  as 

2*  -  [Gi  !  G*  ]  (3-77) 

-Ci  -c2 
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and  the  optimal  control  input  is 


Au* (t . )  =  -G*  6x(t.)  -  G*  6u (t . ) 
—  i  —Ci  l  — c  2  —  l 


(3-78) 


Combining  this  expression  with  the  definition  of  the  con¬ 
trol  difference  given  in  equation  (3-56)  gives 


6^*(ti+l)  = 


[^l' 


•G*  ] 
2 


6x (ti) 
6u(t±) 


(3-79) 


3.4.3  Achieving  Integral  Control .  These  results 
do  not  yet  provide  the  desired  integral  characteristic  for 
the  controller.  Such  can  be  achieved  by  manipulation  to 
a  form  emulating  that  of  a  continuous-time  PI  controller: 

u*(t)  =  -Kxx(t)  +  Kz  f  t[yd-Z(x)]dT  (3-80) 

where  is  the  desired  output  and  y(t)  is  the  actual 
system  output  defined  in  equation  (3-5)  but  with  the  dis¬ 
turbance  term  deleted.  On  the  analogy  of  Euler  integra¬ 
tion  of  the  tracking  error  over  each  sample  period,  a 
discrete-time  equivalent  would  be 

i-'l 

u*(t.)  =  -Kxx(t.)  +  Kz  y^fra-gU.)]  (3-81) 

This  can  be  expressed  equivalently  in  incremental  form, 
wherein  the  input  at  time  t^+^  is  obtained  as  an  update 
on  the  input  at  time  t^: 
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uMt.+1>  =  u*(t.)  -  Kx[x(ti+1)-x(t.)] 


+  (3-82) 

Returning  to  the  perturbation  state  notation  and  noting 
that  the  perturbation  in  output  ^  is 


6Z<ti)  =  ^  (3-83) 

and 

6x(ti)-  =  C6x(ti)  +  Dy6u(ti)  (3-84) 

then 

fiu*<ti+1)  «  6u*(ti)  -  Kx[6x(ti+1)-6x(ti)  ] 

-  K  (C6x(t.)  +Dv6U*(t.)]  (3-85) 

z  i  y  x 


Employing  the  expression  for  the  forward  difference  of 


6x(ti+1)  -  Sxtt.^)  =  [£-I_]  6x  (t^) 
+  B<a6u*(ti) 

allows  equation  (3-85)  to  be  rewritten  as 


(3-86) 


fiu*  ( t 


i+1 


<5u*  (ti) 


(I-I) 

fix (ti) 

1 - 

in 

fiu*(t±) 

(3-87) 


56 


Since  the  expressions  for  6u*(t^+^)  in  equations  (3-79) 
and  (3-87)  are  to  be  equal  and  employ  constant  gains,  it 
follows  that 


(K  K  ] 
—x  — z 


(1-1) 

C  D 

L  ~  -ZJ 


[G*  G*  ] 
"®1  ~C2 


(3-88) 


and  thus  that 


[K  K  ]  = 
— x  — z 


(*-I) 
_  £ 


-1 


I6J  ^  1 

-c2 


(3-89) 


Remembering  the  matrix  composed  of  partitions  defined 
by  equation  (3-47)  and  writing  the  partitioned  equations 
explicitly  gives  values  for  the  feedback  matrices  of 


■  sj/n  + 


(3-90a) 


and, 


K„  =  G*  £,  9  +  G*  jr 
z  1  c  ^  x  ^  c  2  «  * 


(3-90b) 


The  final  equation  for  the  PI  regulator  implemented 
in  incremental  form  and  for  the  sample  time  t^  is, 

«u*<t±)  =  5u*(ti_1)  -  Kx[<Sx(ti)-5x(ti_1)] 


+  .KzCZd(ti)  -  1(^)1 


(3-91) 
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where  now  the  desired  output  is  allowed  to  be  changing 
and  the  time  indices  for  and  y  differ  since  the  pseudo¬ 
integral  of  the  error  must  include  new  error  introduced  at 
time  t^  by  the  changed  value  of  desired  output  (Ref  32) . 
Note  that  this  is  directly  related  to  the  lower  limit  on 
the  sum  being  -1  for  this  formulation. 

3.5  Closed-Loop  CGT/PI 

The  developments  of  Sections  3.3  and  3.4  may  now  be 
combined  in  a  closed-loop  CGT/PI  controller.  It  will  drive 
the  system  so  as  to  achieve  matching  of  the  actual  system's 
outputs  with  the  command  model  outputs  in  steady  state 
despite  possible  errors  in  the  models  used  for  design  and 
despite  unmodeled  constant  disturbances  which  may  drive 
the  system  in  addition  to  those  for  which  rejection  was 
designed. 


3.5.1  Perturbation  Regulator .  Returning  to  the 
concept  of  the  ideal  trajectory,  define  the  control  differ¬ 
ence  for  u^t^)  as 


-  HI(ti+1>  -  vv 


(3-92 


which,  using  equation  (3-52),  can  be  rewritten  as 


Note  that  henceforth  the  assumption  of  equation  (3-40) , 


wherein  the  command  model  input  u  is  assumed  constant,  is 

— m 

not  needed. 


Define  the 

set  of 

perturbation  variables  as 

6x(ti) 

=  x(ti) 

"  -I<ti) 

(3-94a) 

<Su  (t^ 

=  u(ti) 

-  ^(t.) 

(3-94b) 

6y(ti) 

=  *<t±) 

-  zI(ti) 

(3-94c) 

and,  remembering  the  definition  of  Au(t^)  given  by  equa¬ 
tion  (3-56) , 


6Au(t±)  =  Au(ti)  -  Au  (t±)  (3-94d) 

The  augmented  system  perturbation  state  equation  is 
the  same  as  in  equation  (3-59)  but  with  6Au  replacing  Au 
in  that  equation.  Similarly  the  steady-state  optimal  con¬ 
trol  solution  is  in  equation  (3-78)  with  the  Au  substitu¬ 
tion: 


<5Au*(t.)  =  -G*  6x(t.)  -  G*  <5u(t.)  (3-95) 

i  -c1  i  c  2  ~  i 

Substituting  the  expression  of  equation  (3-56)  for  Au(ti) 
into  equation  ( 3— 94 d)  gives 

6Au(ti)  =  u(ti+1)  -  u(t±)  -  AuI(ti)  (3-96) 

Using  equation  (3-96)  to  replace  6Au(t^)  in  equation 
(3-95),  shifting  the  time  argument  backward  one  sample, 
and  making  substitutions  for  6x(t^_^)  and  6u(t^_^)  with 
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time-shifted  versions  of  equations  (3-94a)  and  (3-94b) 
respectively  gives  an  equivalent  expression  for  equation 
(3-95)  as 


uMt^  -  u*(ti_1)  -  Au^t^j.) 

-  -<S1l*(ti-l)-2l(ti-l,> 

-G^2[^(ti_1)-uI(t._1)]  (3-97a) 

or 

u*(ti)  =  u* (ti_1)  + AuI(ti_1) 

-  G*  (3-97b) 

2 

3.5.2  Achieving  Integral  Control .  Making  substitu¬ 
tions  for  AuI#  Xj.,  and  Uj  from  equations  (3-93)  and  (3-36) 
respectively  and  applying  appropriate  shifts  in  time  argu¬ 
ment  yields 


+  A-,[u  (t.)-u  (t.  .)] 
—22  — m  i  — m  l-i 


+  -23tlld(ti)“I1d(ti-l)  1 


-  ^[x(ti_1)-A11xin(ti_1) 


“  — 12— m^i-l^  "  ^13^d(ti-l,] 


60 


which  is  in  the  desired  incremental  form,  but  still  lacks 
the  integral  property  characterizing  the  PI  controller.  As 
in  Section  3.4  this  property  can  be  achieved  using  a  con¬ 
troller  of  form 


<Su*  (t^)  =  6u*(ti_1) 


[£x  £»' 


6x(ti)  -  6x(ti_1) 


6Z(ti-l> 


(3-99) 


Employing  the  expressions  of  equations  (3-94)  and  remember¬ 
ing  that  Yj  =  ym  by  definition  in  equation  (3-34),  equation 
(3-99)  can  also  be  written  as 


u*(t.)  -  ul(t.)  =  «‘(tw)  -  u^t.^) 


-  KjjKxft^-Xjfti) )  -  (x(ti_1) -Xjtt^^  j^)  )  ] 

-i.H'W  ,3-100a) 


or,  after  rearranging  terms. 


uMt^)  =  u*(ti_1)  +  IuI(ti)-uI(ti>_1)  ] 


-  Kxlx(ti) -x(ti_1) ]  +  KxlxI(ti)-xI(ti_1)] 


(3-100b) 


\ 


Using  equation  (3-36)  and  differencing  and  Uj 
at  time  arguments  and  t^_^  leads  to 


r~"  — * 

x  (t.)-x  <t.  .r 

Si'V'Si'Vi' 

-11 

—12 

—13 

— m  x  -in  i-1 

W-aa(W_ 

(3-1011 

*21 

-22 

-23_ 

Equation  (3-101)  allows  substitutions  for  two  terms 
of  equation  (3-100b)  to  be  nade  as  follows: 


+  -23  1 


(3-102a) 


and  the  right-hand  side  can  be  rearranged  to  yield 


'  I£jc-11+-211 


+  (^12+^221 


+  t-x-13+-23*  (nd<til-nd(ti_il) 


(3-102b) 


To  ensure  consistency  in  the  equations  for  the 
ideal  state  trajectory  in  the  face  of  changes  in  the  value 
of  the  model  input  and  to  improve  the  initial  transient 
performance,  it  is  necessary  to  apply  direct  feedthrough 
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of  to  both  x  and  v  so  that  equations  (3-26)  and 
— in  — in  --m 

(3-27)  are  actually  implemented  as  (Ref  5) 


Sm'Vl1 


»X  ( t . )  +  B  u  (t.  .) 
- m— m  i  — m^—m  l+l 


(3-26') 


v  (t.)  =  C  x  (t.)  +  Du  (t.  .) 
m  i  — m— m  i  — rti-m  l+l 


(3-27') 


where  the  time  argument  of  u  has  been  advanced  by  one 

— m 

sample  period.  Thus  the  model  states  and  outputs  cannot  be 
updated  between  sample  times,  but  must  await  the  new  command 
input  sample. 


3.5.3  CGT/PI  Control  Law.  Using  equations  (3-27') 
and  (3-16)  with  backward  time  shifts  for  both  and  deletion 
of  the  disturbance  term  for  the  latter  gives  replacements 
for  ym  and  y  in  equation  (3-100b) .  Combining  these  with 
the  replacements  of  equation  (3-102)  transforms  equation 
(3-100b)  to  the  final  incremental  form  CGT/PI  control  law: 
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In  implementing  the  CGT/PI  control  law  in  this 
form,  it  is  necessary  to  consider  the  initial  conditions  to 


provide  proper  start-up:  if  any  of  the  system,  model,  or 
disturbance  states  are  non- zero  at  start-up  then  the  initial 
control  to  be  applied  u*(t^)  must  be  based  upon  the  total 
values  of  the  state  variables  at  the  initial  sampling 
instant.  This  can  easily  be  achieved  by  employing  the 
incremental  form  of  equation  (3-103)  with  the  states  corres¬ 
ponding  to  index  t^  set  to  the  appropriate  initial  values 
and  the  states  at  index  t^_^  and  all  control  input  vectors 
must  be  set  to  their  previous  steady  state  values.  At  all 
future  sample  times,  the  model  states  are  updated  as  in 
equation  <3  26')  followed  by  an  update  of  the  control  input 
according  to  equation  (3-103) ,  assuming  that  system  and 
disturbance  state  updates  are  also  available  as  obtained 
either  directly  from  measurements  or  as  estimates  from  a 
Kalman  filter  update. 

3.6  Kalman  Filter  (Ref  31) 

The  CGT/PI  control  law  of  equation  (3-103)  entails 
knowledge  of  the  entire  system  and  disturbance  state  vectors 
employed  in  the  design  model  definition  of  equation  (3-1) . 
Typically,  of  course,  not  only  are  these  not  all  available 
from  direct  measurements  taken  of  the  system,  but  the  mea¬ 
surements  which  are  available  are  corrupted  by  noise.  In 
order  to  provide  estimates  of  the  states  a  Kalman  filter 
is  employed. 
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Although  from  the  outset  of  the  theoretical  develop¬ 
ment  of  this  chapter  the  stochastic  nature  of  the  design 
problem  has  been  recognized,  it  has  not  formally  impacted 
any  of  the  designs  thus  far  presented.  By  formulating  the 
problem  definition  according  to  the  so-called  LQG  assump¬ 
tions  ("Linear-Quadratic-Gaussian") — a  stochastic  optimal 
control  design  for  a  system  represented  by  a  Linear  system 
model,  employing  a  Quadratic  cost  criterion,  and  with 
Gaussian  driving  noises,  it  has  been  possible  to  invoke 
the  "Certainty  Equivalence"  for  the  controller/filter 
design  (Ref  2) .  Certainty  equivalence  states  that  the 
optimal  stochastic  controller  for  a  system  designed  accord¬ 
ing  to  the  LQG  assumptions  consists  of  an  optimal  feedback 
controller  designed  independently  of  the  stochastic  nature 
of  the  system  and  which  is  equivalent  to  the  corresponding 
deterministic  optimal  LQ  controller  as  already  discussed, 
and  a  Kalman  filter  for  the  system,  independent  of  the  con¬ 
troller  design,  to  provide  the  needed  state  estimates.  The 
controller  and  filter  may  thus  be  designed  independently 
using  techniques  and  computations  appropriate  to  each,  then 
combined  and  each  design  further  tuned  to  achieve  best 
overall  response  for  the  ensemble  configuration. 

The  development  here  begins  with  the  definition  of 
the  system  model  for  the  Kalman  filter  design.  It  then 
proceeds  to  the  time  propagation  and  measurement  update 
equations  for  the  filter  covariance  and  gain,  to  determina¬ 
tion  of  the  steady  state  Kalman  filter  gain  matrix,  and 
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finally  to  the  corresponding  state  propagation  and  update 
equations . 


3.6.1  Design  Model .  The  Kalman  filter  design  for 
use  with  the  CGT/PI  controller  of  equation  (3-103)  must 
provide  state  estimates  for  both  the  system  and  the  dis¬ 
turbance  of  equation  (3-1) .  Therefore  the  appropriate 
system  model  for  design  is  based  upon  the  augmented  design 
model  of  equation  (3-9) .  The  state  transition  equation  is 
that  of  equation  (3-14)  while  the  measurement  equation 
becomes 


z(t.)  =H  x  (t.)  +  v  ( t . )  (3-104) 

—  l  —a— a  i  —  l 

where  z  and  v  are  as  defined  in  equation  (3-6),  x  is  as 

^  d 

defined  in  equation  (3 -9a) ,  and  H_  is  given  by 

'  'a 

H  =  [H  |  Hi  (3-105) 

—a  —  1  — n 

where  H  and  H  are  as  defined  in  equation  (3-6)  . 

—  — n 

3.6.2  Steady-State  Kalman  Filter  (Ref  31).  As  for 
the  controller  discussed  previously,  the  Kalman  filter 
desired  for  implementation  is  to  employ  constant  gains. 

If  the  augmented  design  model  defined  above  is  stabilizable 
and  detectable  then  a  steady-state  solution  to  the  filter 
Riccati  equation  can  be  achieved  which  determines  a  corres¬ 
ponding  set  of  filter  gains  to  be  used  in  a  constant  gain 
filter.  Further,  this  constant-gain  steady-state  Kalman 
filter  is  asymptotically  stable  (Ref  32) .  A  constant-gain 
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steady-state  filter  to  be  used  in  conjunction  with  a 
constant-gain  steady-state  controller  is  motivated  here 
because  the  system  is  quasi-static;  the  filter  and  con¬ 
troller  have  short  transients  at  the  beginning  and  end, 
respectively,  of  the  control  period  of  interest  compared 
to  the  long  time  of  essentially  steady-state  performance; 
and  for  aircraft  applications  minimal  complexity  of  the 
controller/ filter  is  desirable. 

In  all  equations  to  follow  superscripts  of 
and  "+"  on  the  time  argument  for  a  matrix  indicate  the 
value  for  that  matrix  at  the  given  time  argument  is  either 
before  or  after  a  measurement  update,  respectively. 

The  filter  covariance  is  propagated  forward  in 
time  using  the  equation  (Ref  31) 


where  R  is  as  defined  in  equation  (3-7) . 

The  measurement  update  equation  for  the  filter 
covariance  is 
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(3-108) 


—a ^ -  WSA'V 


In  steady  state  the  values  of  P  (t.)  and  of 

d  X 

Pa(t^_^)  will  be  equal  and  so,  combining  equations  (3-106) 

and  backward  time  shifted  versions  of  equations  (3-108) 

and  (3-107)  gives  the  steady  state  expression  for  the 

covariance  P  as 
—a 


4 


P  HT  [H  P  H*+Rry  P  }4>T  +  Q 
— a — a  — a — a — a  —  — a — a  — a  —a 


d 

(3-109) 


and  the  steady -state  Kalman  filter  gain  K_  is 

— F 

IF  -  Ia5a  ISjafi^Rl  *1 


(3-110) 


3.6.3.  State  Estimates.  Employing  the  steady- 
state  Kalman  filter  gain  of  equation  (3-110) ,  the  augmented 
state  vector  propagation  and  measurement  update  equations 
are 


*a(ti> 


*  x  (tT  .) 
—a— a  l-l 


♦  5adS*(ti-l> 


(3-mi 


and 


A  + 
*a(tt> 


x  (t.) 
—a  l 


(3-112) 


where  B  is  as  defined  in  equation  (3-12)  and  u*  is  the 
ad 

CGT/PI  control  input  of  equation  (3-103)  .  Initialization 
is  achieved  by  setting 


x,(°) 

—a 


x(0) 

Sa<°» 


(3-113) 
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where  x(0)  and  n^tO)  are  the  expected  values  of  the  system 
and  disturbance  state  vectors  at  the  initial  time. 

3.7  CGT/PI/KF  Control  Law 

The  developments  of  this  chapter  provide  the  neces¬ 
sary  design  equations  for  the  CGT/PI/KF  controller  and 
also  the  equations  for  implementation. 

The  elements  of  the  design  process  consist  of 

(1)  designing  a  CGT  controller  providing  feedforward  gains, 

(2)  designing  an  optimal  PI  regulator  for  inner-loop  con¬ 
trol,  and  (3)  designing  a  Kalman  filter  for  state  estima¬ 
tion  needed  for  application  of  the  controller's  gains. 

While  the  mutual  separability  of  these  three  designs  is  a 
crucial  aspect  of  the  methodology  and  serves  to  make  a 
successful  CGT/PI/KF  design  feasible,  it  must  be  realized 
that  the  quality  of  the  final  design  depends  on  the 
behavior  of  the  three  design  elements  acting  in  concert. 
Thus  some  design  tuning  is  both  appropriate  and  to  be 
expected  based  upon  the  final  integrated  system-controller- 
filter  response. 

The  final  form  of  the  CGT/PI/KF  controller  is  best 
represented  by  specifying  the  equations  needed  during  a 
typical  control  input  generation  cycle.  Assume  that  the 
controller  is  transitioning  from  its  having  just  generated 
a  control  input  u*  at  time  t^  through  the  computations 
necessary  to  yield  the  next  control  input  at  time  t^+^. 
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Prior  to  sample  time  t^+^: 

Propagate  the  augmented  state  vector  forward  in 
time  using  equation  (3-111) : 


x  (t:.,,)  ■  $  x  (t.)  +  B  u*(t.) 
—a  i+i  —a— a  1  ~ad  1 


(3-114a) 


At  sample  time  t^+^: 

Incorporate  the  new  measurement  using  equation 

(3-112) : 


£alti+l>  *  £a(ti+l>  +  *F  U(tt+i>-Sa35a(ti+1>] 


(3-114b) 


Propagate  the  command  model  state  vector  forward 
in  time  using  equation  (3-26^): 


x  (t..,)  =  $ _x  (t.)  +  B  u  (t.  ,) 
— m  i+l  — m-m  i  — m^—m  i+i 


(3-115) 


Compute  the  new  control  input  using  the  partitioned 
state  estimates  of  equation  (3-114)  and  the  command  model 
states  of  equation  (3-115)  in  equation  (3-103) : 


~  /  ,  +* 


S*«W  = 


+  <£x^l2^2211%(ti)-Hm(ti-l)) 


+  x— 13+— 23^  ^(ti+l*  _^d(ti^  ^ 


*z 


tC  D  ] 
— m  — m 


— m( ti+l) 


-[CD] 

— y 


A  + 

x(t  +) 
u* ( t ^) 


(3-llb) 
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IV.  CGT/PI/KF  Design  Evaluation 

4 . 1  Introduction 

It  was  shown  in  the  course  of  the  theoretical 
development  of  Chapter  III  that  the  CGT/PI/KF  design 
actually  consists  of  three  separate  designs:  the  regulator 
with  PI  control  action,  the  open-loop  CGT  or  closed-loop 
CGT/PI  controller,  and  the  Kalman  filter.  Each  of  these 
designs  is  best  evaluated  according  to  criteria  specifically 
related  to  the  task  that  each  is  to  perform. 

For  the  controller,  relevant  considerations  include 
the  closed-loop  system  poles,  the  values  of  the  feedforward 
and/or  feedback  gains,  and  the  time  response  of  the  con¬ 
trolled  system  states,  outputs,  and  control  inputs  in 
either  unforced  or  forced  input  conditions.  For  the 
filter,  relevant  considerations  include  the  poles  of  the 
filter,  the  values  of  the  filter  gains,  and  the  filter's 
estimation  error  behavior.  The  specific  elements  of  these 
evaluations  are  discussed  in  the  following  sections. 

An  overall  evaluation  of  the  system  performance 
for  the  CGT/PI/KF  closed-loop  controller  is  necessary  to 
tune  the  entire  design  properly  and  judge  its  ultimate  per¬ 
formance.  In  particular,  the  evaluation  elements  suggested 
here  will  demonstrate  controller  performance  alone  (CGT/PI) , 
as  if  the  needed  system  and  disturbance  feedback  and 
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feedforward  states  are  available  immediately  and  with 
perfect  accuracy.  The  true  CGT/PI/KF  controller  will 
suffer  some  degree  of  degraded  performance  due  to  the 
dynamics  of  the  Kalman  filter's  state  estimation  and  also 
a  slightly  increased  delay  in  control  generation  due  to 
the  needed  filter  computations  (Ref  2) .  These  effects  may 
be  evaluated  and  final  integrated  tuning  of  the  design 
achieved  by  employing  a  performance  evaluation  structure 
of  the  type  discussed  in  Reference  32.  Although  this 
thesis  effort  did  not  generate  the  computer  coding  to  per¬ 
form  such  a  full  controller/filter  performance  evaluation, 
future  efforts  will  extend  this  work  to  accomplish  this 
objective. 

4.2  PI  Regula tor  Evaluation 

The  discrete-time  poles  of  the  closed-loop  system 
incorporating  the  optimal  gain  G*  of  equation  (3-76)  and 

C 

assuming  perfect  state  knowledge  are  computed  from  the 
matrix 


4«CL=[is  (4-1» 

with  and  as  defined  in  equation  (3-66) .  The  equiva¬ 
lent  continuous- time  poles  are  then  computed  using  the 
inverse  of  the  relation  between  the  z  and  e  transforms 

eT 
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z  =  e 


(4-2) 


in  which  a  and  a  are  respectively  the  discrete  and 
continuous-time  complex  poles  and  T  is  the  controller 
sample  period.  This  mapping  is  for  the  primary  strip  in 
the  s -plane  only  and  does  not  consider  any  possible 
aliasing  effects  (Ref  28) . 

An  important  consideration  in  evaluating  the  con¬ 
troller  design  is  the  magnitudes  of  the  feedback  gains  Kx 
and  Kz  of  equation  (3-90)  .  Gains  which  are  relatively 
very  small  for  specific  states  may  indicate  that  reduction 
techniques  may  be  usefully  applied  to  the  design  model. 

On  the  other  hand,  gains  that  would  tend  to  cause  large 
control  inputs  may  indicate  an  unsatisfactory  design  due 
to  the  possibility  of  control  saturation  or  rate  limiting. 

Although  knowledge  of  the  closed-loop  poles  and 
examination  of  the  feedback  gain  matrices  does  provide 
insight  into  the  system  behavior  to  be  expected,  often  the 
most  illuminating  information  can  be  determined  from  the 
system's  time  response  to  initial  conditions  on  the  states. 
The  time  response  can  be  readily  simulated  using  the  dis¬ 
cretized  deterministic  state  transition  equations  of  either 
the  design  or  the  truth  model  and  driven  by  the  control 
input 


u*  ( t±) 


u.*  (t±_1)  -  Kx  [x(tA)  -x  (ti_1)  ] 


(4-3) 
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obtained  from  equation  (3-103)  by  deleting  terms  due  to 
the  CGT  feedforward  control.  Note  that  this  controller 
form  is  valid  only  for  the  PI  regulator  with  zero  refer¬ 
ence  input  and  in  response  to  non- zero  initial  conditions 
on  the  states.  Plots  of  the  time  histories  of  the  states, 
the  outputs,  and  the  generated  control  inputs  then  may  be 
made.  These  allow  evaluation  of  the  quality  of  the  regula¬ 
tion  achieved — speed  and  damping  of  the  state  and  output 
response,  and  the  magnitudes  and  rates  of  the  control 
inputs  actually  required. 

4.3  CGT  or  CGT/PI  Evaluation 

As  for  the  PI  regulator  evaluation  discussed 
above,  consideration  of  the  magnitudes  of  the  feedforward 
gains  may  provide  useful  insights  into  the  CGT  design 
result.  For  the  open-loop  CGT  the  relevant  gains  are 
A21  acting  on  the  command  model  states  and  A23  acting  on 
the  disturbance  states  in  the  open-loop  CGT  control  law 


u(ti) 


H'Vi' 


+  A-. [x  (t . )-x  (t .  .)] 
—21  — m  l  — m  l-l 


+  A,_[u  (t.)-u  (t.  .)] 
—22  — m  i  — m  i-i 


+  A23[nd<ti)-nd(t._1)] 


(4-4) 


obtained  from  equation  (3-103)  by  setting  the  PI  controller 

gains  K  and  K  to  zero.  For  the  CGT/PI  controller  of 
—x  —  z 

equation  (3-103)  the  relevant  gains  are 


(4-5a) 


ts 


K  -  SAl  +  —21 

in 


-xu  “  — x— 12  +  -22 


and 


KXr  «  KxA13  +  A23 


(4-5b) 


(4-5c) 


acting  on  the  coiranand  model  states  and  inputs#  and  distur¬ 
bance  states#  respectively. 

Evaluation  of  the  CGT  or  CGT/PI  control  system  time 
response  behavior  due  to  step  inputs  on  the  command  model 
inputs  is  crucial  to  judging  the  quality  of  the  controller 
design.  Here#  plots  of  the  responses  of  the  system's 
states,  outputs,  and  inputs  along  with  the  outputs  of  the 
command  model  allow  the  designer  to  evaluate  the  merits 
and/or  deficiencies  of  the  design. 

The  time  response  of  the  system  may  be  readily 
simulated  using  the  deterministic  state  transition  equa¬ 
tions  for  the  design  or  truth  models,  the  command  model, 
and  the  control  law  of  equation  (3-103)  with  K  and  K  set 
to  zero  for  the  open-loop  CGT  or  to  their  PI  design  values 
for  the  CGT/PI. 


4.4  Kalman  Filter  Evaluation 

The  state  equations  of  the  Kalman  filter  algorithm 
may  be  obtained  by  rewriting  equation  (3-112)  as  (Ref  31) 


x(tT)  «  U-KFHaJ  x(ti)  +  KFz(ti) 


(4-6) 
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l 


/V  — 

then  substituting  the  expression  x(t^)  of  equation  (3-111) 
to  yield 

5<tf  *  (I-SF»al  «aS(ti-l' 

+  II-KpHa]  Ba  uMt^)  +  KpZtt^  (4-7) 
d 

From  equation  (4-7)  it  can  be  seen  that  the  state 
transition  matrix  for  the  filter  estimate  propagation  is 

— kf  *  ‘I-SfS*1  ~a  (4-8» 

The  filter's  discrete-time  poles  are  computed  as  the  eigen¬ 
values  of  $KF  and  their  continuous-time  equivalents  com¬ 
puted  by  the  method  mentioned  in  Section  4.2.  The  magni¬ 
tudes  of  the  Kalman  filter  gains  may  be  evaluated  from  the 
Kp  matrix  of  equation  (3-110) .  While  the  filter  poles  and 
the  magnitudes  of  the  filter  gains  may  provide  some  useful 
insights  into  the  filter  performance,  often  the  greatest 
insight  is  achieved  and  the  most  useful  filter  tuning  tool 
is  provided  by  a  covariance  analysis  (Ref  31) . 

In  the  covariance  analysis  the  covariance  of  the 
estimation  errors  of  the  Kalman  filter  when  applied  to  the 
system  truth  model  is  propagated  forward  in  time  from  the 
initial  conditions  on  the  covariances  of  the  truth  model 
states.  In  parallel  with  this  estimate-error  covariance 
propagation,  the  filter  covariance  is  itself  propagated 
forward  in  time.  The  true  and  filter  computed  estimation 
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error  covariances  may  then  be  compared  since  the  truth 
model  state  estimation  error  covariances  can  be  transformed 
to  errors  for  the  design  states  using  equation  (3-21) . 

The  designer  may  then  modify  the  dynamics  noise  and  mea¬ 
surement  noise  strengths  to  achieve  the  desired  filter 
performance:  the  duration  of  the  initial  estimation 

transient  and  the  steady  state  covariance  obtained.  The 

0 

development  of  Reference  31  is  summarized  here  to  present  the 
necessary  equations  for  the  covariance  analysis. 


Define  the  augmented  state  vector  x  by 


*c<t  )  = 

1  x(t.) 


(4-9) 


where  xt  are  the  truth  model  states  and  x  are  the  filter 
state  estimates.  For  the  filter  implemented  without  impul¬ 
sive  feedback  (Ref  31) ,  the  augmented  state  vector  satis¬ 
fies  the  time  propagation  relation 


Sc'V  *  2c5c(ti-l)  +  Sca(ti-l' 


(4-10) 


where 
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and  w„  is  zero-mean  white  Gaussian  discrete-time  noise  of 
^d 

discrete-time  noise  covariance 


with 


"V*! ,^d,tj,)  ' 


(4-llb) 


and 


(4-llc) 


(4-lld) 


In  these  equations,  J>  and  are  from  equations  (3-23a) 
and  (3-13a)  respectively,  while  Qt  and  are  from  equations 
(3-18)  and  (3-17)  respectively. 

The  measurement  update  can  be  represented  by 


5c(tt’  - 


A  x 


<t.) 


in  which 


(4-12) 


A_  = 


Mt 


I-K„H 
- F— a 


(4-13a) 


(4-13b) 


and  K„  is  the  steady  state  Kalman  filter  gain,  H  is  as 
defined  in  equation  (3-105) ,  and  Ht  and  vfc  are  as  defined 
in  equation  (3-19)  for  the  truth  model.  Initial  conditions 
are  provided  by 


¥c(to' 


xt(t0) 


(4-13c) 


The  covariance  of  the  augmented  truth  state  and 
estimation  error  is  propagated  forward  from  the  initial 
covariance 


by 


and 


Et<V 


Sc'V 


(4-14a) 


(4-14b) 


Pc(<)  = 


A_P. 


(t 7)  +  K  R  kT 

1  — c  — c— t— c 


(4-14c) 


for  which  Pt(tQ)  is  the  initial  covariance  matrix  for  the 
truth  model  state  initial  conditions  and  Rt  is  the  strength 
of  the  discrete-time  measurement  noise  for  the  truth  model 
(equation  3-20) ) . 

Since  it  will  be  desired  ultimately  to  consider  the 
estimation  error  of  the  design  model's  system  and  dis¬ 
turbance  states,  it  is  appropriate  to  define  the  filter 
estimation  error  as 


sc(ti)  - 


(4-15) 


with 


Sc  -  t-ct  !  I] 


(4-16a) 
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where  TDT  and  T^  are  defined  in  equation  (3-21)  . 

At  each  sample  period,  the  estimation  error  covari¬ 
ance  is  thus 

=  SA(t^  (4~17) 

and  whose  diagonal  elements  are  the  variances  of  the  esti¬ 
mation  error  for  each  system  and  disturbance  state. 

The  Kalman  filter  was  designed  as  a  constant-gain 
filter  according  to  the  assumption  that  its  computed  ’ 
covariance  attains  a  constant  value  in  steady-state.  The 
filter's  computed  covariance  is  given  by  equation  (3-109) 

and  is  denoted  P_ . 

—a 

Finally,  define  the  vectors  of  standard  deviations 
of  the  true  estimation-error  and  the  filter's  computed 
estimation-error  for  the  design  model  system  and  dis¬ 
turbance  states  as  s  and  s  ,  respectively.  The  jth 

C 

element  of  the  time -varying  vector  s  is, 

C 


while  the  jth  element  of  the  constant  vector  Sp  is, 

i - 


(4-18a) 


Plots  of  both  s  and  s  for  each  state  with 

°i  pj 

respect  to  time  gives  information  about  the  Kalman  filter 
performance  ideally  suited  to  the  needs  of  achieving  a 
well-tuned  design.  They  demonstrate  both  the  filter's 
transient  and  steady  state  response  for  each  state  esti¬ 
mated.  Appropriate  changes  in  the  state  and  measurement 
noise  strengths  (£,  Qn,  and  R  of  equations  (3-2),  (3-4), 
and  (3-7) )  can  then  be  made  and  the  corresponding  perform¬ 
ance  re-evaluated  (Refs  2,  31) . 
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V.  CGT/PI/KF  Design  Computer  Program 

5.1  Introduction 

The  primary  objective  of  this  thesis  effort  is  to 
create  a  computer  program  with  which  to  design  CGT/PI/KF 
controllers.  This  chapter  presents  a  general  description 
of  the  program  \diich  has  been  developed — hereafter  to  be 
referred  to  as  "CGTPIF."  More  complete  descriptions  of 
CGTPIF  are  given  in  various  appendices  to  this  thesis. 

While  the  specific  test  application  for  the  pro¬ 
gram  in  the  context  of  this  thesis  has  been  related  to 
aircraft  control  design,  CGTPIF  is  written  to  be  appli¬ 
cable  to  a  wide  variety  of  control  design  problems.  It 
has  the  following  attributes: 

1.  CGTPIF  executes  interactively 

2.  The  program  utilizes  efficient  array  allocation 

a.  Initial  memory  allocation  easily  set 

b.  Dynamic  array  allocation  within  total 
memory  allocated 

3.  Various  modes  of  entry  are  possible  for  the 
dynamics  models 

4.  Design  paths  are  automatically  followed,  with 
user  prompts  at  necessary  decision  points 

5.  Requests  for  input  include  informative  prompts 


6.  Copious  program  output  is  provided 

a.  Output  roost  relevant  to  design  decisions 
are  provided  directly  to  the  terminal 

b.  Additional  detailed  output  is  provided  to 
a  separate  output  file 

7.  Information  relevant  to  design  iteration  is 
preserved 

8.  Error  checking  is  performed,  and  messages  given 
as  appropriate 

CGTPIF  employs  computational  routines  available 
in  a  library  of  matrix  computer  routines  described  in 
Reference  24.  Exclusive  of  the  library  routines,  the 
program  has  a  length  of  about  2500  lines  of  source  code. 

The  programming  language  employed  is  ANSI  Standard 
FORTRAN  IV.  Although  the  resulting  source  code  is  highly 
portable,  local  memory  utilization  limits  for  interactive 
execution  may  impose  constraints.  In  use  on  a  Control  Data 
Corporation  CYBER  machine,  the  necessary  load  size  was 
achievable  with  no  impact  on  the  source  code.  Thus  the 
existing  source  code  is  in  a  pure  form  for  whatever  system/ 
implementation  motivated  modifications  may  be  required  to 
achieve  interactive  load  size  limits. 

5.2  Program  Operating  Principles 

and  Organization 

CGTPIF  has  three  design  paths:  (1)  design  of  a  PI 
regulator;  (2)  design  of  an  open-loop  CGT  or  closed-loop 
CGT/PI  controller,  and  (3)  design  of  a  Kalman  filter. 
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Corresponding  to  the  first  two  design  options  is  a  con¬ 
troller  evaluation  set  of  routines,  and  corresponding  to 
the  third  design  option  is  a  set  of  routines  for  filter 
evaluation.  The  evaluation  routines  perform  the  computa¬ 
tions  discussed  in  Chapter  IV. 

A  general  flowchart  of  CGTPIF  is  given  in  Figure 
5-1,  showing  the  main  execution  paths  and  design  entries. 
The  controller  sample  period  is  entered,  the  design  model 
is  established,  and  then  the  desired  design  path  can  be 
followed.  The  elementary  design  path  choice  is  between 
controller  and  filter  designs.  The  CGT  or  CGT/PI  and  the 
PI  regulator  design  paths  are  then  options  within  the 
controller  design.  If  the  PI  controller  design  is  pursued 
prior  to  the  CGT  design  path,  the  CGT  design  will  auto¬ 
matically  be  of  the  CGT/PI  controller.  If  the  PI  is  not 
already  determined  during  the  current  execution  of  the 
program  then  the  designer  may  elect  to  design  either  a 
CGT  or  CGT/PI  controller.  However,  the  CGT  controller 
design  is  not  pursued  if  the  open-loop  design  model  is 
unstable.  The  controller  design  path  is  followed  auto¬ 
matically  by  the  appropriate  controller  evaluation  path. 
Similarly,  the  filter  design  path  leads  automatically  to 
the  filter  evaluation.  When  the  evaluation  is  complete  the 
designer  is  given  the  opportunity  to  loop  on  the  design 
path,  choose  a  different  design  path,  or  terminate  program 
execution.  More  detailed  flowcharts  and  functional  dia¬ 
grams  are  given  in  the  appendices. 
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5.2.1  Interactive  Execution .  Execution  of  these 
options  entails  both  a  great  deal  of  program  code  and  ] 

memory  usage  for  array  storage.  The  inherent  memory 
requirements  are  so  large  that,  despite  careful  coding  and 
efficient  array  storage  techniques,  the  final  program  size  j 

was  much  greater  than  the  65000  octal  word  limitation  of  < 

! 

the  CYBER  interactive  system.  Thus,  in  order  to  achieve 
interactive  operation  and  provide  sufficient  free  memory 
for  array  allocation  so  that  problems  of  large  and  vari¬ 
able  dimensions  could  be  treated,  a  CYBER  Loader  option 
referred  to  as  "Segmentation"  (Ref  13)  was  employed. 

While  more  detailed  information  about  segmentation 
is  presented  in  the  appendices,  a  few  general  comments 
about  it  are  appropriate  here.  An  attractive  feature  of 
segmentation  distinguishing  it  from  other  methods  for  ( 

achieving  selective  loading  of  program  elements  (such  as 
"Overlays") ,  is  that  no  modification  to  the  source  code 
is  involved  in  achieving  segmentation.  A  segmented 
executable  object  file  is  produced  from  a  job  run  which 
executes  job  control  and  segmentation  directives  operating 
on  the  program's  compiled  object  file  and  any  needed  object 
libraries.  The  job  control  sequence  and  segmentation  direc¬ 
tives  are  invariant  for  the  kinds  of  program  modifications 
which  may  be  required  to  apply  CGTPIF  to  specific  design 
problems.  The  resulting  program  executable  file  then  may 
be  run  interactively  and  all  program  loading  and  unloading 
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occurs  dynamically  in  the  course  of  execution  and  is  fully 
transparent  to  the  user. 

CGTPIF  was  written  specifically  to  run  inter¬ 
actively.  Requests  for  input,  while  intentionally  brief, 
tell  the  user  what  is  expected  of  him — what,  how  many, 
and  in  what  units,  as  appropriate.  Output  of  information 
relevant  to  design  decisions  to  the  terminal  is  compact 
and  automatically  provided.  Also,  the  user  can  determine 
the  amount  and  category  of  output  to  the  terminal  in  some 
cases,  and  according  to  need.  These  and  other  characteris¬ 
tics  pertaining  to  CGTPIF1 s  interactive  operation  are  dis¬ 
cussed  in  more  detail  in  succeeding  subsections  of  this 
chapter . 


5.2.2  Array  Allocation.  It  is  intended  that 
CGTPIF  will  not  require  modification  in  order  to  be  applied 
to  specific  design  problems.  Named  Common  blocks  are  used 
corresponding  to  various  computational  elements  of  the 
program.  For  each  Common  block  an  equation  is  given 
in  Appendix  A  to  allow  computation  of  the  minimum  total 
memory  allocation  needed  as  a  function  of  the  dimensions 
of  the  various  dynamics  models  described  in  Chapter  III: 
the  "design,"  "truth,"  and  "command"  models.  As  currently 
written,  CGTPIF  handles  problems  of  orders  15,  20,  and  10 
for  the  design,  truth,  and  command  models,  respectively. 

In  many  cases  the  existing  allocations  will  be  satisfac¬ 
tory  since  any  actual  usage  less  than  that  already 
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allocated  for  will  execute  properly.  If  necessary ,  alloca¬ 
tions  can  be  computed  according  to  the  above  mentioned 
equations  defining  storage  requirements.  The  changed 
allocations  need  only  be  set  in  the  "Main"  program  routine. 
All  other  routines  adapt  accordingly  with  no  need  for 
modification  to  the  source  code. 

For  each  Common  block  a  single  array  is  allocated. 
The  various  subroutines  of  CGTPIF  then  partition  the  given 
allocation  as  needed  to  accommodate  the  specific  arrays 
employed  in  execution.  Each  array  is  thus  stored  in  the 
minimum  memory  needed  to  contain  it  fully.  This  method  of 
storage  is  in  contrast  to  the  common  technique,  wherein 
each  array  is  individually  and  explicitly  allocated.  This 
dynamic  allocation  technique  instead  allows  full  use  of  the 
total  allocation  for  each  Common  block  as  needed  according 
to  the  specific  set  of  dimensionalities  employed,  and  takes 
no  more  storage  for  any  individual  array  than  absolutely 
required  to  contain  it.  Thus,  for  any  specific  set  of 
dimensions,  memory  usage  is  minimized,  and  a  large  variety 
of  different  dimensionalities  can  be  accommodated  within 
the  same  total  allocation. 

5.2.3  Entry  of  Dynamic s  Models.  Any  of  the  three 
dynamics  models  described  in  Chapter  III  may  be  established 
in  any  of  three  ways:  (1)  the  dimensions  and  matrix  ele¬ 
ments  may  be  entered  directly  from  the  terminal  with  input 
prompting  from  CGTPIF;  (2)  the  dimensions  and  matrix 


elements  may  be  entered  from  a  "DATA"  file  on  which  the 
dynamics  model  from  a  previous  run  of  the  program  was 
written  (the  writing  of  such  a  file  entry  for  each  model 
is  a  program  option) ;  or  (3)  the  dimensions  and  matrix 
elements  may  be  established  by  user  provided  subroutines. 

When  entered  from  the  terminal  directly  or  set  by 
subroutines,  only  non-zero  elements  of  the  various  matrices 
need  be  given.  In  many  cases,  this  substantially  simpli¬ 
fies  establishing  the  matrices  and  reduces  the  probability 
of  erroneous  entries.  Subsequent  design  runs  for  the  same 
problem  can  then  simply  read  the  models  from  the  "DATA” 
file  previously  created. 

5.2.4  Design  Paths .  Rather  than  require  the  user 
to  specify  step-by-step  the  computations  to  be  performed 
in  the  design  process,  the  program  follows  paths  auto¬ 
matically  and  gives  prompts  at  points  where  options, 
changes,  or  inputs  are  required.  Thus,  no  elaborate  or 
possibly  coded  list  of  directives  is  needed  to  execute  a 
design. 

Each  design  sequence  is  followed  by  entry  into  the 
appropriate  design  evaluation  section  of  code.  After 
evaluation  of  the  specific  design,  any  of  the  design  paths 
may  again  be  selected,  or  program  execution  terminated. 

5.2.5  Input  Prompts.  Each  request  for  input 
from  the  user  includes  a  prompt  by  CGTPIF  of  information 
defining  the  input  desired.  When  an  option  is  offered,  the 
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prompt  succinctly  describes  the  option  and  specifies  the 
appropriate  form  of  response  (e.g.,"(Y  or  N) "  signifying 
yes  or  no,  respectively) .  When  a  dimension  is  to  be 
entered,  the  dimension  is  identified.  When  a  matrix  is  to 
be  entered,  the  iratrix  is  identified,  as  is  its  dimension¬ 
ality,  and  the  appropriate  form  of  input.  Similarly,  all 
other  requests  for  input  identify  the  input  expected  and 
the  form  the  entry  is  to  take. 

However,  since  writing  extensive  prompts  to  the 
terminal  would  impede  streamlined  execution  of  the  program, 
the  messages  are  brief  and  require  that  the  user  have  some 
understanding  of  what  is  involved  in  achieving  each  of  the 
designs.  The  prompts  are  intended  to  assist  users  familiar 
with  the  elements  of  the  PI,  CGT,  CGT/PI,  and  KF  design 
methods  and  with  the  terminology  used  in  this  thesis  to 
enter  the  necessary  information  for  such  designs  into  the 
program. 

5.2.6  Program  Output .  In  the  computations 
involved  in  the  various  design  paths,  a  great  deal  of  infor¬ 
mation  is  generated.  While  all  the  information  generated 
may  be  relevant  to  the  design,  in  the  usual  case  only  a 
small  fraction  of  the  information  is  needed  to  pursue  itera¬ 
tion  of  the  design  paths. 

Information  most  relevant  to  achieving  the  various 
designs  is  output  to  the  terminal  either  automatically  or 
at  the  user's  option.  The  same  information,  along  with  all 
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other  potentially  useful  information  generated  by  the  pro¬ 
gram,  is  also  output  to  a  "LIST"  file.  After  program  exe¬ 
cution  is  complete,  the  user  may  then  look  through  this 
file  at  the  terminal  (as  with  "PAGE"  command,  for  example) 
or  may  "ROUTE"  it  to  a  line  printer  for  a  complete  listing. 
All  output  is  given  an  identifying  name  or  description 
which  is  consistent  with  the  terminology  of  this  thesis. 

Plots  are  available  as  options  of  the  design  evalua 
tion  routines.  The  plots  produced  are  of  the  "line-printer 
type.  For  the  controller,  plots  of  selected  variables  may 
be  output  to  the  terminal,  in  addition  to  the  full  set  of 
plots  output  to  the  "LIST"  file. 

5.2.7  Preservation  of  Design  Information .  As  dis¬ 
cussed  previously,  the  dynamics  models,  once  defined,  may 
be  written  to  a  file  from  which  they  may  be  extracted 
during  subsequent  runs  of  the  program.  In  addition,  the 
feedback  gains  of  the  PI  regulator  may  be  written  to  the 
same  file  for  later  reuse. 

During  program  execution,  user  entries  in  the 
design  iteration  are  preserved.  Thus,  for  each  iteration 
only  those  entries  to  be  changed  need  be  given  as  input, 
making  design  iteration  both  fast  and  easy  in  terms  of  the 
simple  mechanics  of  the  process.  Also,  computations  which 
are  not  modified  by  design  iteration  are  performed  only 
once  and  the  results  are  preserved  for  reuse  within  the 
current  run,  as  needed. 
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5.2.8  Error  Checking.  A  common  problem  encountered 
in  executing  computer  programs  with  variable  dimension  array 
storage  is  the  tin  intentional  (and  often  unknown)  over¬ 
running  of  the  allocated  storage  area.  When  this  happens, 
the  program  may  fail  (due  to  the  over-writing  of  program 
code,  for  example)  or,  even  worse,  the  program  may  appear 
to  run  properly  but  provide  erroneous  results.  To  avoid 
these  difficulties,  before  using  each  of  the  Common  arrays, 
OGTPIF  computes  the  allocation  needed  for  the  arrays  it 
will  generate  and  compares  it  with  the  number  of  words  of 
memory  actually  allocated.  If  more  memory  is  needed  than 
has  been  provided,  a  message  is  written  indicating  the 
problem,  the  Common  in  question,  and  the  minimum  allocation 
needed.  Execution  is  then  aborted.  As  described  in 
Appendix  A,  allocations  may  be  changed  in  the  Main  routine 
and  a  new  program  created  to  achieve  the  necessary  array 
storage . 

Error  checks  for  array  entry  from  the  terminal  are 
also  performed.  Not  only  are  the  array  dimensions  identi¬ 
fied  in  prompts  requesting  entry  of  arrays,  but  each  entry 
is  checked  to  verify  that  it  is  within  the  row,  column 
bounds  for  the  array.  If  not  within  bounds,  the  entry  is 
not  accepted  and  a  message  identifying  the  problem  is  given. 
Also,  for  matrices  requiring  special  properties  (e.g., 
positive  semi-def initeness) ,  entries  which  clearly  violate 
these  requirements  (e.g.,  negative  diagonal  elements)  are 
not  accepted  and  a  message  is  written  to  the  terminal. 
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Various  other  error  checks  are  performed  for  each 
input  entered  from  the  terminal.  These  checks  ensure  that 
no  invalid  entries  are  accepted  (e.g.,  non-positive  con¬ 
troller  sample  time) ,  and  no  storage  out  of  array  bounds 
occurs.  Obviously  it  is  not  possible  to  guard  against 
valid  yet  erroneous  entries ,  but  in  most  cases  the  program 
provides  opportunities  to  correct  mistaken  entries  made 
and  discovered  by  the  user. 

Additional  tests  are  performed  to  ensure  dimen¬ 
sional  consistency  of  the  dynamic  models.  These  conditions 
are  discussed  in  Chapter  III  as  each  model  description  is 
defined.  For  example,  CGTPIF  checks  that  the  numbers  of 
inputs  and  outputs  defined  for  the  design  model  are  equal, 
as  well  as  checking  each  of  the  other  dimensional  condi¬ 
tions  as  each  model  is  established.  If  the  condition  is 
not  met,  a  message  identifying  the  problem  is  written  and 
the  program  execution  is  aborted. 

5.3  Program  Usage 

CGTPIF  was  written  with  the  intention  of  providing 
a  CGT/PI/KF  controller  design  program  that  is  efficient 
in  memory  and  execution  time,  streamlined  for  the  user  in 
input  and  output  requirements  in  interactive  operation, 
and  applicable  to  a  variety  of  problems.  These  objectives 
have  been  achieved. 

Preparation  for  use  of  the  program  consists  pri¬ 
marily  in  determining  the  dimensions  of  the  models  to  be 
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employed  and  the  specific  coefficients  of  the  matrices 
comprising  those  models.  It  is  appropriate  for  the  PI 
design  that  an  initial  set  of  quadratic  weights  be  estab¬ 
lished  also  in  order  to  begin  the  design  iteration  process. 

Assuming  that  the  user  understands  the  nature  and 
requirements  of  the  various  design  elements,  and  is 
familiar  with  terminology  related  to  such  designs,  the 
actual  execution  of  the  program  involves  straightforward 
response  to  input  prompts.  The  specific  meaning  of  all 
program  prompts  employed  is  delineated  in  Appendix  B  of 
this  thesis.  Information  most  useful  to  the  design  itera¬ 
tions  is  available  directly  as  output  to  the  terminal, 
while  additional  information  is  provided  in  a  separate 
output  file.  The  terminology  employed  in  the  output  is 
defined  in  Appendix  B  and  related  to  specific  program  com¬ 
putations  in  Appendix  A  of  this  thesis. 

Results  from  the  use  of  CGTPIF  are  presented  in 
the  next  chapter.  More  detailed  information  in  the 
appendices  (A  through  E)  include:  a  "Programmer's  Guide," 
a  "User's  Guide,"  sample  input  and  output,  a  program  list¬ 
ing,  and  a  listing  of  the  job  card  sequence  needed  to 
obtain  a  segmented  executable  object  file  of  CGTPIF. 
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VI.  CGT  Design  Results 


6 .1  Introduction 

The  program  developed  in  this  study  has  been 
applied  to  a  variety  of  control  design  problems  in  order 
to  verify  its  proper  functioning  and  also  to  evaluate  the 
design  process  and  the  characteristics  of  the  controller 
designs  achievable.  This  chapter  discusses  several 
CGT/PI  and  CGT/PI/KF  designs.  The  purpose  is  to  demon¬ 
strate  the  systematic,  logical  design  process  and  to  show 
the  capabilities  achieved  by  the  designs.  Although  the 
design  iterations  are  not  discussed  in  detail,  a  summary 
of  the  iterative  process  is  given  in  Section  6.1.2  below. 

6.1.1  Design  Examples.  The  designs  discussed  here 
are  for  two  different  systems.  A  simple  design  example 
uses  a  lightly-damped,  second-order,  single-input-single- 
output  (SISO)  system  model,  with  two  constant  disturbances 
also  driving  the  dynamics.  The  second  system  is  an  air¬ 
craft  longitudinal  dynamics  model  with  unstable  dynamics 
and  with  two  control  inputs.  For  both  systems,  CGT  designs 
based  upon  two  different  command  model  dynamics  descrip¬ 
tions  are  developed. 

For  the  simple  second-order  system,  open-loop  CGT 
and  closed-loop  CGT/PI  designs  are  developed.  Command 
model  dynamics  representing  a  first-order  system  and  a 
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well-damped  second-order  system  are  used.  Specifically, 
the  following  designs  are  discussed:  (1)  design  of  open- 
loop  CGT  and  closed-loop  CGT/PI  controllers  using  the 
first-order  command  model  without  disturbances  to  demon¬ 
strate  the  improvement  in  model-following  achievable  with 
the  closed-loop  design;  (2)  two  different  representations 
of  the  second-order  command  model  are  used  to  demonstrate 
the  invariance  of  controller  response  with  alternative 
representations  of  command  model  dynamics;  (3)  open-loop 
CGT  and  closed-loop  CGT/PI  designs  are  developed  for  the 
second-order  command  model  dynamics;  (4)  these  are  evalu¬ 
ated  for  the  cases  of  modeled  constant  disturbances  of  zero 
or  non- zero  unknown  magnitude  acting  on  the  system  dynamics 
to  demonstrate  the  capabilities  of  the  open-loop  implementa¬ 
tion  and  the  improvements  possible  with  the  closed-loop 
implementation;  and  (5)  these  same  controllers  are  then 
applied  to  a  second-order  system  of  different  dynamics 
than  that  used  for  design  to  show  the  affects  of  modeling 
error  on  the  performance  of  the  CGT  controllers. 

For  the  aircraft  longitudinal  dynamics  system, 
three  different  design  models  are  used.  All  three  design 
models  employ  simple  three-state  dynamics  models,  while  one 
includes  actuator  states  also,  and  the  third  includes  both 
actuator  and  clear-air  turbulence  states.  These  are 
derived  as  approximations  of  a  truth  model  which  includes 
a  four-state  dynamics  model  and  actuator  and  clear-air 
turbulence  states.  Several  of  these  simplified  models  then 
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ar- 


serve  as  design  and  truth  models  for  specific  design  cases. 
The  following  CGT/PI  designs  are  discussed:  (1)  design  of 
decoupled  pitch  controllers  for  all  three  design  models 
demonstrate  the  performance  of  the  CGT/PI  controller  in 
achieving  decoupled  response  characteristics  and  the 
effects  of  inclusion  or  non-inclusion  of  actuator  dynamics 
in  the  design  model;  (2)  the  decoupling  CGT/PI  controller 
for  the  simplest  design  model  is  evaluated  with  respect 
to  truth  models  with  different  parameter  values  to  demon¬ 
strate  controller  performance  when  subject  to  modeling 
errors  or  parameter  variation;  (3)  design  of  conventional 
pitch  controllers  for  the  two  simplest  design  models  to 
demonstrate  the  ability  of  the  CGT/PI  design  in  the  case 
in  which  fewer  system  outputs  are  to  be  of  constrained 
dynamics  than  there  are  independent  controls  available; 
and  (4)  as  in  item  (2)  immediately  above,  the  conventional 
pitch  controllers  for  the  simplest  design  model  are  evalu¬ 
ated  when  subject  to  model  errors  and  parameter  variation. 
Finally,  a  single  Kalman  filter  design  and  covariance 
analysis  is  discussed.  The  filter  design  is  based  upon 
the  design  model  which  includes  turbulence  states  and  is 
evaluated  with  respect  to  the  system  truth  model. 

All  controller  and  filter  designs  are  based  upon  a 
sample  time  of  0.02  seconds.  Such  a  sample  period  is  repre¬ 
sentative  of  the  controller  sampling  times  currently 
employed  in  digital  flight  control  systems.  Since  a 
direct-digital  design  is  effected,  the  controller  sample 
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period  may  be  set  at  any  value  appropriate  to  the  design 
problem  under  consideration.  This  is  in  contrast  to  digi¬ 
tal  designs  based  upon  discrete  approximations  of  continu¬ 
ous  controller  designs  for  which  the  sample  period  is  con¬ 
strained  by  the  validity  of  the  approximations  employed. 

In  order  to  facilitate  reference  to  gains  employed 
by  the  feedforward  CGT  controller,  the  following  defini¬ 
tions  of  matrices  developed  in  Chapter  III  are  given: 
feedforward  gain  matrices  applying  to  the  command  model 
states,  command  model  inputs,  and  design  model  disturbance 

states  are  referred  to  as  K  ,  K  ,  and  K  respectively. 

"*u  “*n 

Thus  for  the  open-loop  CGT  these  matrices  are,  from  equa¬ 
tion  (3-52) 


and 


*x  -*21 
m 


*xu  *22 


-xn  “  *23 
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(6-lb) 
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For  the  closed-loop  CGT/PI  the  feedforward  gain  matrices 
are,  from  equation  (3-116) 


and 
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Similarly,  gains  of  the  PI  regulator  are  referred 
to  repeatedly  in  this  chapter.  The  feedback  gain  matrices 


are  as  given  by  equations  (3-90a)  and  (3-90b)  for  Kx  gains 
on  system  states  and  K  gains  on  system  output  errors. 

The  gain  matrix  for  the  simple  regulator  alone  is  G*  as 

C 

given  in  equation  (3-76). 

6.1.2  PI,  Regulator  Design  for  CGT/PI  Application. 
The  design  of  PI  regulators  in  the  context  of  the  CGT/PI 
controller  proves  to  be  straightforward,  and  is  not  dis¬ 
cussed  in  detail  for  the  specific  design  examples.  The 
PI  regulator  is  used  to  null  errors  between  the  reference 
feedforward  inputs  provided  by  the  CGT  and  the  system's 
true  state  and  output  response.  The  PI  design  objective 
for  the  inner-loop  implementation  is  primarily  motivated 
by  the  need  to  achieve  rapid,  well-damped  zeroing  of  these 
errors.  The  actual  control  magnitudes  and  rates  required 
to  achieve  such  response  need  not  be  constrained  in  the  PI 
design  itself  to  design  limit  values.  The  CGT/PI  controller 
generates  input  magnitudes  and  rates  that  often  are  more  a 
function  of  the  CGT  feedforward  solutions  than  of  the  PI 
feedback  commands  alone.  Typically,  in  the  closed-loop 
configuration  the  most  severe  input  magnitude  and  rate 
commands  are  due  to  the  feedforward  gains  on  the  command 
model  inputs.  As  can  be  seen  in  equation  (6 -2b) ,  the  feed¬ 
forward  gain  K  on  the  command  model  inputs  depends  on  the 
xu 

PI  feedback  gain  K  (defined  in  equation  3-90a))  and  also 


on  the  partitions  and  Ajj  of  the  CGT  equations  (given 
by  equations  (3-49b)  and  (3-49e)). 

The  procedure  which  was  found  useful  in  the  PI- 
CGT/PI  design  iteration  is  described  by  the  following  steps 

1.  Determine  initial  estimates  of  appropriate 
quadratic  weights  (for  the  cost  formulation  of  equation 
(3-61))  by  selecting  a  set  of  output  errors  and  correspond¬ 
ing  input  magnitude  and  rate  maximum  admissible  values  to 
null  the  output  errors.  Output  weights  are  then  chosen  as 
the  reciprocal  of  the  square  of  each  output  error.  Simi¬ 
larly,  the  input  magnitude  and  rate  weightings  are  chosen 
as  the  reciprocals  of  the  squares  of  the  maximum  magnitude 
and  rate  input  values.  These  weightings  may  be  scaled  to 
be  symmetrical  about  unity,  with  smallest  value  unity,  or 
however  desired  for  convenience  since  only  the  relative 
values  of  the  weights  affect  the  design  results. 

2.  Using  the  quadratic  weights  determined  in 
step  1,  a  PI  regulator  is  designed  and  the  system  response 
to  state  initial  conditions  is  evaluated. 

3.  Often  the  PI  gains  are  significantly  greater 
than  the  gains  of  the  simple  regulator  for  the  same  problem 
and  the  speed  and  damping  of  the  response  is  inferior  for 
the  given  input  maximum  admissible  values.  But  since  the 
inputs  generated  in  command  following  with  the  CGT/PI  con¬ 
troller  depend  in  large  measure  on  the  feedforward  gains, 
initial  iteration  of  the  PI  design  is  directed  at  achieving 
good  speed  and  deunping  of  the  regulator  response  despite 


possibly  exceeding  the  maximum  admissible  input  values. 

Thus  the  weights  on  output  errors  are  increased  and  the  out¬ 
put  behavior  in  response  to  initial  conditions  is  evaluated 
iteratively.  At  the  same  time  it  is  often  useful  to 
include  additional  weights  on  the  first  derivatives  of  the 
outputs  to  improve  the  damping.  Approximate  output  deriva¬ 
tive  weighting  can  be  achieved  easily  in  many  cases  by 
introducing  weights  on  states  included  in  the  output  rate 
equations.  Although  the  input  magnitudes  and  rates  are 
not  specifically  constrained  by  design  limits  at  this  stage 
of  the  design,  it  is  necessary  that  judicious  adjustment 
of  the  relevant  quadratic  weights  be  pursued  with  the 
objective  of  achieving  the  desired  output  response  with 
minimum  control  power.  Thus  it  is  not  appropriate  to 
increase  output  weightings  by  orders  of  magnitude  arbi¬ 
trarily.  The  designer  must  instead  selectively  apply 
quadratic  weights  on  states  and  outputs  that  optimize 
output  regulation  with  respect  to  control  power. 

4.  When  good  regulation  has  been  achieved  by 
iteration  on  the  PI  design  alone,  compute  the  CGT/PI  con¬ 
troller  and  evaluate  its  response  to  command  inputs. 
Difficulties  with  input  magnitudes  and  rates  may  then  be 
apparent,  and  will  be  most  severe  during  the  first  few  con¬ 
troller  sample  periods  after  application  of  the  command  due 
to  direct  feedforward  of  the  command  input. 

5.  If  necessary,  iterative  design  of  the  CGT/PI 
controller  then  proceeds  by  computing  new  PI  gains  to  give 
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improved  output  regulation  and/or  reduced  control  inputs, 
followed  by  recomputation  of  the  CGT/PI  gains  and  reevalua¬ 
tion  of  the  system  response  to  command  inputs.  The  goal  is 
to  achieve  a  final  design  which  provides  good  model- 
following  without  exceeding  input  constraints.  This  may 
not  be  achievable  in  all  cases:  sometimes  the  CGT  feed¬ 
forward  gain  on  the  command  input  is  large  in  itself,  and 
the  contribution  of  the  PI  gain  to  it  cannot  achieve  accept¬ 
able  final  feedforward  gain  (see  equation  (6-2b)). 

6.1.3  Determining  Command  Models.  For  single- 
input-single-output  (S1S0)  system  command  models  for  spe¬ 
cified  dynamics  are  readily  determined  (e.g.,  from  a  trans¬ 
fer  function  defining  the  desired  response,  a  state  model 
may  be  obtained  directly  in  the  "standard  controllable 
form"  (Ref  31)).  Moreover,  any  state  description  yielding 
the  specified  command  model  dynamics  will  yield  identical 
response  to  that  of  any  other  state  description  with  the 
same  input/output  characteristics,  for  either  the  open-loop 
CGT  or  the  closed-loop  CGT/PI  controller.  This  is  demon¬ 
strated  in  cases  to  follow. 

For  multi-input-multi-output  (MIMO)  systems, 
determining  appropriate  command  models  can  be  more  diffi¬ 
cult.  In  some  cases  a  command  model  can  be  obtained  from 
a  generalization  of  the  standard  controllable  form  which 
can  be  written  for  a  matrix  of  transfer  functions  (Ref  19) . 
If  the  desired  output  response  is  defined  as  decoupled  and 
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first-order  between  each  output  and  input,  then  command 
models  are  easily  determined  and  need  include  only  diagonal 
matrices.  Determination  of  such  first-order,  decoupled 
command  models  is  shown  for  the  aircraft  controller  design 
in  Section  6.4.5.  As  in  the  SISO  case,  any  model  which 
achieves  the  desired  dynamics  may  be  used  in  the  controller 
with  equivalent  results. 

The  design  examples  discussed  here  demonstrate  the 
determination  of  command  models.  For  the  flight  control 
example  (M1M0)  only  first-order  decoupled  response  command 
models  are  used. 

6.1.4  Need  for  Complete  Evaluation  Software .  The 
designs  discussed  in  this  chapter  are  evaluated  as  open- 
loop  CGT  or  closed-loop  CGT/PI  controllers  or  as  indepen¬ 
dent  Kalman  filters.  Thus,  the  results  demonstrated  for 
the  CGT/PI  controller  are  valid  strictly  for  the  case  of 
perfect  and  complete  state  knowledge.  Some  degradation  in 
performance,  especially  with  regard  to  robustness  (Ref  32) , 
may  accrue  due  to  inclusion  of  the  state  estimation  pro¬ 
vided  by  a  Kalman  filter.  As  discussed  in  Chapter  IV, 
in  order  to  evaluate  the  final  controller  designs  for 
implementation  as  CGT/PI/KF  controllers  requires  that 
additional  software  be  written.  These  results  provide  the 
basis  for  final  design  iteration  and  demonstrate  the  ulti¬ 
mate  model -following  performance  vdiich  may  be  achieved. 
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6.2  Using  the  CGTPIF  Design  Program 

6.2.1  Introduction.  Many  different  controller 
designs  have  been  developed  using  the  CGTPIF  design  com¬ 
puter  program.  The  program  has  proven  to  be  readily  appli¬ 
cable  to  these  different  problems  despite  the  many  and 
varied  dimensions  needed  to  define  the  dynamics  models 
relevant  to  each.  The  preservation  of  design  information 
both  between  separate  executions  and  within  a  single  execu¬ 
tion  of  the  program  substantially  aids  efficient  design 
iteration.  Also,  the  information  available  directly  at 
the  user's  interactive  computer  terminal  during  program 
execution  has  been  found  appropriate  and  adequate  for 
intelligent  selection  of  design  parameters  necessary  to 
pursue  designs  to  successful  conclusion. 

This  chapter  discusses  results  obtained  by  using 
CGTPIF  for  specific  design  problems,  but  does  not  detail 
the  program  entries  these  designs  entailed.  An  example 
of  the  input  and  output  for  such  a  design  is  included  in 
Appendix  C. 

6.2.2  Operating  Considerations .  A  few  simple 
considerations  of  the  operation  of  the  design  program  par¬ 
ticularly  impact  its  ease  of  application  to  specific  design 
problems.  These  and  other  more  specific  considerations 
are  discussed  more  fully  in  context  in  the  "Programmer's 
Manual"  (Appendix  A)  and  the  "User's  Manual”  (Appendix  B) , 
and  the  reader  is  encouraged  to  refer  to  those  manuals 


104 


before  attempting  to  use  the  program.  The  following  should 
be  noted: 

1.  The  design  model  may  be  entered  only  once  dur¬ 
ing  a  single  program  execution.  Thus,  in  order  to  modify 
the  design  model,  it  is  necessary  to  re-execute  the  pro¬ 
gram. 

2.  The  truth  and  command  models  may  be  redefined 
as  often  as  desired  during  a  single  program  execution. 

This  facilitates  evaluation  of  performance  when  subject 
to  design  model  errors,  and  allows  design  iteration 
involving  changes  to  the  command  model . 

3.  The  "SAVE"  and  "DATA"  files  should  be  used  to 
preserve  the  design,  truth,  and  command  models  between 
distinct  program  executions.  The  two  files  are  different 
according  to  their  local  file  names:  the  "SAVE"  file  is 

a  write-only  file  created  by  CGTPIF  during  program  execu¬ 
tion,  the  "DATA"  file  is  a  read-only  file  which  the  user 
provides  (if  desired)  by  assigning  a  previously  created 
"SAVE"  file  the  local  file  name  "DATA."  Manipulation 
of  local  file  names  does  not  occur  during  program  execu¬ 
tion. 

4.  Models  may  be  read  from  the  "DATA"  file  repeat¬ 
edly  during  program  execution.  Thus,  if  specific  elements 
of  the  matrices  defining  the  command  and/or  truth  models 
are  to  be  modified  during  program  execution,  it  is  most 
convenient  to  retain  nominal  representations  on  a  "DATA" 
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file.  The  model  may  then  be  read  and  modified  readily 
during  subsequent  program  execution. 

5.  Changes  in  model  dimensions  require  that  all 
relevant  defining  matrices  be  entered  anew.  Such  change 
cannot  be  achieved  by  modifying  existing  models  of  differ¬ 
ent  dimension  as  through  augmentation.  Note  also  that 
there  is  some  interdependence  among  the  model  definitions 
implicit  in  the  specific  dimensions  of  each.  Thus, 
intended  changes  in  dimensions  for  a  single  model  may 
require  changes  in  definition  for  another  dynamics  model 
also.  Examples  of  such  implied  constraints  are  equality 
in  number  of  outputs  for  the  design  and  command  models, 
and  of  inputs  and  of  measurements  for  the  design  and  truth 
models.  Also,  because  of  the  need  to  define  the  corres¬ 
pondence  between  the  states  of  the  design  and  truth 
models,  changes  in  the  design  model  may  entail  changes 

in  the  truth  model's  TDT  and  matrices  (equations 
3-21a,b)  . 

6.  Open-loop  CGT  designs  can  only  be  pursued  in 
a  given  execution  of  the  program  if  a  PI  regulator  design 
has  not  previously  been  accomplished  and  PI  gains  have 
not  been  read  from  the  "DATA"  file.  The  program  logic 
assumes  the  usual  design  will  be  of  a  CGT/PI  controller 
and  pursues  such  a  design  automatically  in  the  CGT  design 
path  whenever  PI  gains  exist  within  program  storage. 

7.  Any  of  the  three  design  paths  (PI  regulator, 
CGT  controller,  or  Kalman  filter)  may  be  executed  in  any 


order  and  as  often  as  desired  during  a  single  program 
execution.  Results  from  each  design  path  are  preserved 
independently  throughout  execution. 

8.  In  addition  to  the  output  available  at  the 
user's  terminal,  the  same  and  much  additional  information 
is  output  to  the  file  named  "LIST."  Following  execution, 
it  is  good  practice  to  route  it  to  a  line  printer  for 
listing . 

9.  Plots  of  system  time-response  at  the  user 
terminal  are  the  most  useful  evaluation  tool.  Each  plot 
requires  about  2  minutes  to  be  printed.  In  evaluating  CGT 
controller  response,  it  is  most  convenient  to  include  cor¬ 
responding  system  and  model  outputs  in  the  same  plot  for 
direct  comparison. 

10.  If  it  is  desired  to  include  results  for  all 
of  the  initial  controller  sample  times,  a  value  for  plot 
duration  that  is  less  than  50  times  the  controller  sample 
period  should  be  specified. 

11.  Regardless  of  the  plot  duration  specified,  the 
plot  will  include  50  time  samples  and  will  therefore 
require  a  fixed  amount  of  time  to  be  printed.  These 
samples  are  uniformly  distributed  in  time  over  the 
entire  duration.  The  time  duration  of  the  response  may 

be  adjusted  by  CGTPIF  automatically  to  achieve  50  evenly 
spaced  time  samples  coincident  with  controller  sample  times 

12.  In  order  to  bypass  the  time-response  evaluation 
one  need  only  specify  that  no  plots  are  wished.  For  the 


CGT  evaluation,  specifying  a  zero  index  for  the  command 
input  also  bypasses  the  time-response  evaluation. 

6.2.3  Interpreting  Plots  to  the  Terminal .  The 
plots  of  system  response  produced  by  CGTPIF  are  of  the 
"line-printer"  type.  As  printed  at  the  terminal,  the  time 
axis  runs  vertically  down  the  page  with  the  initial  time 
at  the  top-left  margin  and  time  points  printed  for  each 
sample  down  the  left  margin.  The  dependent  variables  are 
plotted  in  the  horizontal  sense  from  left-to-right  for 
increasing  magnitude.  Each  variable  plotted  is  marked  by 
a  distinct  number  (1  through  5)  at  each  sample  time.  In 
the  event  that  two  variables  occupy  the  same  location  in 
the  plot  field  at  a  sample  time,  only  the  plot  symbol  of 
largest  value  will  be  marked  at  the  point  in  question 
(e.g.,  if  the  variables  represented  by  symbols  1  and  3  both 
quantize  to  the  same  print  location  at  time  tjj,  the  print 
position  will  be  filled  with  symbol  3) .  Such  coincidences 
of  position  can  be  inferred  from  the  behavior  of  the  vari¬ 
able  with  missing  symbol  at  proximate  time  samples.  Note 
that  a  special  case  of  this  is  when  two  or  more  variables 
to  be  plotted  actually  represent  the  same  variable  (as  may 
occur,  for  example,  when  an  output  for  a  system  is  simply 
a  state  and  both  are  plotted) ;  in  such  a  case  only  the  cor¬ 
responding  plot  symbol  of  largest  value  will  appear  in  the 
plot  field. 
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The  horizontal  width  of  the  plot  field  is  50  print 
positions.  The  scale  of  each  plot  is  printed  along  the 
bottom  margin  of  the  plot  and  includes  the  values  of  the 
grid  at  each  multiple  of  10  print  positions.  In  the  usual 
case,  each  variable  is  plotted  on  its  own  unique  scale  and 
thus  a  scale  range  for  each  is  given  corresponding  to  the 
plot  symbol  used.  However,  to  facilitate  comparisons  of 
system  and  model  output  responses,  any  plot  which  includes 
a  model  output  variable  applies  a  single  scale  encompassing 
the  full  range  of  all  variables. 

The  plots  included  in  this  chapter  are  those  pro¬ 
duced  by  CGTPIF  at  the  user  terminal  in  the  course  of 
achieving  the  designs  discussed.  Although  the  plots 
originally  included  identifying  titles,  these  have  been 
removed.  As  presented,  the  plots  are  rotated  180°  from 
their  usual  orientation.  By  rotating  the  bound  thesis  90° 
clockwise  the  independent  axis  (time)  is  horizontal  and  the 
dependent  axis  is  vertical;  the  initial  time  is  at  lower- 
left  and  the  minimum  value  of  each  scale  range  is  along  the 
bottom  grid  line;  scale  identification  is  along  the  right 
vertical  margin. 

6 .3  Simple  Design  Example 

6.3.1  Design  Model.  The  design  model  for  the 
simple  design  example  includes  a  system  state  differential 
equation  in  modified  canonical  form  (Ref  31) 
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I 

<  f 


y  =  [1.  0.]  x 


(6-30 


Note  that  the  disturbances  are  simply  constants.  Responses 
for  the  controlled  system  for  disturbances  acting  or  not 
acting  on  the  system  are  run  by  setting  initial  conditions 
on  the  disturbance  states  as 


nj<°>  - 


or  as 


— d (0)  = 


0, 

0. 

-1. 

2, 


(6-4a) 


(6-4b) 


respectively.  Thus  the  disturbances  are  of  significant 
magnitude  for  the  responses  run  "with  disturbances." 

The  transfer  function  relating  the  output  to  the 
input  is. 


Z  = 
u 


+  4e  +  29 


(6-5) 
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which  has  a  natural  frequency  of  about  5.4  radians  per 
second  and  a  damping  ratio  of  about  0.37.  This  model  has 
a  frequency  and  damping  typical  of  aircraft  "short-period" 
dynamic  s . 

6.3.2  Truth  Model .  The  nominal  truth  model  repre¬ 
sents  the  same  system  dynamics  as  the  design  model, 
incorporating  the  disturbance  states  into  the  truth  model 
state  vector: 


St  = 


-2. 

-5. 

0. 

0. 


5. 

-2. 

0. 

0. 


3. 

1. 

0. 

0. 


-1. 

5. 

0. 

0. 


0. 

1. 

St + 

0. 

0. 

_  — 

u. 


(6-6a) 


The  matrices  relating  the  design  and  truth  model  states 
are,  for  the  design  model's  system  states 


T 

— DT 


1,  0.  0.  0. 

0.  1.  0.  0, 


(6-6b) 


and  for  the  design  model  disturbance  states 


T 


0.  0.  1.  0. 

0.  0.  0.  1. 


(6 -6c) 


An  alternate  dynamics  matrix  (A^)  for  the  truth 
model  is  used  to  evaluate  the  degradation  in  performance 
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due  to  misrepresentation  of  system  dynamics  in  the  design 
model: 


1 


5.  3. 


(6-7) 


0. 


0.  0. 


This  corresponds  to  a  transfer  function  for  system  output 
to  input  of 


Z  =  _ 5 _ 

u  a2  +  2e  +  26 


(6-8) 


which  has  a  natural  frequency  of  about  5.1  radians  per 
second  and  a  damping  ratio  of  about  0.20. 

6.3.3  Command  Models  . 

6. 3. 3.1  First-Order  Command  Model  CM(Ol) . 
For  a  simple  transfer  function  relating  the  desired  out¬ 
put  response  to  a  command  input  given  by 


u  8+5 


m 


(6-9) 


an  appropriate  state  model  is  readily  determined  and  is 
represented  as 

*  *  -5x  +  u  (6-10a) 

m  mm 
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with  output  equation 

ym  -  5xm  (6-10b) 

For  this  model,  the  output  response  is  a  simple  exponen¬ 
tial  with  a  0.2  second  time-constant.  The  command  model 
defined  by  equations  (6-10a,b)  is  hereafter  referred  to  as 
"CM (01) ." 


6.3. 3. 2  Second-Order  Command  Models  CM(02C) , 
CM(02C)  ** .  A  second-order  command  model  with  complex  poles 


of 


8 


1,2 


-5.  ±  j5. 


(6-11) 


and  transfer  function  relating  model  output  response  to 
command  input  of 


50 

2 

8  +  10 8  +  50 


(6-12) 


is  chosen.  This  gives  a  good  representation  of  desirable 
short-period  dynamics,  having  a  natural  frequency  of  about 
7.1  radians  per  second  and  damping  ratio  of  0.707.  State 
differential  equations  are  obtained  corresponding  to  the 
modified  canonical  form  and  standard  controllable  form 
representations.  The  command  models  for  each  form  are 
hereafter  referred  to  as  "CM(02C) "  and  "CM(02C) respec¬ 
tively  . 


L3 


The  state  differential  equation  for  CM(02C)  is, 


Figure  6-1  shows  the  system  output  response  to  a  unit  step 
input  to  the  command  model  with  no  disturbances  acting  on 
the  system,  and  using  the  nominal  truth  model.  Plot 
symbol  1  is  the  system  output  and  plot  symbol  2  is  the 
command  model  output.  Note  that  the  system  response  is 
slow  and  oscillatory  with  a  peak  overshoot  of  about 
14%.  Other  time-response  runs  (not  shown  here)  verified 
that  model -following  was  achieved  in  steady-state.  Clearly 
for  this  controller  a  closed-loop  CGT/PI  is  appropriate 
even  for  the  case  of  no  disturbances  and  perfect  system 
modeling. 

A  PI  regulator  was  designed,  employing  quadratic 
weights  of  100.  on  the  output,  500.  on  state  2,  and  1.  on 
the  control  input  magnitudes  and  rates.  Equal  weights  for 
all  inputs  was  arbitrary  initially  and  a  value  of  unity 
was  used  as  a  scale  basis  for  the  quadratic  weights. 

Weights  on  the  output  and  state  2  approximately  provide 
output  and  output  rate  weighting,  and  the  values  selected 
are  simply  relative  to  those  for  the  inputs.  The  resulting 
feedback  gain  matrix  for  the  simple  regulator  is  (from 
equation  (3-77)) 

G*  =  (-0.1028  0.6440  !  0.4352]  (6-16) 

“—"C 

and  the  PI  regulator  gains  are  (from  equation  (3-90)) 

Kx  =  [16.54  21.41]  (6-17a) 

and  K  =  2.679  (6-17b) 

— z 


115 


Fig.  6-1.  Open-Loop  CGT  for  CM(Ol) ,  (No  Disturbances,  Nominal  Truth  Model) 


Note  that  the  PI  gains  are  much  larger  than  the  gains  of 
the  simple  regulator.  Figure  6-2  shows  the  closed-loop 
response  for  initial  conditions  of  1.  on  both  system  states 
and  zero  on  both  disturbance  states  for  the  nominal  truth 
model.  Plot  symbols  of  1  and  2  signify  the  system  states 
1  and  2,  respectively  (the  system  output  is  state  1) . 

The  response  achieved  is  rapid  and  well-damped. 

The  corresponding  CGT/PI  controller  feedforward 
gains  then  are 


3 

II 

50.91 

(6-18a) 

K 

—v 

20.72 

(6-18b) 

xu 

and 

[-15.05  -.3172] 

(6-18c) 

As  can  be  seen  from  equations  (6-2a,b,c) ,  the  feedforward 
gains  depend  in  part  on  the  feedback  gain  Kx;  in  general, 
increases  in  quadratic  weights  result  in  increased  values 
in  Kx  and  thence  in  the  feedforward  gains. 

The  time  response  for  a  command  input  step  of 
unity,  with  no  disturbance  and  nominal  truth  model  is 
shown  in  Figure  6-3.  Plot  symbols  1  and  2  are  the  system 
and  command  model  outputs,  respectively.  The  resulting 
response  demonstrates  the  benefit  achievable  by  employing 
inner-loop  PI  regulation  to  null  errors  between  the  actual 
system  response  and  the  response  commanded  by  the  CGT 
feedforward  control  laws.  Although  the  model -following  is 
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Fig.  6-2.  PI  Regulator  Response  (Simple  Design 
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CGT/PI  for  CM (01 ) ,  (No  Disturbances,  Nominal  Truth  Model) 


not  perfect,  the  system  output  response  is  clearly  a  good 
representation  of  the  desired  dynamics. 

6.3.5  Design  for  CM(02C)  and  CM(02C) " .  Results 
obtained  using  CM(02C)  are  discussed  in  detail  for  open- 
and  closed-loop  CGT  designs  with  and  without  disturbances 
and  for  nominal  and  alternate  truth  models.  As  expected, 
different  representations  of  given  command  model  dynamics 
yielded  identical  controller  performance — only  the  gains 
on  command  model  states  differed.  Thus,  gains  for  designs 
based  upon  CM(02C)  and  CM(02C) "  are  presented  together, 
but  only  plots  for  CM(02C)  designs  are  shown. 

Gains  for  an  open-loop  CGT  based  upon  CM(02C)  are 


K  =  [20. 

m 

-56.13] 

(6-19a) 

K  =  9.413 
'Xu 

(6-19b) 

and 

K  =  [-2.2 

xn 

-4.6] 

(6-19c) 

For  the  design  based  upon  CM(02C)',  gains  K  and  K  are 


unchanged  while 

the  model 

-xu 

state  gains  become 

~n 

K 

=  [-180.6 

-56.13] 

(6-19a' 

"xm 

Note  that  the  feedforward  gain  for  the  command  model 
CMC02C) "  is  the  product  of  the  gain  of  equation  (6-19a) 
for  CM(02C)  and  the  transformation  matrix  relating  the  two 


state  variable  sets: 


(6-20a) 


x  *  Tx" 
—in  — in 


in  which  x„  and  xl  are  the  state  vectors  of  CM(02C)  and 
— in  -in 

CM(02C)',  respectively  and  the  transformation  matrix  is 


5.  0. 

5.  1. 


and  is  obtained  readily  since  CM(02C) "  is  in  standard 
controllable  form  (Ref  19) .  Relating  the  feedforward  gains 
for  the  two  state  descriptions  leads  to 

K  x  =  K  (Tx')  (6-20c) 

— — m  — x  ' — m 
m  m 

and  finally  gives  the  simple  result  that 

K'  «  Kx  T  (6-20d) 

m  m 

which  is  verified  for  the  values  of  K'  ,  K  ,  and  T  given 

m  m 

by  equations  (6-19a')  ,  (6-19a)  ,  and  (6-20b)  above. 

Figures  6-4  and  6-5  show  the  response  of  the  system 
using  the  open-loop  CGT  controller  for  the  nominal  truth 
model  (equation  (6-6)) ,  without  and  with  disturbances,  respec¬ 
tively.  The  case  of  no  disturbance  but  alternate  truth 
model  (equation  (6-7))  is  given  in  Figure  6-6.  A  unit  step 
on  the  command  input  is  applied  in  all  three  responses. 

Plot  symbol  1  is  the  system  output  and  symbol  2  is  the  com¬ 
mand  model  output. 
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Fig.  6-4.  Open-Loop  CGT  for  CM(02C) ,  (No  Disturbances,  Nominal  Truth  Model) 
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Fig.  6-5.  Open-Loop  CGT  for  CM(02C) ,  (Disturbances,  Nominal  Truth  Model) 
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As  seen  in  Figure  6-4,  for  this  design  model- 
conunand  model  pair,  the  open-loop  CGT  controller  achieves 
essentially  perfect  model-following  when  no  disturbances 
affect  the  dynamics  and  the  design  model  is  an  exact  repre¬ 
sentation  of  true  system  dynamics.  Figures  6-5  and  6-6 
demonstrate  the  deterioration  in  transient  performance 
when  disturbances  are  present  or  design  model  errors 
exist.  However,  the  modeled  disturbances  are  actually  suc¬ 
cessfully  overcome  by  the  open-loop  CGT  controller  in 
steady-state;  the  bounds  on  system  output  response  are 
-0.25  and  +1.15.  In  the  time-response  of  Figure  6-6  the 
actual  response  reflects  the  lower  damping  of  the  alternate 
truth  model  and  has  an  overshoot  of  about  38%. 

In  direct  contrast  to  these  open-loop  results, 
Figures  6-7,  6-8 ,  and  6-9  show  the  corresponding  cases 
employing  a  closed-loop  CGT/PI  controller.  The  same  PI 
regulator  gains  determined  in  Section  6.3.4  above  are  used 
in  the  inner-loop.  For  the  design  based  upon  CM(02C) , 
the  feedforward  gains  are  determined  as 


K  =  [56.76  157.7] 

m 

K  =9.456 
-xu 

and  Kx  =  [-15.05  -.3172] 


(6-21a) 

(6-21b) 

(6-21c) 


The  gains  for  the  design  based  upon  CM(02C) "  are  identical 

to  these  for  K  and  K  ,  while  the  feedforward  gains  on 
xu  xn 

the  command  model  states  become 


125 


Fig.  6-7.  CGT/PI  for  CM(02C) ,  (No  Disturbance,  Nominal  Truth  Model) 
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Fig.  6-8.  CGT/PI  for  CM(02C) ,  (Disturbances,  Nominal  Truth  Model) 
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Fig.  6-9.  CGT/PI  for  CM(02C) ,  (No  Disturbances,  Alternate  Truth  Model) 


pipss 


K'  «  [1072.  157.7]  (6-22) 

m 

Note  that  these  closed-loop  values  for  K'  and  K  are 

ID  ID 

related  as  before  through  equation  (6-20d) . 

The  plots  of  Figure  6-7  show  that  the  system  and 
command  model  outputs  are  coincident  (within  the  resolu¬ 
tion  of  the  printer  plots)  using  the  CGT/PI  controller. 
More  signif icantly,  the  responses  of  Figure  6-8  and  6-9 
demonstrate  considerable  improvement  in  the  transient 
response  of  the  controlled  system  in  the  face  of  non-zero 
modeled  disturbances  and  of  modeling  errors,  respectively. 
With  disturbances  the  response  is  within  the  bounds  of 
-0.1  to  +1.07  and  the  steady-state  model-following  is 
quickly  achieved,  while  even  during  the  transient  rise 
portion  the  matching  is  quite  good.  Similarly  for  design 
model  dynamics  errors,  the  transient  behavior  is  much 
improved,  an  overshoot  of  only  about  11%  occurs,  and  the 
steady-state  model -following  is  quickly  achieved. 

6.3.6  Discussion  of  Results  for  Simple  Design 
Example .  Comparing  the  various  feedforward  gain  matrices 
given  earlier  in  this  section,  it  can  be  seen  that  the 
gains  applied  to  the  disturbance  states  are  independent  of 
the  command  model.  Also,  for  given  command  model  dynamics, 
the  gains  on  the  model  states  vary  according  to  the  spe¬ 
cific  state  representation  employed,  but  the  gains  on  com¬ 
mand  model  inputs  do  not,  as  seen  in  Section  6.3.5. 
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Designs  for  CGT  controllers  based  on  command  models 
CM(02C)  and  CM(02C)'  demonstrated  the  invariance  of  con¬ 
trolled  system  behavior  with  alternative  state  representa¬ 
tions  of  specified  dynamics.  Thus,  particularly  for  SISO 
systems,  desirable  dynamic  behavior  may  be  formulated 
readily  into  command  model  state  representations,  as 
through  translation  from  transfer  function  representation 
to  a  state  representation  in  standard  controllable  form. 
Design  iteration  then  primarily  involves  the  inner-loop  PI 
regulator  or  in  some  cases  may  entail  modifications  in 
the  proposed  dynamics  but  not  in  the  representations 
achieving  such  dynamics. 

As  expected,  the  closed-loop  CGT/PI  controller 
shows  significant  advantages  over  an  open-loop  CGT  con¬ 
troller.  Transient  response  is  improved,  modeling  errors 
are  more  readily  overcome,  steady-state  exact  model¬ 
following  is  obtained  more  quickly,  and  rejection  of  modeled 
disturbances  in  the  transient  phase  of  response  is  superior. 
Additional  tests  aimed  at  exploring  the  qualities  provided 
specifically  by  the  PI  regulator  (as  opposed  to  standard 
regulator  designs)  could  be  performed.  Since  the  designs 
employed  in  this  study  assumed  use  of  PI  regulation,  such 
additional  tests  have  not  been  examined.  However,  for 
example,  it  is  clear  that  the  CGT/PI  controller  response 
when  subject  to  unmodeled  disturbances  would  be  superior 
(would  achieve  zero  steady-state  error)  to  either  the 
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open-loop  CGT  or  closed-loop  CGT  employing  simple  inner- 
loop  regulation. 

6.4  Aircraft  Flight  Control  Design 

6.4.1  Introduction.  The  designs  discussed  in  this 
section  are  based  upon  the  dynamics  corresponding  to  the 
Advanced  Fighter  Technology  Integration  aircraft,  an  F-16 
aircraft  (AFTI/F-16)  specially  modified  for  the  purpose  of 
advanced  controls  research.  Two  distinct  controllers  are 
designed  based  upon  longitudinal  dynamics  about  straight  and 
level  flight  at  Mach  0.8  and  an  altitude  of  10000.  feet. 

This  particular  flight  condition  was  chosen  as  representa¬ 
tive  of  an  operating  point  within  the  usual  air-to-air 
combat  flight  regime.  Since  the  AFTI/F-16  is  unstable  in 
pitch  in  this  regime,  open-loop  CGT  designs  are  infeasible. 
Therefore,  only  CGT/PI  designs  are  pursued. 

One  of  the  concepts  the  AFTI/F-16  aircraft  is  to 
demonstrate  is  the  feasibility  and  usefulness  of  uncon¬ 
ventional  decoupled  pitch  modes  of  motion  compared  with 
conventional  pitch  control.  Decoupled  pitch  control  is 
possible  because  of  two  independently  controlled  longitudi¬ 
nal  control  surfaces:  horizontal  tail  and  trailing-edge 
flap.  The  designs  discussed  here  are  for  a  so-called  pitch¬ 
pointing  controller  and  a  pitch  rate  controller. 

The  pitch-pointing  controller  is  in  the  class  of 
unconventional,  decoupling  controllers.  The  objective  is 
to  achieve  direct  control  of  aircraft  pitch  attitude  while 


constraining  the  flight-path  angle  to  remain  constant. 

Thus  the  aircraft  may  be  pointed  in  pitch  through  small 
angles  without  changing  the  flight  trajectory.  Such  con¬ 
trol  is  useful,  for  example,  in  air-to-air  target  tracking 
to  null  gun-aiming  error  briefly  when  the  trajectory  match¬ 
ing  is  good  but  small  errors  in  the  aim-point  persist  in 
the  pitch  plane.  The  pitch-pointing  control  mode  may  be 
practiced ly  achieved  through  appropriate  control  input  com¬ 
mand  to  the  two  independent  control  surfaces.  However, 
such  designs  are  difficult  to  achieve  using  conventional 
control  design  techniques  since  the  problem  inherently 
requires  effective  multi- input-multi-output  (MIMO)  design 
methodology.  The  CGT/PI  controller  design  technique  is 
readily  applicable  to  such  MIMO  decoupling  design  problems, 
systematically  determining  appropriate  crossfeeds  while 
assuring  closed-loop  stability. 

The  pitch  rate  controller  is  conventional  in  air¬ 
craft  flig..-  control,  except  that  for  this  case  appropriate 
utilization  of  both  control  surfaces  is  problematical  when 
using  conventional  control  design  methods.  For  the  CGT/PI 
control  design  technique,  such  conventional  control  objec¬ 
tives  represent  a  special  case  in  which  model -following  is 
desired  for  fewer  outputs  than  the  number  of  available  con¬ 
trol  inputs.  The  CGT/PI  design  described  subsequently 
abides  by  the  requirement  (imposed  in  Sections  3.2.1  and 
3.3.1)  that  the  numbers  of  outputs  and  inputs  of  the  design 
model  and  outputs  of  the  command  model  all  be  equal ,  but 


treats  only  one  output  of  each  model  as  a  non-zero  combina¬ 
tion  of  states.  Although  the  resulting  augmented  matrix 
inverse  n  of  equation  (3-47)  does  not  exist,  the  matrix 
pseudo-inverse  is  computed  instead  (Ref  32) .  The  final 
CGT/PI  controller  achieves  the  desired  model -following 
in  the  single  output  channel  specified.  For  this  case  in 
which  more  controls  are  available  than  the  number  of  out¬ 
puts  to  have  specified  dynamics,  the  pseudo- inverse  pro¬ 
vides  an  exact  controller  solution  (Ref  32) .  However,  in 
the  opposite  case  of  fewer  inputs  than  outputs,  the  use  of 
the  pseudo- inverse  would  give  an  approximate  solution  only. 

Data  describing  the  AFTI/F-16  was  obtained  from 
References  3  and  18.  Reference  18  gives  the  dimensional 
stability  derivatives  in  the  aircraft  body-axis  coordinate 
frame  in  angular  units  of  radians,  but  does  not  include 
data  for  the  trail ing-edge  flap.  The  corresponding  deriva¬ 
tives  for  the  trailing-edge  flap  were  computed  from  the 
data  of  Reference  3  which  gives  the  data  in  the  stability- 
axis  frame  and  in  angular  units  of  degrees.  The  state- 
space  representations  employed  in  describing  the  aircraft 
longitudinal  dynamics  are  given  below.  Reference  30  details 
the  derivation  of  the  linear  perturbation  model  of  aircraft 
longitudinal  dynamics  while  the  model  of  clear-air  tur¬ 
bulence  is  obtained  from  Reference  21. 

Various  models  are  used  in  specific  design  examples 
presented  below.  After  describing  the  full  truth  model 
used  for  the  longitudinal  dynamics,  several  simplifications 


are  proposed.  These  simplifications  are  used  to  demon¬ 
strate  the  effects  of  parameter  variation  and  of  unmodeled 
and  modeled  control  surface  actuator  lags  on  the  CGT/PI 
time -response.  Finally,  one  of  these  simplified  models  is 
used  as  the  design  model  for  a  Kalman  filter  design.  The 
covariance  analysis  is  with  respect  to  the  full  truth  model 


6.4.2  AFTI/F-16  Truth  Model .  The  full  truth  model 
for  the  AFTI/F-16  aircraft  employs  four  system  states,  two 
actuator  states,  and  three  turbulence  (gust)  states.  Noise 
corrupted  measurements  are  assumed  available  for  three  of 
the  system  states.  In  addition,  one  of  the  measurements 
includes  a  gust  state.  The  model  is  hereafter  referred  to 
as  AFTI(S4 ,A2,G3) ,  i.e.,  the  model  based  on  4  system  states 
2  actuator  states,  and  3  gust  states.  Computations  needed 
to  determine  the  model  parameters  were  formulated  in  sub¬ 
routines  and  included  in  the  final  object  file  of  CGTPIF. 

The  state  vector  is  defined  to  consist  of  the  fol¬ 


lowing  elements: 


(6-23) 
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in  which  the  perturbation  states  are  as  follows:  6  is  the 
aircraft  pitch  angle,  a  is  the  aircraft  angle -of -at tack, 
q  is  the  aircraft  pitch  rate,  u  is  the  forward  velocity, 

6ht  is  the  deflection  angle  of  the  horizontal  tail,  STEF 
is  the  deflection  angle  of  the  trailing-edge  flap,  a'  and 
are  angle-of-attack  gust  states,  and  qg  is  a  pitch  rate 
gust  state.  Since  the  aircraft  state  description  is 
obtained  in  two  steps,  it  is  convenient  to  first  describe 
the  actuator  and  gust  models. 

Both  control  surfaces  are  modeled  as  driven  through 
first-order  lags  representing  the  actuator  dynamics.  .  The 
transfer  function  for  the  actuator  lag  is  of  the  form 


(6-24) 


in  which  6  and  6  are  the  surface  deflection  and  com- 

S  S  j 

mand  to  the  surface  input  (actuator) ,  respectively.  The 

value  of  T  is  20.,  implying  an  actuator  time -constant  of 
s 

0.05  seconds  for  both  actuators.  The  corresponding  state 
model  is 


6 

s 


=  -  T  6  +  T  6 

s  s  s  s. 


(6-25) 


The  state  equations  for  the  gust  model  are  as 


given  in  Reference  21: 


6T  =  a"  a ' 
"9  g^g 

*  =» 

(6-26a) 

ct_  =  a.'  a' 
-9  92~9 

+  +  anPa 

(6-26b) 

q  =  a  d 

%  gd~g 

+  V® 

(6 -2  7c) 

in  which  there  is  only  one  noise  source  (na) ,  modeled  as  a 
zero-mean  white  Gaussian  noise  process  with  unit  variance. 
The  various  coefficients  are  computed  as 


(6-28a) 


(6-28b) 


a 


g 


(6-28c) 


(6-28d) 


(6-28e) 


( 6—28  f ) 


and  VT  is  the  aircraft  total  velocity  in  feet  per  second 
and  b  is  aircraft  wing  span  in  feet.  Both  and  Lw  are 
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determined  from  specifications  given  in  MIL-SPEC  8785B 
(Reference  12) :  is  the  root-mean -square  intensity  for 

the  clear  air  turbulence  (as  a  function  of  "Level"  and 
altitude)  in  feet  per  second  and  L^  is  the  scale  length 
(as  a  function  of  altitude)  in  feet. 

It  is  convenient  to  define  the  system  dynamics 
natrix  for  the  truth  model  in  two  steps.  First,  the  full 
matrix  will  be  specified  based  on  zero  values  for  stability 
derivatives  related  to  the  rate  of  change  of  angle-of- 
attack  (a)  and  neglecting  the  term  involved  dg  in  the 
pitch  rate  gust  equation.  The  resulting  matrix  is  then 
modified  to  incorporate  the  non-zero  d  stability  deriva¬ 
tives  and  the  d  term  in  the  q  state  equation.  Thus, 
g  g 

define  the  initial  dynamics  matrix  as 
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(6-29) 


A  new  matrix,  Afc,  is  next  determined  by  incorporating 
rows  1,  5,  6,  7,  and  8  intact  from  matrix  A^  and  modifying 
the  remaining  rows  as  follows: 
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At(2,I)  -  A^ (2,1) 

/ 

(l.-Z.) 

a 

(6-30a) 

A t (3,1)  =  A'(3,I) 

+ 

M.  •  A. (2,1) 
a  — t 

(6-30b) 

Afc(4,I)  =  A^ (4,1) 

+ 

X.  •  A. (2,1) 
a  — t 

(6-30c) 

At(9,I)  =  A' (9,1) 

+ 

a  •  A'(8,I) 

yd 

(6-30d) 

where  the  index  "I"  ranges  from  1  to  9  and  a  is  as  defined 
in  equation  (6-28e) .  In  both  equations  (6-29)  and  (6-30) 
the  stability  derivatives  denoted  by  "X",  "Z",  and  "M" 
are  related  to  the  forces  acting  along  the  aircraft  body 
X-  and  Z-axes  and  the  pitching  moment  (about  the  body  Y- 
axis) ,  respectively.  Subscripts  indicate  the  state  with 
respect  to  which  the  derivative  is  determined.  The  sub¬ 
scripts  "SHT"  and  "<Stefm  xefer  to  the  control  surface 
deflections  of  horizontal  tail  and  trailing-edge  flap. 

The  term  "g"  is  the  acceleration  of  gravity.  Finally, 
aQ,  Up,  and  WQ  are  the  trim  values  of  ang le-of -attack , 
forward  velocity  component  along  the  body  X-axis,  and  the 
downward  velocity  component  along  the  body  Z-axis.  Here, 
"trim"  is  the  nominal  operating  point  about  which  lineari¬ 
zation  is  accomplished  and  is  a  condition  of  wings-level, 
constant  velocity,  and  constant  altitude  flight  at  the 
Mach  number  and  altitude  specified  for  the  flight  condition 
(0.8  and  10000.  feet  respectively). 

The  values  of  all  the  parameters  of  equations 

(6-25),  (6-28),  (6-29),  and  (6-30)  are  listed  in  Table  6-1. 

The  value  of  o  listed  is  for  a  "Level  1"  turbulence,  i.e., 
w 
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-0.295  ( (1/SBC) /RAD)  Mr  =-5.88  ( (1/SEJC  ) /RAD)  X.  =1.46  ( (FT/SBC  ) /RAD) 


in  practice,  both  otf  and  Lw  are  computed  by  a  subroutine 
and  their  values  are  independent  of  the  specific  aircraft 
modeled. 

The  form  of  the  state  equation  is 

kt  =  At£t  +  BfcUt  +  Gt^t  (6-31) 

where  the  input  vector  is 


and  w  is  actually  the  scalar  n  of  equation  (6-26) ,  and 
of  strength  Qt=l . 

The  matrices  Bt  and  are  simply 


and 


St  - 


(6-34) 


Note  the  element  "q  "  of  G.  is  due  to  the  a  term  in  the 

^n  — t  g 

pitch  gust  equation  (equation  (6-27c) )  and  has  the  value 


q  »  a  •  a 

"n  9d  « 


(6-35) 


where  a  and  a  are  defined  by  equations  (6-28e)  and 
gd  n 

(6-28d) ,  respectively. 

The  measurement  equation  is  of  the  form 


— t(ti)  =  5t4t(ti)  +  ?t(ti) 


(6-36) 


where  it  is  assumed  that  states  0,  a,  and  q  are  measured 

and  the  angle-of-attack  measurement  includes  effects  of 

the  gust  state  ag,  i.e.,  measured  angle-of-attack  is 

a  ®a+a  +noise  (modeled  as  white) .  Values  for  the  measure- 
m  g 

roent  noises  for  the  AFTI/F-16  sensors  were  not  available. 


Sensor  noises  reported  in  Reference  16  were  used  instead 


142 


(as  suggested  by  the  author  of  that  paper,  the  values 
determined  from  ground  and  flight  test  for  each  sensor  were 
averaged)  to  establish  Rt,  the  covariance  of  in  equation 
(6-36)  above. 

The  values  of  the  matrices  defining  the  truth  model 
AFTI (S4 ,A2,G3)  are  (from  equations  (6-29)  through  (6-36), 
the  data  of  Table  6-1,  and  the  gust  model  employing  Level  1 
characteristics) 
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(6-37a) 


and 
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The  indicated  TDT  matrix  (see  equation  (3-21a) )  is  with 
respect  to  the  design  model  AFTI (S3,A2,G3)  described  in 
the  next  subsection. 


6.4.3  AFTI/F-16  Model  Simplifications :  Design 
Models.  Several  different  levels  of  simplification  were 
employed  in  defining  reduced  models  of  the  relevant  dynamics. 
Simplif ications  included  deletion  of  a  system  state,  elimina¬ 
tion  of  control  surface  actuator  models,  and  deletion  of 
the  gust  states.  None  of  these  simplifications  were  neces¬ 
sary — they  were  effected  in  order  to  evaluate  specific 
aspects  of  the  design.  For  the  various  reduced  models  used 
for  design  (design  models) ,  other  models  were  used  as 
"truth  models”  according  to  the  design  attribute  under 
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consideration  (e.g.,  effect  of  actuator  lags,  robustness 
of  response  when  subject  to  design  model  errors) .  The 
models  are  discussed  below  in  order  of  increasing  simplifi¬ 
cation.  Output  equations  for  each  model  are  given  in  the 
next  subsection  according  to  the  controller  to  be  designed 
(pitch-pointing  or  pitch  rate) .  All  models  in  which  the 
forward  velocity  and  gust  states  are  absent  are  considered 
to  be  in  angular  units  of  degrees  for  convenience  in  presen¬ 
tation  of  time  responses.  Models  AFTI (S4 ,A2,G3)  and 
AFTI (S3,A2,G3)  employ  angular  units  of  radians. 

6. 4. 3.1  Model  AFTI(S3,A2,G3) .  The  first  reduc¬ 
tion  of  the  truth  model  AFTI (S4 ,A2,G3)  is  achieved  by 
deleting  the  system  state  equation  for  the  perturbation  in 
forward  velocity  (u)  and  by  not  incorporating  stability 
derivatives  in  the  angle-of-attack  rate  (ct)  .  Deletion  of 
the  forward  velocity  state  from  the  design  model  is  desir¬ 
able  since  in  the  usual  case  for  flight  control  designs 
feedback  of  the  u  state  is  not  appropriate.  The  input 
vector  and  noise  process  scalar  are  as  for  AFTI (S4 ,A2,G3) . 
The  state  vector  is  simply 


This  reduction  is  equivalent  to  employing  the 
matrix  of  equation  (6-29)  with  row  4  and  column  4 
deleted  and  with  equation  (6-30d)  computed  as 


A'(8,I)  *  A^ (8,1)  +  a£ 


A'(7,I) 


(6-39) 


with  the  index  "I"  ranging  from  1  to  8.  In  addition, 
row  4  of  the  matrices  and  is  deleted  and  column  4  of 
the  matrix  Ht  is  deleted.  As  noted  in  the  previous  sub¬ 
section,  matrix  TDT  of  equation  (6-37g>  gives  the  corres¬ 
pondence  of  states  between  the  truth  model  AFTI (S4 ,A2,G3) 
and  this  design  model. 

The  parameters  defining  AFTI (S3,A2,G3)  are  computed 

in  subroutines  incorporated  into  the  object  file  of  CGTPIF. 

The  values  of  the  matrices  are  given  below.  Note  that 

since  disturbances  are  included  in  the  system  state, 

matrices  E  .  A  ,  G  ,  Q  ,  E  ,  and  H  of  the  design  model  do 

—x  — n  — n  — n  — y  — n 

not  exist.  Also,  the  matrix  TNT  of  the  truth  model 
(equation  (3-21b) )  does  not  exist,  since  the  equivalences 
for  the  disturbance  states  are  included  in  TDT  (equation 
6-37g) ) , 
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Design  model  AFTI(S3 ,A2,G3)  is  used  for  design  of 
a  pitch-pointing  CGT/PI  controller  and  the  Kalman  filter. 

For  this  purpose,  deletion  of  the  u  state  is  a  desirable 
simplification  in  practice,  since  this  ensures  that  it  will 
not  be  used  as  a  feedback  variable.  However,  elimination  of 
coefficients  in  a  is  not  actually  necessary  but  was  done  so 
that  the  design  and  truth  models  differed  both  in  dimension 
(u  state  deletion)  and  in  specific  matrix  elements. 
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6. 4. 3. 2  Model  AFTKS3.A2)  .  The  model  AFTI(S3,A2) 
is  used  only  in  CGT/PI  design  and  assumes  full  state  feed¬ 
back  is  available.  Since  a  Kalman  filter  is  not  designed 
for  this  case  the  matrices  of  the  measurement  equation 
(equation  (3-6))  are  not  needed.  Also,  the  gust  states 
and  associated  driving  noise  are  deleted  from  the  system 
state  equation  to  facilitate  consideration  strictly  of 
aircraft  and  actuator  dynamics.  Since  the  model  entails 
few  matrices,  it  is  convenient  to  enter  it  directly  from 
the  terminal.  The  state  vector  is 


(6-41) 


and  the  input  vector  is  the  same  as  for  AFTI (S4 ,A2,G3) 
and  AFTI (S3 , A2 ,G3) . 

The  values  of  the  matrices  are 
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6. 4. 3. 3  Model  AFT I  (S3)  ♦  The  simplest  model  used 
was  AFTI(S3),  which  includes  only  the  three  system  states 
of  the  simplified  models  discussed  above.  This  reduction 
provides  the  capability  to  evaluate  the  design  for  the 
simplest  aircraft  representation  and  judge  the  effects  of 
actuator  lags.  The  state  vector  is 


0 


x  = 


a 


(6-43a) 


and  the  input  vector  is 


u  = 


HT 


TEF 


(6-43b) 


The  matrices  are  conveniently  entered  directly  at 
the  computer  terminal  and  are 
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and 


B  = 


0.  0. 

-0.179  -0.295 

-25.3  -5.88 


(6-44b) 


Design  model  AFTI(S3)  is  used  for  the  CGT/PI  con¬ 
troller  designs.  Effects  of  actuator  errors  in  modeling 
(or  parameter  variations  in  flight)  are  evaluated  by 
defining  three  different  "truth  models"  of  the  same  order 
as  AFTI(S3)  but  with  different  matrix  coefficients.  The 
design  based  upon  AFTI(S3)  is  not  modified  in  performing 
the  evaluations  with  respect  to  the  various  "truth  models." 
Two  of  the  models  were  obtained  by  modifying  several  of 
the  stability  derivatives  at  the  given  flight  condition, 
while  t.ie  third  represents  a  markedly  different  flight  con¬ 
dition.  These  will  now  be  delineated  in  detail. 


6. 4. 3. 3.1  Truth  Model  TM(S3) +.  Truth  model 

TM(S3)  +  was  obtained  from  model  AFTI(S3)  by  increasing  the 

values  of  the  derivatives  Z  .  M  .  and  M  by  about  20%. 

a  a '  <3 

The  resulting  dynamics  matrix  is 
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The  control  matrix  Bt  is  identical  to  B  of  equation 
(6-43b) . 
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6. 4. 3. 3. 2  Truth  Model  TM(S3)  .  Truth  model 

TM(S3) ”  was  obtained  from  model  AFTI(S3)  by  decreasing  the 

values  of  the  derivatives  Z  ,  M  ,  and  M  by  about  20%. 

ot  ot  g 

The  resulting  dynamics  matrix  is 
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(6-46) 


The  control  matrix  is  identical  to  B  of  equation 
(6-43b) . 


6 . 4 . 3 . 3 . 3  Truth  Model  TM(S3) " .  Truth  model 
TM(S3)  "*  was  obtained  using  the  same  simplifications  as  in 
AFTI (S3)  but  employing  data  for  a  flight  condition  of 
Mach  0.6  at  an  altitude  of  30000.  feet.  Substantial 
changes  in  both  the  stability  and  control  derivatives 
result  from  such  an  extreme  change  in  nominal  flight  con¬ 
dition.  Useful  implementations  would  employ  gain  schedul¬ 
ing  in  order  to  accommodate  such  large  changes  in  aircraft 
operating  point.  The  dynamics  and  control  matrices  become 
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0.997 

0. 

0.500 

-0.274 

0. 

0. 

St  - 

-6.06E-2 

-0.112 

-5.82 

-0.219 

(6-4 7a) 


(6-4 7b) 
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6.4.4  Model  Output  Equations 
6. 4. 4.1  Output  Equations  for  Pitch-Pointing 
Control .  Pitch-pointing  control  entails  output  response 
in  pitch  attitude  (6)  while  constraining  flight-path  angle 
at  a  constant  value.  Referring  to  the  flight-path  angle 
as  y  and  employing  the  usual  approximation  derived  for 
wings-level  flight  (Ref  30)  that 

Y  ~  6  -  a  (6-48) 


yields  the  following  output  matrices  (see  equation  (3-5)): 
For  AFTI (S3 ,A2 ,G3) , 


1.000  0.  0. 

• 

O 

• 

o 

c  - 

1 

1.000  -1.000  0. 

o 

• 

o 

• 

For  AFTI (S3  # A2) , 

1.000  0.  0. 

o 

. 

o 

C  = 

;  ! 

1.000  -1.000  0. 

•  0.  O.j 

and  for  AFTI (S3) , 

1.000  0.  0. 

C  = 

1.000  -1.000  0. 

The  matrices 

D  and  E  of  equation  (3-5)  c 
— y  — y 

matrices  for 

all  models. 

0.  0. 

0.  0. 


0. 

0. 

(6-4*93) 


(6-49b) 


(6-49c) 
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Pitch  rate  control  simply  entails  output  response  in  pitch 
rate  with  no  constraints  on  other  potential  output  vari¬ 
ables.  Thus  the  output  matrix  actually  specifies  only  one 
output  (q) ,  while  the  second  output  row  is  included  only 
to  maintain  dimensional  compatibility  with  the  control 
input  dimension: 


For  AFTI (S3,A2,G3) , 


C  = 


0.  0, 

0.  0, 


For  AFTI (S3, A2), 


C  * 


and  for  AFTI (S3) , 


C  = 


1.000 

0. 


0.  0.  1.000 

0.  0.  0. 


0.  0.  1.000 

0.  0 .  0 . 


0.  0. 

0.  0. 


0. 

0. 


0. 

0. 


0.  0. 

0.  0. 


0. 

0. 

(6-50a) 


(6-50b) 


(6-50c) 


The  matrices  and  of  equation  (3-5)  are  both  zero 
matrices  for  all  models. 

6.4.5  Command  Models  (Equations  (3-24)  and  (3-25)) 
6. 4. 5.1  Command  Model  for  Pitch-Pointing  Control. 


Pitch-pointing  control  is  readily  specified  as  requiring 
decoupled  output  response  with  respect  to  each  of  the 
individual  command  input  channels.  Such  response 


characteristics  are  conveniently  formulated  entirely  in 
terms  of  diagonal  command  model  matrices  to  achieve  first' 
order  response  trajectories.  The  general  form  in  this 
case  (of  two  inputs  and  outputs)  is. 


0 .  a 


(6-51a) 


B_  = 


b1  0. 

0.  b_ 


(6-51b) 


— m  = 


0.  c 


(6-5lc) 


The  matrix  D  is  the  zero  matrix. 

For  this  command  model  the  transfer  functions  for 

outputs  Ym(l)  and  ym(2)  are: 

For  input  u  (1)  , 
m 


y_(l)  c. 
*  m  _  _ 1 

U_(l)  ”  8 

m 


s  -  a. 


(6-52a) 


yJ2' 

*  m 

srnr 

m 


(6-52b) 


and  for  input  u  (2), 
m 


y_u> 

m 


(6-53a) 


(6-53b) 


ym(2)  ,  c2 :  b2 

HTT2T  8  -  a. 

The  values  of  a^ ,  a2,  b2,  c^,  and  c2  may  all 

be  chosen  according  to  the  specific  responses  desired. 

Of  course,  a 1  and  a2  must  both  be  chosen  as  negative  for 
stable  output  response.  The  products  (c^-b^)  and  (c2*b2) 
will  affect  the  values  of  the  CGT  feedforward  gains  on  the 
command  model  inputs  um(l)  and  um(2),  respectively  (equa¬ 
tions  (6-52a)  and  6-53b) ) .  The  specific  choices  of  c. , 
b1 ,  and  c2,  values  impact  the  CGT  gains  and  the  system 
response  only  through  their  respective  products,  i.e., 
scaling  similarity  transformations  do  not  affect  results; 
thus  their  individual  values  are  arbitrary  and  selection 
according  to  convenience  is  appropriate.  Finally, 
decoupling  is  generally  not  perfectly  achievable.  But  in 
any  CGT  design  in  which  the  PI  gains  are  fixed  in  value, 
the  ratio  of  the  steady- state  magnitude  of  the  commanded 
output  to  the  maximum  transient  magnitude  achieved  in  the 
nominally  zero  output  channel  is  constant  for  all  values 
of  c1,  b^,  c2,  and  b2  since  the  feedforward  gains  for  each 
command  input  channel  are  independent  for  command  models 
defining  decoupled  output  response. 

From  the  definition  of  the  output  matrix  of  equa¬ 
tion  (6-49),  output  1  is  the  desired  pitch-pointing  output. 
The  desired  form  of  response  is  thus  defined  by  equations 
(6-52a,b) ;  command  input  um(l)  then  commands  pitch-pointing. 


The  values  of  a^,  b^,  and  must  be  selected  to  give 
desirable  response  while  the  values  of  a2,  b 2,  and  c2  are 
arbitrary . 

The  value  of  a^  was  selected  as  -5.,  giving  a 
first-order  pitch-pointing  command  model  response  with  a 
rise-time  of  about  0.6  seconds  for  a  step  command  input. 
When  used  in  conjunction  with  design  models  AFTI (S3)  and 
AFTI (S3, A2)  the  numerator  of  equation  (6-52a)  was  selected 
to  give  an  output  to  input  ratio  of  unity  at  low  frequency 
(Ci’bi*5!  ajJ )  .  However,  since  the  angular  units  for  design 
model  AFTI(S3,A2,G3)  are  radians  to  correspond  to  truth 
model  AFTI (S4 ,A2,G3) ,  by  design  choice  a  command  model 
output  to  input  ratio  of  0.02  was  selected  for  designs  with 
AFTI  (S3,A2,G3)  (ci  *^*0 . 02 •  |  a^^  | )  .  The  values  of  a2,  b2, 
and  c2  were  arbitrarily  set  equal  to  the  values  of  a^,  b^, 
and  c^ ,  respectively. 

For  all  pitch-pointing  designs  the  command  model's 
dynamics  matrix  is 


(6-54) 


In  conjunction  with  design  models  AFTI (S3)  and  AFTI(S3,A2) 
the  matrices  B  and  C  are 


9 


and 


—  — 

1.  0. 

i  c  ■ 

I  “* 

0.  1. 

1 

—  — 

(6-55b) 


In  conjunction  with  design  model  AFTI (S3 , A2,G3) ,  the 


matrices  B  and  C  are 
— m  — m 


2*  = 


0.1  0. 

0.  0.1 


(6-56a) 


and 


— m  ~ 


1.  0. 

0.  1, 


(6-56b) 


Matrix  D  is  the  zero  matrix. 
— m 


6 . 4 . 5 . 2  Command  Model  for  Pitch  Rate  Control .  The 
same  command  model  form  described  above  was  also  used  for 
the  pitch  rate  controller  designs.  In  this  case  the  pro¬ 
duct  of  c2  and  b2  was  selected  as  zero  since  only  the  first 
of  the  two  outputs  is  actually  specified  in  the  design 
model  output  equation.  The  pitch  rate  controller  was 
designed  using  only  the  AFTI (S3)  and  AFTI (S3, A2)  design 
models  and  an  output  to  input  response  ratio  of  unity  was 
chosen.  The  matrices  defining  the  pitch  rate  command  model 
are 


*m 


-5.  0. 

0.  -5, 


(6-57a) 
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2* 


5.  0, 

0.  5, 


and 


1. 

0. 


0. 

0, 


Matrix  D_  is  the  zero  matrix. 
— m 


(6-57b) 


(6-57c) 


6.4.6  CGT/PI  Pitch-Pointing  Controller  Design . 
CGT/PI  controller  designs  achieving  pitch -pointing  con¬ 
trol  are  descr.-* ...  M  briefly  below.  Designs  are  determined 
for  design  models  AFTI (S3) ,  AFTI (S3, A2),  and  AFTI (S3 ,A2,G3) . 
An  important  observation  is  that  designs  which  do  not 
account  for  the  lags  introduced  by  the  control  surface 
actuators  are  unstable  when  required  to  control  the  air¬ 
craft  through  actuators  with  the  lags  of  the  AFTI/F-16 
aircraft.  Designs  which  did  include  actuator  models  were 
stable  but  did  show  response  behavior  somewhat  inferior  to 
that  achievable  without  actuators.  When  the  actuators  are 
modeled  as  having  lag  time-constants  several  times  smaller 
than  those  of  AFTI/F-16  (x=0.05  seconds),  designs  which 
neglected  actuator  dynamics  behaved  very  well  when  con¬ 
trolling  the  aircraft  through  the  actuators. 


6. 4. 6.1  AFTI (S3)  Pitch-Pointing  CGT/PI.  Consistent 
with  usual  experience  in  aircraft  flight  control  designs 
for  pitch  motion,  the  utility  of  pitch  rate  damping  was 
quickly  determined  in  the  PI  and  CGT/PI  designs  for  the 
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pitch-pointing  mode.  Figure  6-10  illustrates  the  response 
of  an  initial  PI  regulator  design  with  quadratic  weights 
of  200.  on  .both  outputs  and  1.  on  input  magnitudes  and 
rates.  The  initial  conditions  were  1.  on  output  1  and  2. 
on  output  2.  Plot  symbols  1  and  2  are  outputs  1  and  2, 
respectively.  The  response  in  pitch  attitude  (output  1) 
is  clearly  poorly  damped.  The  CGT/PI  controller  response 
with  the  same  PI  regulator  is  shown  in  Figure  6-11.  Plot 
symbols  1  and  2  are  the  aircraft's  actual  and  desired  pitch 
attitude  response  and  symbols  3  and  4  are  the  aircraft's 
actual  and  desired  flight-path  angle  deviations.  Compared 
to  the  desired  response,  the  actual  pitch-pointing  response 
initially  lags,  then  overshoots,  and  finally  oscillates 
about  the  final  value.  The  alternate  channel  response  is 
small  but  oscillatory.  While  the  response  is  good  overall, 
clearly  improved  pitch  damping  is  needed. 

Figures  6-12  and  6-13  show  the  corresponding 
responses  for  a  new  PI  regulator  design.  This  final  PI 
design  employs  identical  weightings  to  that  discussed 
above,  but  also  includes  a  weight  of  50.  on  the  pitch  rate 
state.  Note  the  dramatic  improvement  in  the  PI  regulator's 
response  to  initial  conditions  both  in  speed  of  response 
and  damping  of  oscillations.  The  pitch-pointing  response 
of  the  CGT/PI  controller  shown  in  Figure  6-13  is  nearly  an 
exact  match  with  the  desired  response.  Also,  the  alternate 
output  channel  response  is  very  small  (its  maximum  magni¬ 
tude  is  about  0.025  as  in  the  previous  case)  and  is 
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Fig.  6-10.  AFTI(S3)  PI  Regulator — No  Weight  on  Pitch  Rate 
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Fig.  6-11.  AFTI (S3)  Pitch-Pointing  CGT/PI — No  Weight  on  Pitch  Rate 
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Fig.  6-12.  AFTI (S3)  PI  Regulator — Weight  on  Pitch  Rate 
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rapidly  damped  (about  0.14  seconds  compared  to  more  than 

1.0  seconds  in  the  previous  case) .  Figure  6-14  shows  the 

responses  of  the  control  surfaces  for  the  CGT/PI  response 

of  Figure  6-13.  Plot  symbols  1  and  2  are  the  horizontal 

tail  and  trailing-edge  flaps,  respectively.  For  a  steady- 

state  pitch-pointing  command  of  l.°  steady-state  surface 

deflections  of  about  1.7°  for  the  tail  and  -6.9°  for  the 

flap  are  required.  Note  that  the  tail  surface  has  an 

initially  large  deflection  at  the  first  controller  sample 

time.  As  mentioned  earlier  in  Section  6.1.2,  this  initial 

commanded  deflection  of  large  magnitude  is  typically  the 

most  important  design  constraint  for  the  CGT/PI  controller 

find  is  due  to  the  feedforward  gain  K  on  the  command  model 

~*u 

input.  Reduction  of  this  initial  input  command  may  require 
that  PI  gains  be  reduced  and/or  that  command  model  dynamics 
be  changed. 

The  set  of  quadratic  weights  of  this  PI  regulator 
was  used  for  all  the  pitch-pointing  designs.  The  PI  gains 
for  design  model  AFTI(S3)  are 


and 


K  = 
—x 


-41.72  19.47  -1.361 

76.73  -83.11  0.4885 


-0.1603  -0.8697 

-2.963  3.069 


(6-58a) 


(6-58b) 
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Fig.  6-14.  AFTI (S3 )  CGT/PI  Control  Surface  Deflections 


and  the  feedforward  CGT/PI  gains  are 


-14.68  -16.85 

-14.67  70.88 


(6-59a) 


and 


-5.736  -4.451 

1.419  19.11 


(6-59b) 


Note  that  for  a  unit  step  input  on  command  model  input  1 

(pitch-pointing  command  input)  element  (1,1)  of  matrix 

Kx  is  the  value  of  the  command  input  during  the  first 

sample  period  following  the  command.  Thus  in  order  to 

modify  this  initial  magnitude,  designs  should  be  sought 

which  give  different  values  for  K  (1,1). 

xu 

The  next  three  figures  demonstrate  the  response  of 
this  same  CGT/PI  controller  when  the  design  model  includes 
parameter  errors.  Figure  6-15,  6-16,  and  6-17  are  the 
CGT/PI  responses  evaluated  using  truth  models  TM(S3)+, 
TM(S3)~,  and  TM(S3)',  respectively.  As  before,  plot 
symbols  1  and  2  are  actual  and  desired  pitch-pointing 
response  and  symbols  3  and  4  are  actual  and  desired  alter¬ 
nate  channel  response.  The  first  two  figures  show  that 
for  significant  errors  (±20%)  in  key  model  parameters  the 
controlled  response  is  virtually  identical  to  that  for 
nominal  parameter  values.  Figure  6-17  demonstrates  that 
even  for  very  great  differences  in  parameter  values  due 
to  a  large  change  in  flight  condition  the  CGT/PI  response. 
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Fig.  6-15.  AFTI (S3)  Pitch-Pointing  CGT/PI  WRT  TM(S3) 
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Fig.  6-16.  AFTI (S3)  Pitch-Pointing  CGT/PI  WRT  TM(S3) 
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Fig.  6-17.  AFTI (S3 )  Pitch-Pointing  CGT/PI  WRT  TM(S3) 


■ 


1 


while  seriously  degraded,  is  not  unusable.  Clearly, 
gain  scheduling  is  appropriate  (and  is  typically  employed) 
in  order  to  accommodate  such  extremes  in  operating  condi¬ 
tions  . 


6. 4. 6. 2  AFT I (S3 ,A2)  Pitch-Pointing  CGT/PI . 

Figure  6-18  shows  the  response  of  the  PI  regulator  deter¬ 
mined  using  the  quadratic  weightings  described  in  the  sub¬ 
section  above.  Model  AFTI(S3,A2)  is  used  both  for  design 
and  evaluation  in  this  case.  The  response  shown  is  the 
analog  of  the  response  plotted  in  Figure  6-12.  Note  that 
inclusion  of  actuator  dynamics  leads  to  slightly  slower 
and  less  well-damped  pitch  attitude  response. 

The  corresponding  CGT/PI  response  is  shown  in 
Figure  6-19,  which  is  the  analog  of  the  response  of  Figure 
6-13.  The  pitch  attitude  response  is  degraded  in  compari¬ 
son  to  what  was  achieved  when  the  inputs  were  assumed  to 
be  the  control  surface  deflections  directly.  The  maximum 
magnitude  of  the  deviation  in  flight  path  angle  remains 
small  (about  0.05). 

Although  the  controller's  inputs  to  the  actuators 
were  nearly  identical  to  those  generated  to  drive  the  con¬ 
trol  surfaces  in  the  AFTI(S3)  design,  the  surface  deflec¬ 
tions  actually  attained  in  the  initial  phase  of  response 
were  quite  different.  Figure  6-20  shows  the  deflections 
of  the  horizontal  tail  and  trailing-edge  flaps  (plot  symbols 
1  and  2,  respectively) .  These  deflections  may  be  compared 
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Fig.  6-18.  AFTI(S3,A2)  Pitch-Pointing  PI  Regulator 
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Fig.  6-20.  AFTI(S3,A2)  Pitch-Pointing  CGT/PI  Control  Surface  Deflections 


to  Figure  6-14,  for  which  actuators  were  not  modeled.  The 
lags  introduced  by  the  control  actuators  account  for  the 
differences  in  the  control  surface  deflections  and  in  the 
CGT/P1  controller's  response  during  the  initial  transient 
phase  of  response. 

The  PI  gains  for  the  AFTI(S3,A2)  pitch-pointing 
design  model  are 

1.282  5.712E-2 

5.502E-2  1.045 

(6-60a) 

(6-60b) 


-38.87  19.69  -1.687 

K  = 

x  78.07  -83.61  0.4227 

and 

7.00E-3  -0.8716 

K  = 

2  -2.929  3.068 


and  the  CGT/PI  feedforward  gains  are 


-9.487 

-18.44 


-14.92 

61.13 


and 


-5.908 

-1.059 


-8.554 

36.43 


(6-61a) 


(6 -6  lb) 


6. 4. 6. 3  AFT I (S3 , A2 ,G3)  Pitch-Pointing  CGT/PI. 
The  response  of  the  CGT/PI  controller  designed  based 
upon  model  AFTI (S3 ,A2 ,G3) ,  using  the  same  quadratic 
weights  as  previously,  and  evaluated  with  respect  to  the 
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truth  model  AFTI (S3,A2,G3)  is  shovm  in  Figure  6-21.  It  is 
the  analog  of  the  previous  plots  of  Figure  6-13  and  6-19 
for  design  models  AFTI (S3)  and  AFTI(S3,A2),  respectively. 
The  response  shovm  is  in  units  of  radians.  Clearly,  there 
is  no  significant  degradation  in  response  compared  to  that 
shovm  in  Figure  6-19  for  AFTI (S3, A2).  Thus  deletion  of  the 
forward  velocity  perturbation  state  and  neglect  of  effects 
due  to  d  derivatives  is  of  less  impact  than  failure  to 
model  actuator  dynamics. 

The  gains  for  the  PI  regulator  are 


-38.87 

19.69 

-1.687  |  1.282 

1 

5.710E-2 

78.03 

-83.61 

0.4227  ! 

> . 500E-2 

1.045 

0.1851 

-62.16 

0.0294 

(6-6 2a) 

-0.8324 

280.0 

7.690E-3 

and 


= 


7.040E-3  -0.8717 

-2.929  3.068 


(6-6 2b) 


and  the  CGT/PI  feedforward  gains  are 


K 

—x 


m 


and 


K 

— x. 


u 


-9.482  -14.92 

-18.47  61.13 

-0.1181  -0.1711 

-0.0215  0.7287 


(6-6 3a) 


(6-63b) 
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The  gains  on  system  and  actuator  states  in  the 

matrix  K  and  the  gains  of  K  and  K  are  essentially 
x  z  X 

in 

identical  to  those  of  equations  (6-60a,b)  and  (6-61a) 

determined  for  model  AFTI(S3,A2).  Also,  the  input  gains 

K  of  equation  (6-6 3b)  above  are  identical  to  those  of 
xu 

K  in  equation  (6 -6 lb)  scaled  by  a  factor  of  0.02.  Note 
u 

that  this  is  the  same  factor  as  between  output  to  input 

ratios  for  the  respective  command  models  (see  Section 

6. 4. 5.1  above).  Finally,  the  large  gains  on  the  angle-of- 

attack  gust  state  (a  )  in  the  K  matrix  must  oe  noted. 

g  -x 

These  imply  that  the  regulation  of  the  state  will  be 
very  tight  and  accurate.  Timely  estimation  of  that  state 
in  the  Kalman  filter  will  be  crucial  to  achieving  desirable 
response  to  turbulence. 

6.4.7  CGT/PI  Pitch  Rate  Controller  Design. 

CGT/PI  pitch  rate  controllers  determined  for  design  models 
AFTI(S3)  and  AFTI(S3,A2)  are  described  briefly  below.  The 
effects  due  to  actuator  lags  for  this  controller  type  are 
as  noted  previously  (Section  6.4.6)  for  the  pitch-pointing 
controller.  Although  the  n  matrix  of  equation  (3-47) 
proved  to  be  rank  deficient  (non- invertible)  because  of  the 
degenerate  output  matrix  C  defined  in  equation  (6-50) , 
the  matrix  pseudo-inverse  was  computed  (automatically  in 
CGTPIF)  and  the  designs  were  successfully  accomplished 
(Ref  32) .  An  alternative  would  be  to  specify  desired 
dynamics  for  a  second  output.  However,  for  the  pitch  rate 
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controller,  it  is  appropriate  to  specify  only  a  desired 
pitch  rate,  allowing  the  natural  aircraft  dynamic  relation¬ 
ships  to  determine  the  responses  in  the  remaining  outputs. 


6. 4. 7.1  AFTI  (S3)  Pitch  Rate  CGT/PI.  Quadratic 
weights  providing  desirable  response  qualities  in  the  PI 
regulator  and  CGT/PI  controller  were  readily  determined. 
Weights  of  150.  on  the  output  (output  1  only)  and  50.  on 
states  6  and  a,  and  of  1.  on  all  control  input  magnitudes 
and  rates  (equal  weighting  for  both  inputs  was  arbitrarily 
selected  and  proved  adequate)  yielded  PI  feedback  gains  of 


-0.0357  -0.4087  -1.925 

-0.00909  0.4542  -0.4875 


and 


K 


— z 


-0.8266 

-0.2055 


(6-64a) 


(6-64b) 


Figures  6-22  and  6-23  show  the  CGT/PI  output 
response  and  commanded  control  deflections,  respectively. 
In  Figure  6-22  plot  symbols  1  and  2  are  the  actual  and 
desired  pitch  rate  responses;  in  Figure  6-23  plot  symbols 
1  and  2  are  the  horizontal  tail  and  trailing-edge  flap 
deflections.  The  pitch  rate  response  achieves  essentially 
perfect  matching  to  the  command  model  output  response. 

The  commanded  pitch  rate  of  1°  per  second  entailed  maximum 
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Fig.  6-23.  AFTI(S3)  Pitch  Rate  CGT/PI  Control  Surface  Deflections 


deflections  of  about  -0.37°  and  +0.5°  for  the  tail  and  flap, 
respectively. 

The  CGT/PI  feedforward  gains  are 


(6-6  5a) 


and 


(6 -6  5b) 


Figures  6-24,  6-25,  and  6-26  illustrate  the  per¬ 
formance  of  the  CGT/PI  pitch  rate  controller  when  subject 
to  parameter  errors  in  the  design  model.  They  show  the 
response  of  the  CGT/PI  design  described  above  when  evalu¬ 
ated  with  respect  to  truth  models  TM(S3)+,  TM(S3)~  and 
TM(S3) *  .  Response  of  the  controlled  system  when  subject 
to  substantial  errors  (±20%)  in  the  values  of  several  key 
parameters  yields  response  virtually  indistinguishable 
from  that  for  nominal  design  values.  In  the  case  of  the 
response  of  Figure  6-26,  even  with  gross  errors  in  most  of 
the  design  model  parameters  (due  to  a  large  change  in 
flight  conditions)  the  response  is  still  usable  despite 
some  degradation. 

6. 4. 7. 2  AFT I (S3 ,A2)  Pitch  Rate  CGT/PI.  Using  the 
same  quadratic  weights  as  discussed  immediately  above,  the 
following  PI  regulator  gains  were  computed: 
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Fig.  6-24.  AFTI (S3)  Pitch  Rate  CGT/PI  WRT  TM(S3) 
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Fig.  6-25.  AFTI(S3)  Pitch  Rate  CGT/PI  WRT  TM(S3) 
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-26.  AFTI (S3)  Pitch  Rate  CGT/PI  WRT  TM(S3) 
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-9.440E-3  0.9006  -0.5936j 0.1894  0.6637 

(6-66a) 


-0.6578  0. 

-0.1803  0. 


(6-66b) 


The  corresponding  CGT/PI  feedforward  gains  are 


K 

— x_ 
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-1.698  0. 

-0.8371  0. 


and 


-1.151  0. 

0.8768  0. 


(6-6 7a) 


(6-6 7b) 


The  response  achieved  by  this  CGT/PI  controller 
is  shown  in  Figure  6-27.  Note  that  a  short-lived  and  small- 
magnitude  reversal  in  output  response  occurs  at  about 
mid-value  in  the  rise  portion  of  the  transient  response. 
Additional  investigation  of  the  state  and  control  surface 
responses  would  be  necessary  in  order  to  determine  the 
cause  of  this  observed  reversal  in  pitch  rate  response. 

The  PI  regulator  was  redesigned  to  overcome  the 
observed  reversal  in  output  response.  Increasing  the  out¬ 
put  quadratic  weight  from  100.  to  200.  and  the  weights 
on  6  and  a  from  50.  to  100.  gave  much  improved  response. 

The  alternate  PI  gains  are 
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-27.  AFTI (S3, A2)  Pitch  Rate  CGT/PI 
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The  corresponding  CGT/PI  feedforward  gains  are 
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and 
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(6-6 9a) 


(6-69b) 


The  pitch  rate  response  achieved  by  this  alternate 

CGT/PI  design  is  shown  in  Figure  6-28.  Although  the 

response  is  very  good,  this  particular  design  would  be 

unsatisfactory  due  to  excessive  control  deflection.  The 

flap  is  deflected  a  naximum  of  about  1.3°  during  an  early 

transient  for  a  pitch  rate  command  of  one  degree  per 

second.  Element  (2,1)  of  matrix  K  determines  the  initial 

-xu 

command  input  to  the  flap  actuator,  so  design  iterations 

should  be  directed  at  achieving  a  smaller  value  for  K  (2,1) 

u 

without  increasing  K  (1,1)  to  the  point  that  would  cause 

xu 

difficulty  with  tail  deflections  (deflection  limits  for 
the  horizontal  tail  are  ±25.°  while  those  for  the 


189 


-7 


-j 


_ 

m 

tit 

0000 *T  0008* 

0009* 

000** 

0002* 

0000*0 

♦2 

9 

9 

9 

9 

9 

+2  *  *  * 

*9* 

t  t  t 

*9* 

t 

i  t 

*9* 

t  t  «  i+t 

1 

t  t 

*9 

92  . 

9  . 

9 

9 

9 

9 

•  +z 

•9 

9 

9 

9 

9 

92 

9 

9 

9 

9 

9 

92 

9 

9 

9 

9 

9 

92 

9 

9 

9 

9 

9 

92 

9 

9 

9 

9 

9 

92 

9 

9 

9 

9 

9 

92 

9 

9 

9 

9 

9 

92 

9 

9 

9 

9 

9 

912*  *  * 

*9* 

t  t  : 

*9* 

• 

t  : 

*9* 

t  t  t  t+: 

t 

:  : 

*9 

9  2 

9 

9 

9 

9 

9 

9  2 

9 

9 

9 

9 

9 

9  2 

9. 

9 

9 

9 

+ 

9  2 

9 

9 

9 

9 

9 

9  2 

9 

9 

9 

9 

9 

9  2 

9 

9 

9 

9 

9 

9  T2 

9 

9 

9 

9 

9 

9  z 

9 

9 

9 

9 

9 

9  2 

9 

9 

9 

9 

9 

9*  £  *  * 

*9* 

t  :  : 

*9* 

t 

:  t 

*9* 

:  :  t 

5 

:  t 

*9 

9  2 

9 

9 

9 

9 

9 

9  2 

+ 

9 

9 

+ 

9 

9  2 

9 

9 

9 

9 

9 

9  2 

9 

9 

+ 

9 

9 

9  2 

9 

9 

9 

9 

9 

9  12 

9 

9 

9 

9 

9 

9  2 

9 

9 

9 

9 

9 

9  12 

9 

9 

9 

9 

9 

9  2 

9 

9 

9 

9 

+ 

+ :  :  :  : 

2*9* 

t  :  : 

*9* 

: 

s  : 

*9* 

t  :  :  :+: 

t 

s  : 

*9 

'  + 

129 

9 

9 

9 

9 

9 

2 

9 

9 

9 

9 

9 

92 

9 

9 

9 

9 

9 

9 

2 

9 

9 

9 

9 

9 

9 

2T 

9 

9 

9 

9 

9 

9 

21 

9 

9 

9 

9 

9 

9 

2 

9 

9 

9 

9 

9 

9 

2 

9 

+ 

9 

9 

9 

9 

T2+ 

9 

9 

9 

+  :  :  j  : 

*9* 

•  •  • 

«  •  • 

192 

• 

• 

•  • 

*9* 

:  :  :  :+: 

• 

• 

•  • 

•  • 

*9 

;  9 

9 

91 

2 

9 

9 

9 

9 

9 

9 

T 

2 

9 

9 

9 

9 

9 

9 

12  9 

9 

9 

,  ♦ 

’  9 

9 

92 

9 

9 

9 

9 

9 

9 

2  T  9 

9 

9 

9 

9 

9 

2  9 

I 

9 

♦  9 

9 

9 

9 

92 

T 

9 

9 

9 

9 

9 

9 

2 

T9 

9 

9 

9 

9 

9 

2 

3TV08 

00*1 
86* 

96* 

96* 

26  * 

06* 

88* 

98* 

98* 

ZB' 

08* 

BL  * 

9  2* 

9Z.* 

22* 

0L * 

89* 

99* 

99* 

29* 

09* 

8S  * 

9S  * 

95* 

2S  * 

OS* 

89* 

99* 

99* 

29* 

09* 

BE¬ 
OS  * 

9E* 

2E  * 

OE* 

82* 

92* 

92* 

22* 

02* 

81* 

9T  * 

91* 

21* 

01* 

80* 

90* 

90* 

20* 

00*0 


1 


I i 


Fig.  6-28.  AFTI(S3,A2)  Pitch  Rate  CGT/PI— Alternate  PI  Regulator 


trailing-edge  flap  are  +20.°  and  -23.°).  Different  values 
of  quadratic  weights  for  q,  a,  and  6  in  the  range  of  those 
used  in  the  PI  designs  discussed  above  would  undoubtedly 
achieve  a  satisfactory  design. 


6.4.8  Kalman  Filter  Design.  A  Kalman  filter 
design  for  the  AFTI/F-16  longitudinal  flight  controllers 
was  readily  achieved  using  the  design  model  AFTI (S3 ,A2,G3) 
described  above  in  Section  6. 4. 3.1  (based  on  measurements 
of  0,  a,  and  q  system  states) .  The  computed  filter  gain 
matrix  is  (see  Equation  3-110) , 
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3.477E-2 

2.068 

-0.1628 

(6-70) 


Note  the  two  rows  of  zero  elements  are  the  filter 
measurement  gains  in  updating  the  actuator  states.  This 
is  due  to  the  design  model's  representation  of  the  actuator 
states  as  independent  of  driving  noises  and  therefore 
having  response  that  depends  on  the  control  inputs  only 
with  no  uncertainty.  Thus,  as  determined  here,  estimation 
of  the  actuator  states  would  be  based  solely  on  the  filter's 
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internal  dynamics  model.  This  is  clearly  inappropriate 
since  the  actuator  dynamics  are  not  known  perfectly  and 
are  actually  subject  to  non-linearities  (rate-  and  position 
limiting) ,  external  dynamic  forces  (air  loads  affect 
achieved  surface  deflections  according  to  flight  condition) 
and  atmospheric  turbulence.  Although  it  is  desirable  that 
flight  controllers  not  include  control  surface  deflection 
angles  as  feedbacks,  it  is  clear  from  the  earlier  results 
of  Sections  6.4.6  and  6.4.7  that  the  sampled-data  imple¬ 
mentations  of  the  CGT/PI  controller  require  such  feedback. 
Therefore,  measurements  of  actuator  angles  would  be  neces¬ 
sary  and  the  design  model  should  include  artificial  noise 
affecting  the  actuator  dynamics  as  well  as  corruptive 
noise  on  the  measurements,  to  yield  nonzero  elements  where 
zeros  exist  in  equation  (6-70) . 

A  covariance  analysis  was  performed  for  this  filter 
design  with  respect  to  truth  model  AFTI ( S4 , A2 , G3 )  described 
in  Section  6.4.2  above.  The  "true"  and  filter-computed 
estimation  error  covariances  are  computed  for  fifty  filter 
sample  periods  by  the  method  described  in  Section  4.4. 

The  true  and  computed  standard  deviations  for  each  state 
estimate  during  the  fifty  samples  are  output  in  plot  form 
to  the  "LIST"  file.  In  addition,  the  standard  deviations 
at  the  final  sample  time  are  listed  directly  at  the  user's 
computer  terminal.  As  discussed  above,  since  the  actuator 
states  are  modeled  as  determinstic,  the  estimation  error 


is  computed  as  zero  for  all  time.  Figures  6-29  through 
6-34  are  plots  of  the  true  and  filter  computed  standard 
deviations  for  all  design  model  states  excepting  the 
actuator  states.  Plot  symbols  1  and  2  are  the  true  and 
filter-computed  standard  deviations,  respectively;  the 
horizontal  axis  gives  the  time  samples,  and  the  dependent 
axes  are  scaled  in  units  of  radians  or  radians  per  second, 
as  appropriate.  The  corresponding  final-sample  standard 
deviations  (printed  at  the  terminal)  are  (in  units  of 
radians  or  radians  per  second) 
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Fig.  6-29.  True(l)  and  Computed (2)  Standard  Deviations 
of  Estimation  Error  in  0 
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Pig.  6-30.  True(l)  and  Computed (2)  Standard  Deviations 
of  Estimation  Error  in  a 
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As  can  be  seen  in  equation  (6-71d) ,  the  estimation  error 

for  the  turbulence  model  state  a'  is  quite  large  (about 

16°)  while  that  for  state  a is  very  good  (about  0.1°) . 

In  the  turbulence  model  developed  in  Reference  21,  the 

angle-of-attack  gust  is  modeled  by  a  second-order  system. 

In  this  representation  (see  equations  (6-26a)  and  (6-26b) ) 

the  state  a  is  the  physical  variable  while  state  a'  is 
g  g 

simply  a  dummy  state  defined  arbitrarily  to  achieve  the 

desired  second-order  dynamics.  Also,  note  from  equation 

(6-62a)  that  the  feedback  gains  on  the  a'  state  are  smaller 

than  those  on  the  physical  state  by  a  factor  of  about 

0.003.  Thus,  the  apparently  large  estimation  error  for 

a'  is  not  a  problem  since  it  is  a  dummy  state  and  the  feed- 
y 

back  gain  on  it  in  the  controller  is  small.  If  desired, 
a  transformation  could  be  applied  to  the  turbulence  model 
so  as  to  achieve  comparable  estimation  errors  for  each  of 
the  two  new  angle-of-attack  gust  states. 

These  results  show  that  the  steady-state  Kalman 
filter  for  this  design  problem  achieves  state  estimation 
of  adequate  quality  for  the  controller  application.  As 
mentioned  above,  this  Kalman  filter  design  assumed  that 
control  surface  deflections  are  correctly  modeled  as  deter¬ 
ministic.  Successful  controller/filter  designs  would  need 
to  model  actuator  states  as  non-deterministic  and  would 
probably  require  that  measurements  of  the  actuator  deflec¬ 
tion  angles  be  taken,  and  made  available  to  the  Kalman 
filter.  Finally,  note  that  these  results  allow  evaluation 


only  of  the  filter  alone.  The  combined  controller /filter 
implemented  as  a  CGT/PI/KF  controller  must  be  evaluated 
as  an  ensemble  in  order  to  judge  the  adequacy  of  the  final 
and  complete  design.  Additional  design  iteration  of  the 
CGT/PI  and/or  Kalman  filter  elements  may  be  necessary  in 
light  of  such  an  evaluation.  Development  of  efficient 
software  for  evaluating  an  ensemble  CGT/PI/KF  is  a  strong 
recommendation  for  extending  the  work  of  this  thesis. 

This  would  allow  the  controller  to  be  evaluated  in  a 
realistic  implementation  and  to  explore  the  impact  on 
system  robustness  in  case  Kalman  filter  estimates  replace 
the  assumed  full-state  feedback  variables. 

6.4.9  Discussion  of  Results  for  FI ight  Control 
Design  Example .  The  designs  discussed  in  this  section 
have  demonstrated  the  applicability  of  the  CGT/PI  and 
Kalman  filter  design  and  performance  evaluation  techniques 
to  aircraft  flight  control  problems.  The  CGT/PI  designs 
achieve  both  conventional  and  decoupled  pitch  control. 
Command  models  providing  the  desired  closed-loop  dynamics 
are  readily  formulated.  The  controllers  achieve  good 
model -following  and,  in  the  case  of  perfect  state  knowledge, 
continue  to  give  good  model-following  performance  even  in 
the  face  of  parameter  error  in  the  design  model  dynamics 
description.  Finally,  the  sampled-data  implementation  of 
the  CGT/PI  controller  must  be  based  upon  a  design  model 
which  includes  models  of  control  actuator  lags  in  order  to 


ensure  stability  when  actuators  with  significant  lags 
actually  effect  control  surface  deflections. 

The  Kalman  filter  design  for  the  aircraft  longi¬ 
tudinal  dynamics  subject  to  clear-air  turbulence  is  readily 
achieved.  The  estimation  errors  determined  from  a  covari¬ 
ance  analysis  prove  to  be  sufficiently  small  to  be  used 
as  feedback  variables  in  the  final  closed-loop  controller. 
Because  of  the  possibility  of  large  gains  in  the  feedback 
and  feedforward  paths,  the  Kalman  filter  should  be  imple¬ 
mented  so  as  to  achieve  minimum  computational  delay.  Also 
it  is  noted  that  design  models  should  include  artificial 
noise  sources  driving  the  actuator  states,  and  that  measure¬ 
ments  of  the  actuator  deflections  would  be  needed. 

However,  the  evaluations  described  above  provide 
insight  into  the  CGT/PI  controller  and  Kalman  filter  per¬ 
formances  only  as  separate  entities.  In  particular,  the 
CGT/PI  evaluation  assumes  that  all  feedforward  and  feedback 
quantities  are  known  perfectly  and  with  no  delay.  The 
final  CGT/PI/KF  controller  requires  state  estimates  from 
the  Kalman  filter  in  order  to  generate  the  control  input 
commands.  In  this  case,  the  filter  would  itself  introduce 
additional  dynamics  and  the  state  estimates  are  available 
only  after  a  short  delay .  A  complete  evaluation  of  the 
full  CGT/PI/KF  controller  requires  additional  software 
implementing  the  controller/ filter  system  evaluation  compu¬ 
tations  as  described  in  Reference  32. 
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VII.  Conclusions  and  Recommendations 


for  Future  Research 


7.1  Conclusions 

The  objectives  of  this  study  have  been  twofold: 

(1)  develop  a  user-oriented  computer  program  aiding  in 
the  design  of  CGT/PI/KF  controllers  and  applicable  to  a 
wide  variety  of  specific  control  design  problems,  and 

(2)  employ  the  design  computer  program  to  develop  control 
designs  relevant  to  aircraft  flight  control  in  order  to 
evaluate  the  design  process  and  qualities  of  the  designs 
achieved.  These  objectives  have  been  accomplished  and  lead 
to  the  following  conclusions: 

1.  T.-j  design  computer  program  (CGTPIF)  developed 
during  this  study  provides  an  effective  means  to  pursue 
design  of  CGT/PI/KF  controllers. 

2.  Command  models  in  state  variable  form  may  be 
readily  formulated  for  desired  dynamic  behavior  described 
in  terms  of  transfer  functions  (e.g.,  direct  translation 
into  standard  controllable  form) .  The  command  models  may 
specify  either  conventional  or  decoupled  input-output 
behavior  for  the  closed-loop  system.  The  model  dynamics 
are  expressed  in  the  continuous-time  domain  and  discretized 
within  the  computer  program  to  determine  the  appropriate 
sampled-data  controller. 
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3.  CGT/PI  designs  achieving  good  model -following 
are  readily  achieved  in  a  systematic,  iterative  design  pro¬ 
cedure.  The  PI  regulator  is  first  designed  independently 
to  provide  rapid  and  well -damped  output  response  to  initial 
conditions.  This  design  is  then  incorporated  into  the 
closed-loop  CGT/PI  control  law  and  evaluated  in  terms  of 
the  model -following  accuracy  and  control  input  magnitudes 
and  rates  for  a  step  command  input.  Further  tuning  of  the 
PI  design  is  done  iteratively  (if  necessary)  to  achieve 
the  desired  CGT/PI  performance. 

4 .  The  CGT/PI/KF  controller  provides  a  design 
technique  ideally  suited  to  achieving  controllers  for 
sampled-data  applications.  The  controller  is  a  direct 
digital  design,  not  employing  approximations  to  transform 
a  continuous  design  to  discrete  implementation. 

5.  The  CGT/PI  controller  can  be  applied  success¬ 
fully  to  realistic  flight  control  design  problems.  The 
design  technique  is  particularly  useful  in  achieving  con¬ 
troller  designs  in  the  case  of  multiple  independent  control 
surfaces  and  multiple  outputs.  In  such  MIMO  design  problems, 
the  CGT/PI  controller  provides  all  single-  and  cross¬ 
channel  feedforward  and  feedback  gains  directly.  By 
employing  optimal  control  design  methods  based  upon  the 
Linear  system-Quadratic  cost  weighting  (LQ)  assumptions, 

a  stable  closed-loop  system  is  assured  and  the  single  and 
cross-channel  gains  are  determined  in  a  direct  and  logical 


manner. 


6.  In  the  case  of  perfect  state  knowledge,  the 


CGT/PI  controller  maintains  good  model -following  perform¬ 
ance  despite  disturbances  acting  on  system  dynamics,  and 
despite  parameter  errors  of  significant  magnitude  in  the 
design  model . 

7.  Design  of  the  corresponding  Kalman  filter  for 
the  controller  application  is  straightforward  and  presents 
no  unique  difficulties.  For  the  flight  control  designs 
considered,  a  covariance  analysis  employing  a  system  truth 
model  demonstrated  that  state  estimates  of  acceptable 
accuracy  for  use  in  the  CGT/PI  controller  are  obtained 
based  upon  a  reduced  design  model.  In  actual  implementa¬ 
tion,  actuator  states  would  need  to  be  measured  and  the 
design  model  should  include  pseudo-noise  driving  the  con¬ 
trol  surface  deflections  in  addition  to  the  deterministic 
control  actuator  dynamics. 

7.2  Recommendations  for  Future  Research 

As  a  result  of  experience  obtained  in  designing 
CGT/PI/KF  controllers  in  this  study,  several  areas  requir¬ 
ing  additional  study  are  apparent.  The  following  recommer 
dations  include  items  directed  at  extending  the  CGT/PI/KF 
theory,  modifying  and  extending  the  capabilities  of  the 
design  computer  program,  and  pursuing  further  evaluation 
of  CGT/PI/KF  designs. 


1.  Develop  software  extending  the  evaluation  capa¬ 
bility  of  the  existing  design  computer  program  to  encompass 
the  complete  closed-loop  CGT/PI/KF  controller.  Such  soft¬ 
ware  is  essential  for  evaluating  the  effects  due  to  Kalman 
filter  estimation  of  states  required  by  the  CGT/PI  and 
will  allow  final  tuning  of  the  entire  design.  In  par¬ 
ticular,  robustness  of  the  CGT/PI  and  CGT/PI/KF  designs 
can  be  conveniently  compared  with  such  software. 

2.  Modify  the  software  in  the  PI  regulator  design 
path  to  allow  quadratic  weights  to  be  applied  directly  to 
derivatives  of  the  output  variables;  these  would  be  manipu¬ 
lated  by  the  program  coding  to  achieve  appropriate  weights 
for  states,  inputs,  and  input  rates.  Improved  damping  of 
the  output  response  can  often  be  achieved  by  weighting 
both  output  and  output  rate  deviations  in  the  quadratic 
cost  function  of  the  optimal  control  formulation. 

3.  Develop  (and  implement)  in  the  existing  soft¬ 
ware)  an  extension  to  the  CGT  theory  to  provide  for 
matching  of  the  derivatives  of  the  design  and  command 
model  outputs  (as  well  as  meeting  objectives  as  discussed 
in  this  research) .  This  is  analogous  to  the  technique 
employed  in  the  implicit  model -following  method  and  will 
tend  to  improve  dynamics  matching  during  the  transient 
phase  of  response. 

4.  Using  the  software  developed  according  to 
Recommendation  1  above,  evaluate  the  flight  control  designs 
determined  in  this  study  as  complete  CGT/PI/KF  closed-loop 
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controllers.  Evaluation  should  focus  on  the  quality  of  the 
model-following  achieved  for  the  designs  in  cases  of 
nominal  truth  model  parameters  and  erroneous  truth  model 
parameters.  These  will  demonstrate  the  influence  of  filter 
dynamics  and  computational  delay  on  the  closed- loop  con¬ 
troller  performance  and  system  robustness  degradation  due 
to  filter  estimation  errors  caused  by  modeling  errors  and 
parameter  variation.  Robustness  recovery  techniques  (as 
described  in  Reference  32)  may  be  employed  in  the  final 
system  formulation,  if  considered  necessary. 

5.  Using  the  software  developed  according  to  Rec¬ 
ommendation  1  above,  compare  the  CGT/PI/KF  controller  per¬ 
formance  with  respect  to  disturbance  rejection  for  flight 
control  designs  based  upon  design  models  in  which  distur¬ 
bances  (gusts)  are  modeled  separately  or  included  among 
the  system  states.  The  former  case  entails  a  distinct 
dynamics  model  for  the  disturbance  states  and  the  resulting 
CGT/PI/KF  controller  would  employ  feedforward  gains  on  the 
estimated  disturbance  states  to  achieve  disturbance  rejec¬ 
tion  (CGT  disturbance  rejection) .  The  latter  case,  in 
which  the  disturbance  states  are  included  in  the  design 
model's  state  vector,  would  result  in  a  CGT/PI/KF  controller 
which  would  employ  feedback  gains  on  the  estimated  distur¬ 
bance  states  to  achieve  disturbance  rejection  (PI  distur¬ 
bance  rejection) .  In  either  case,  the  Kalman  filter  design 
is  invariant. 
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